Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,521,054 A→T 6.2% V671D (GTC→GAC)  yncD ← putative iron outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,521,0540AT 6.2% 172.4 / 2.1 64V671D (GTC→GAC) yncDputative iron outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base T (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCA  >  NC_000913/1520910‑1521198
                                                                                                                                                |                                                                                                                                                
cAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAg                                                                                                                                               >  3:16867/1‑148 (MQ=255)
 aGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGaa                                                                                                                                             <  7:5959/149‑1 (MQ=255)
      tACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAAc                                                                                                                                        >  3:166918/1‑149 (MQ=255)
       aCAATCGTGAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCAGGCCGCCACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGCCTCATTGACAATGACAGAACCCACg                                                                                                                                       >  5:258500/1‑149 (MQ=255)
         aaTCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTa                                                                                                                                     >  6:36923/1‑149 (MQ=255)
                aaaGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCttt                                                                                                                               <  7:102170/148‑1 (MQ=255)
                 accttgaTGCCCGGTATTTGCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACCCCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTAt                                                                                                                             >  5:96123/7‑149 (MQ=255)
                 aaGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTAt                                                                                                                             >  6:345264/1‑149 (MQ=255)
                  aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCCCGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATc                                                                                                                            <  2:164566/149‑1 (MQ=255)
                  aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATc                                                                                                                            <  2:159155/149‑1 (MQ=255)
                  aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATc                                                                                                                            <  4:166918/149‑1 (MQ=255)
                  aGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATc                                                                                                                            <  1:75136/149‑1 (MQ=255)
                    tAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCg                                                                                                                           >  5:13656/1‑148 (MQ=255)
                     aaaTGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGCCACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGaa                                                                                                                         >  2:112092/1‑149 (MQ=255)
                      aaTGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGaa                                                                                                                         <  5:345264/148‑1 (MQ=255)
                           ccGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAAt                                                                                                                    >  8:339053/1‑148 (MQ=255)
                           ccGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAAt                                                                                                                    >  1:130760/1‑148 (MQ=255)
                           ccGGCAAATACAGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAAtt                                                                                                                   >  8:365461/1‑149 (MQ=255)
                                    aCCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGacac                                                                                                           >  7:290438/1‑148 (MQ=255)
                                            tCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCaaat                                                                                                   >  1:394286/1‑148 (MQ=255)
                                              aaGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCaaataaa                                                                                                >  3:131073/1‑149 (MQ=255)
                                                    aTTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCca                                                                                            >  1:156227/1‑147 (MQ=255)
                                                         tCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGTCAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTa                                                                                      >  8:8603/1‑148 (MQ=255)
                                                         tCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGTCAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTa                                                                                      >  8:14281/1‑148 (MQ=255)
                                                                  cACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATtgtg                                                                               <  1:17376/146‑1 (MQ=255)
                                                                  cACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATtgtg                                                                               <  5:36923/146‑1 (MQ=255)
                                                                  cACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTgg                                                                              <  4:81638/147‑1 (MQ=255)
                                                                           ttCATGCCGACACCCTAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTAtt                                                                    <  2:160446/148‑1 (MQ=255)
                                                                            tCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGCCCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTta                                                                  >  6:86113/1‑149 (MQ=255)
                                                                             cATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTta                                                                  <  6:96123/148‑1 (MQ=255)
                                                                             cATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTta                                                                  <  4:32866/148‑1 (MQ=255)
                                                                              aTGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTtata                                                                >  5:142223/1‑149 (MQ=255)
                                                                              aTGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTtata                                                                >  6:369093/1‑149 (MQ=255)
                                                                                   gACAACATAATTTCGTCCGGGCGAAGGTTCGTAATATAGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCgg                                                           <  2:319580/149‑1 (MQ=255)
                                                                                    ccacCATAATTTCGTCCGGGCGAAGGGTCGTAATAGCGCCCGTTGGACTCAGTGACAAGGACAGAACCACCGTAGTCTTTCTCGAATCAATTATCGACACGACCAAATAACTCCACAGTTAACTTGTGGGAATTGTATTTATACCCGGt                                                          <  6:258500/148‑1 (MQ=255)
                                                                                    acacCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGt                                                          >  2:254322/1‑149 (MQ=255)
                                                                                       ccATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTaaa                                                                                    >  4:312462/1‑120 (MQ=255)
                                                                                       ccATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTaaa                                                                                    <  3:312462/120‑1 (MQ=255)
                                                                                        cATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGaa                                                       >  8:358608/1‑148 (MQ=255)
                                                                                        cATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGaa                                                       >  6:224055/1‑148 (MQ=255)
                                                                                          tAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATaa                                                    <  4:131073/149‑1 (MQ=255)
                                                                                            aTTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGc                                                  <  2:130760/149‑1 (MQ=255)
                                                                                             tttCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATAGACACGACCAAATAAATCCACAGTTAAATTGTTGTAAtt                                                                         >  1:71195/1‑125 (MQ=255)
                                                                                                  tCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCcgacga                                            <  7:358608/149‑1 (MQ=255)
                                                                                                    cGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACgaga                                          <  1:254322/149‑1 (MQ=255)
                                                                                                       gCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCGTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACACGGTGAATAAGCCGACGAGAGt                                        >  7:272726/1‑148 (MQ=255)
                                                                                                            ggTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGa                                  <  8:290438/149‑1 (MQ=255)
                                                                                                              ttCGTAATATCGCCCGTTTCACTCATTGACAATGAAAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGATAAATTGTGGTAATTGTAATTAAACCCGCTGAATAAGCAGACCAGAGTATAAGACg                                >  2:274631/1‑149 (MQ=255)
                                                                                                                cGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGgc                              <  3:264944/149‑1 (MQ=255)
                                                                                                                     tatCGCCCGTTTGACTCATTGAGAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATAGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCttt                         <  1:274631/149‑1 (MQ=255)
                                                                                                                     tatCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCttt                         <  2:394286/149‑1 (MQ=255)
                                                                                                                       tCGACCTTTTGACTCATTGACAATTATAGAACCAACGTATTCTTTATCGAATAATTTATCGAGACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCtttt                        <  6:142223/148‑1 (MQ=255)
                                                                                                                         gCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTCTCGACACGACCAAATAAATCGACCGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCGGACGAGAGTATAAGACGGCGCTTTTGcc                     <  2:217/149‑1 (MQ=255)
                                                                                                                             gTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTta                                                                  >  4:195175/1‑100 (MQ=255)
                                                                                                                             gTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTta                                                                  <  3:195175/100‑1 (MQ=255)
                                                                                                                                  aCTCATTGACAATGTCAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTtca            <  7:14281/149‑1 (MQ=255)
                                                                                                                                  aCTCATTGACAATGTCAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTtca            <  7:8603/149‑1 (MQ=255)
                                                                                                                                  aCTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTtc             <  7:339053/148‑1 (MQ=255)
                                                                                                                                  aCTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAAATGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTtc             <  7:365461/148‑1 (MQ=255)
                                                                                                                                         gACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGCCACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTg     >  1:190469/1‑149 (MQ=255)
                                                                                                                                            aaTGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCa  >  5:85796/1‑149 (MQ=255)
                                                                                                                                            aaTGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCa  <  5:224055/149‑1 (MQ=255)
                                                                                                                                               gACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGt                                                          <  2:292072/90‑1 (MQ=255)
                                                                                                                                               gACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGt                                                          >  1:292072/1‑90 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGCAATACAATCGTTAAAGTAAATGCCCGGCAAATACCGGGCATCAAGGTAATTACTCAAATCTCCACGCAATATTCATGCCGACACCATAATTTCGTCCGGGCGAAGGTTCGTAATATCGCCCGTTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCA  >  NC_000913/1520910‑1521198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑