Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 420,945 A→C 5.2% intergenic (+35/‑41) brnQ → / → proY branched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val/proline‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913420,9450AC 5.2% 678.4 / 23.0 193intergenic (+35/‑41)brnQ/proYbranched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val/proline‑specific permease
Reads supporting (aligned to +/- strand):  ref base A (85/98);  new base C (4/6);  total (89/104)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 195 total aligned reads displayed.

TTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCA  >  NC_000913/420796‑421094
                                                                                                                                                     |                                                                                                                                                     
ttaaCCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACAGTTGTTTTAACCACGGGGCTGCGa                                                                                                                                                       >  2:2154207/4‑150 (MQ=255)
  tACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATg                                                                                                                                                     >  1:1294309/1‑150 (MQ=255)
    ccGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGcc                                                                                                                                                   <  1:1560910/150‑1 (MQ=255)
       cTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCg                                                                                                                                                <  2:1150554/150‑1 (MQ=255)
       cTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACAGTTGTTTTAACCACGGGGCTGCGATGCCCCg                                                                                                                                                >  2:237641/1‑150 (MQ=255)
        tGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGt                                                                                                                                               <  2:1001486/150‑1 (MQ=255)
        tGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGt                                                                                                                                               <  2:2993825/150‑1 (MQ=255)
         ggCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTg                                                                                                                                              <  1:1179932/150‑1 (MQ=255)
         ggCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTg                                                                                                                                              <  1:1199077/150‑1 (MQ=255)
         ggCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTg                                                                                                                                              <  1:2102095/150‑1 (MQ=255)
          gCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTgg                                                                                                                                             >  1:3205881/1‑150 (MQ=255)
          gCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTgg                                                                                                                                             >  1:412470/1‑150 (MQ=255)
          gCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTgg                                                                                                                                             <  2:1808132/150‑1 (MQ=255)
          gCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCATCGCAGCCCCGTgg                                                                                                                                             <  2:1568225/150‑1 (MQ=255)
            cGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGtt                                                                                                                                           >  2:721428/1‑150 (MQ=255)
              aaCAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGtttt                                                                                                                                         <  1:1257112/150‑1 (MQ=255)
                cAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGtttttt                                                                                                                                       <  2:1572783/150‑1 (MQ=255)
                 aaGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTa                                                                                                                                      >  1:69144/1‑150 (MQ=255)
                  aGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAt                                                                                                                                     <  2:1292771/150‑1 (MQ=255)
                      cTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATtgtg                                                                                                                                 <  2:2263404/150‑1 (MQ=255)
                      cTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATtgtg                                                                                                                                 <  2:1197353/150‑1 (MQ=255)
                            tGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATg                                                                                                                           >  1:58399/1‑150 (MQ=255)
                              gTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATggg                                                                                                                         >  1:1094762/1‑150 (MQ=255)
                               ttAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGt                                                                                                                        <  2:2091565/150‑1 (MQ=255)
                                 aaTGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTa                                                                                                                      >  1:2543683/1‑150 (MQ=255)
                                 aaTGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTa                                                                                                                      >  1:2242091/1‑150 (MQ=255)
                                 aaTGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTa                                                                                                                      >  1:1622187/1‑150 (MQ=255)
                                     ccAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGaa                                                                                                                  >  1:2130279/1‑150 (MQ=255)
                                          aGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGAt                                                                                                             <  2:2817601/150‑1 (MQ=255)
                                           gtggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGACTTGATg                                                                                                            >  1:2770633/1‑150 (MQ=255)
                                             ggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGa                                                                                                          <  1:1921932/150‑1 (MQ=255)
                                              gtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaa                                                                                                         <  2:3036304/150‑1 (MQ=255)
                                              gtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGaa                                                                                                         <  2:2823235/150‑1 (MQ=255)
                                                 aTGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGt                                                                                                      >  2:3128393/1‑150 (MQ=255)
                                                    gtggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAACTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAg                                                                                                   <  1:2888753/150‑1 (MQ=255)
                                                    gtggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAg                                                                                                   >  2:1336451/1‑150 (MQ=255)
                                                    gtggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAg                                                                                                   <  1:450969/150‑1 (MQ=255)
                                                    gtggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAg                                                                                                   <  1:450788/150‑1 (MQ=255)
                                                     tggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGa                                                                                                  <  1:916842/150‑1 (MQ=255)
                                                        ttCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACa                                                                                               >  1:1297297/1‑150 (MQ=255)
                                                          cTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAg                                                                                             <  1:2420053/150‑1 (MQ=255)
                                                           tGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGc                                                                                            >  1:2345941/1‑150 (MQ=255)
                                                           tGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGc                                                                                            >  2:54669/1‑150 (MQ=255)
                                                            ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTGTTTATTGTGTTGATGGGTTAGGACTTGATGGAAAGTAAGAACAAGCt                                                                                           >  1:1583266/1‑150 (MQ=255)
                                                             gCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTa                                                                                          <  1:2193164/150‑1 (MQ=255)
                                                             gCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTa                                                                                          >  1:1706195/1‑150 (MQ=255)
                                                              ccATTATCTGTGATCGCGCGGCAGCTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTTAACATTTGTTTTAACAACGTGGCTGAGATGACCCGTGGATTTTTATTGTGTTTATGGGGTAGGAATTGATGGAAAGTAAGAAAAAGCTaa                                                                                         >  2:1943701/1‑150 (MQ=255)
                                                                 ttATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGc                                                                                      >  2:605055/1‑150 (MQ=255)
                                                                          aTCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTaa                                                                             <  2:2718651/150‑1 (MQ=255)
                                                                          aTCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTaa                                                                             <  2:1417258/150‑1 (MQ=255)
                                                                               gcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTATATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTAcc                                                                        >  2:1072217/1‑150 (MQ=255)
                                                                                         gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATAc                                                              <  1:1056701/150‑1 (MQ=255)
                                                                                         gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATAc                                                              >  2:2317995/1‑150 (MQ=255)
                                                                                         gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATAc                                                              >  1:2272368/1‑150 (MQ=255)
                                                                                           cAGGTGACCTCCATAGATCTATAAATCACTGAACTTTTGTTTTAACCACGGGGCTGCGATGCCCGGTGTTTTTTTATTGTGTTGATGGGTTAGGATTTGATGTAAAGTAAGAACAAGCTAAAGCGTGGTCTAAGTACCCGCCACATACGc                                                            <  1:721428/150‑1 (MQ=255)
                                                                                               tGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTa                                                        <  1:843218/150‑1 (MQ=255)
                                                                                               tGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTa                                                        <  1:2073490/150‑1 (MQ=255)
                                                                                               tGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTa                                                        <  1:2580327/150‑1 (MQ=255)
                                                                                               tGACCTCCAGCGCTCACTAAATCACTGAACAGTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTGTTTATTGGGTTGCTGGGTTAGGAAGTGAGGGACAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCtttc                                                        >  2:2224616/1‑149 (MQ=255)
                                                                                                 aCCTCCAGCGCTCACTAACTCACTGAACATTTGTGTTAACCACGGGGCTGCGATGCCCCGTGGATTTTTCTTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTATGTACCCGCCACATACGCTTTATg                                                      <  1:2213206/150‑1 (MQ=255)
                                                                                                 aCCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATg                                                      <  1:1408293/150‑1 (MQ=255)
                                                                                                  ccTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATgg                                                     >  2:1627636/1‑150 (MQ=255)
                                                                                                    tCCAGCGCTCACTAAATCACTGAACATTTGTTTTCACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCa                                                   <  2:2875782/150‑1 (MQ=255)
                                                                                                    tCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCGTTATGGCa                                                   >  1:477689/1‑150 (MQ=255)
                                                                                                     ccAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCAc                                                  <  2:2625463/150‑1 (MQ=255)
                                                                                                        gcgcTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTgg                                               <  1:2543300/150‑1 (MQ=255)
                                                                                                           cTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGtt                                            >  2:2420451/1‑150 (MQ=255)
                                                                                                            tCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTc                                           >  2:3129319/1‑150 (MQ=255)
                                                                                                             cACTAACTCACTGCACAGTTGTTTTAACCCCGGGGCTGCGATGCCCCGTGGTTTTTTATTGGGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACACGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                          <  2:57751/150‑1 (MQ=255)
                                                                                                             cACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                          <  2:2542784/150‑1 (MQ=255)
                                                                                                             cACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                          <  1:2979794/150‑1 (MQ=255)
                                                                                                             cACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCa                                          <  1:906760/150‑1 (MQ=255)
                                                                                                              aCTAAATCACTGAACCTTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGAGGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTGTATGGCACTGGGTTCAg                                         >  2:731187/1‑150 (MQ=255)
                                                                                                                   aTCCCTGAACATTTGTTTTAACCTCGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAGCAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAAtt                                    <  2:2736501/150‑1 (MQ=255)
                                                                                                                   aTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAGCAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAAtt                                    <  2:2883105/150‑1 (MQ=255)
                                                                                                                   aTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAAtt                                    <  2:2943286/150‑1 (MQ=255)
                                                                                                                    tCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTg                                   <  2:2066380/150‑1 (MQ=255)
                                                                                                                     cACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTgg                                  <  1:209777/150‑1 (MQ=255)
                                                                                                                       cTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCa                                <  1:1919783/150‑1 (MQ=255)
                                                                                                                       cTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTGTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCAAAAGCGTGGGCTAAGTACCCGCCCCAGACGCGTTATGGCACAGGGTTCAGCACATGGCa                                >  1:2517968/1‑150 (MQ=255)
                                                                                                                         gAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGTTGGGTTTGGAATTGTTGGAAAGTAAGAACAAGCTAAATCGTGGGCGAAGTACCCGCCACAAACCCGGTATGGCACTGGGTTCAGCAAATGGCAcc                              >  1:553206/1‑150 (MQ=255)
                                                                                                                             aattGTTTAAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGGGTTGATGGGATAGGAATTGATGGACAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCACTTGGCACCGGGc                          >  2:2425702/3‑150 (MQ=255)
                                                                                                                              tttGTTTTAACCACGGGGCTTCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCt                         <  2:1147315/150‑1 (MQ=255)
                                                                                                                              tttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGTCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCt                         <  2:1647438/150‑1 (MQ=255)
                                                                                                                              tttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCt                         >  1:745091/1‑150 (MQ=255)
                                                                                                                              tttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCt                         >  1:1822975/1‑150 (MQ=255)
                                                                                                                              tgtGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCt                         <  2:2214742/148‑1 (MQ=255)
                                                                                                                               ttGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTTTGTTGTTTGGTTAGGAATTTATGGAAAGTAAGAAAAAGCTAAAACCTGGGCTAAGTAACCCCCACAAAACCGTTATTGCACAGGGTTAAGCAATTGGCACCGGGCTg                        >  1:1479479/1‑150 (MQ=255)
                                                                                                                               aataaaaaaCCACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTg                        >  2:2788806/8‑150 (MQ=255)
                                                                                                                                    ttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTa                   <  2:616233/150‑1 (MQ=255)
                                                                                                                                    ttAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTa                   <  2:907514/150‑1 (MQ=255)
                                                                                                                                     aaaCCACGGGGCATCGCAGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGTTGGTAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACAAACGCATTATGGCACGGGGTTCATCAATTGGCACCGGGCTGTTTTAc                  >  1:2843081/2‑150 (MQ=255)
                                                                                                                                      aaCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACg                 <  1:1336451/150‑1 (MQ=255)
                                                                                                                                      aaCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACg                 <  1:2004825/150‑1 (MQ=255)
                                                                                                                                         cACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATGGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGtt              >  2:1862602/1‑150 (MQ=255)
                                                                                                                                          aCGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATTGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTAACCGCCACATACGCTTTATGGCACGGGGTTCAGCAATTTGCACCGGGCTGTTTTACGGTTc             >  1:3111819/1‑150 (MQ=255)
                                                                                                                                          aCGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCGTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTc             >  1:2016834/1‑150 (MQ=255)
                                                                                                                                                 tGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGAc      >  2:2351582/1‑150 (MQ=255)
                                                                                                                                                  gCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACg     <  2:1294309/150‑1 (MQ=255)
                                                                                                                                                    gATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGcc   <  2:474106/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCA  >  NC_000913/420796‑421094

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: