Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,363 A→T 6.2% S100T (TCC→ACC)  cyoC ← cytochrome o ubiquinol oxidase subunit III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,3630AT 6.2% 366.3 / 8.9 127S100T (TCC→ACC) cyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base A (57/62);  new base T (5/3);  total (62/65)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 132 total aligned reads displayed.

GTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCA  >  NC_000913/448214‑448512
                                                                                                                                                     |                                                                                                                                                     
gTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGga                                                                                                                                                       <  1:899851/150‑1 (MQ=255)
    aGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATa                                                                                                                                                   <  1:194440/150‑1 (MQ=255)
       ccGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAAcc                                                                                                                                                >  1:13538/1‑150 (MQ=255)
        cgtgcgtgCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCATGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCt                                                                                                                                               <  1:699785/150‑1 (MQ=255)
           gcgtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGc                                                                                                                                            <  1:1059837/150‑1 (MQ=255)
           gcgtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTAAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGTTATCTCCTGGc                                                                                                                                            <  1:116803/150‑1 (MQ=255)
            cgtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCt                                                                                                                                           >  2:1233655/1‑150 (MQ=255)
             gtgcCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCAGATAAATCCGGCACCAAACAACCAGGTCAACTCCAGCCACGAGATAACCTCGCtt                                                                                                                                          >  1:767564/1‑150 (MQ=255)
              ggcCGACCAGCGCAAAGAACGCTGCCAGGAAGCCTCTGCGATCCTGACCCATGCCGTTAACAAGCAGGTGATGGAATTCATAGATTTCCAGCCCGAGAAATCCGGCACCATCCAACCAGGTCACCGCCAGCCAGGAGATAACCTGGcggg                                                                                                                                         <  1:1203192/149‑4 (MQ=255)
                ccGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATAAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGTTATCTCCTGGCtggcg                                                                                                                                       <  2:722750/150‑5 (MQ=255)
                   aCCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgttg                                                                                                                                    <  2:928133/150‑1 (MQ=255)
                   aCCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTttgttg                                                                                                                                    <  1:1378198/150‑1 (MQ=255)
                   aCCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTCATGGAATTCATAGATATCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCATGAGATAACCTGGCTTttgttg                                                                                                                                    <  2:1040990/150‑1 (MQ=255)
                    ccAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTtgttgt                                                                                                                                   >  2:1466570/1‑150 (MQ=255)
                       gcgcAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGtttt                                                                                                                                <  1:1685520/150‑1 (MQ=255)
                        cgcAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGttttt                                                                                                                               <  1:796673/150‑1 (MQ=255)
                          cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGt                                                                                                                             <  1:775555/150‑1 (MQ=255)
                          cAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTa                                                                                                                            >  2:1462503/1‑150 (MQ=255)
                          cAAAGAACGCGGACAGGAAGCCGCGGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGt                                                                                                                             <  2:875980/150‑1 (MQ=255)
                           aaaGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTa                                                                                                                            <  1:308092/150‑1 (MQ=255)
                           aaaGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTa                                                                                                                            <  1:747975/150‑1 (MQ=255)
                               aaCGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCACCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATg                                                                                                                        >  2:199806/1‑150 (MQ=255)
                                aCGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATgg                                                                                                                       <  2:981634/150‑1 (MQ=255)
                                aCGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATgg                                                                                                                       <  2:1535338/150‑1 (MQ=255)
                                 cGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGc                                                                                                                      <  1:12155/150‑1 (MQ=255)
                                    tGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGAt                                                                                                                   >  1:450279/1‑150 (MQ=255)
                                    tGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGAt                                                                                                                   >  1:724120/1‑150 (MQ=255)
                                      aCAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTACATACCGATAAATACGGCACCAACCAACCAGGTCAACGCCAGCCAGGAGAGAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAg                                                                                                                 >  2:257194/1‑150 (MQ=255)
                                        aGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAgcc                                                                                                               <  1:380243/150‑1 (MQ=255)
                                            aGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCa                                                                                                           >  2:590360/1‑150 (MQ=255)
                                                gCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTAAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGcc                                                                                                       <  2:1565405/150‑1 (MQ=255)
                                                 cTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCg                                                                                                      >  1:1024518/1‑150 (MQ=255)
                                                    cGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAg                                                                                                   >  2:1421694/1‑150 (MQ=255)
                                                    cGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCAACCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAg                                                                                                   >  2:27033/1‑150 (MQ=255)
                                                     gATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAgg                                                                                                  >  1:302498/1‑150 (MQ=255)
                                                      aTCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGt                                                                                                 >  2:1439791/1‑150 (MQ=255)
                                                      aTCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGt                                                                                                 >  2:1377418/1‑150 (MQ=255)
                                                         cGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGAt                                                                                              >  1:1019595/1‑150 (MQ=255)
                                                           gACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATgg                                                                                            <  1:585265/150‑1 (MQ=255)
                                                           gACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATgg                                                                                            <  2:1292699/150‑1 (MQ=255)
                                                             cccATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAg                                                                                          >  1:1593593/1‑150 (MQ=255)
                                                                     gTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaacaa                                                                                  >  2:572262/1‑150 (MQ=255)
                                                                     gTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGGTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGaacaa                                                                                  >  1:1520108/1‑150 (MQ=255)
                                                                            aTCAGGTGATGGAGTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAg                                                                           <  2:1663156/150‑1 (MQ=255)
                                                                            aTCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCa                                                                             >  2:962284/1‑148 (MQ=255)
                                                                            aTCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCa                                                                             <  1:962284/148‑1 (MQ=255)
                                                                               aGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGaaa                                                                        >  1:987120/1‑150 (MQ=255)
                                                                                  tGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGtt                                                                     >  2:705957/1‑150 (MQ=255)
                                                                                     tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCt                                                                                        <  1:1022225/128‑1 (MQ=255)
                                                                                     tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCt                                                                                        >  2:1022225/1‑128 (MQ=255)
                                                                                     tGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCa                                                                  >  2:1551415/1‑150 (MQ=255)
                                                                                        aaTTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAAcc                                                               <  2:1402948/150‑1 (MQ=255)
                                                                                         aTTCATAGAATTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACATGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGACCAACAGCAAGAAAGTTTCAACCa                                                              >  1:307781/1‑150 (MQ=255)
                                                                                          ttCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAg                                                             <  1:845084/150‑1 (MQ=255)
                                                                                           tCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAga                                                            <  2:1661842/150‑1 (MQ=255)
                                                                                            cATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaa                                                           >  1:1154381/1‑150 (MQ=255)
                                                                                             aTAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaac                                                          <  2:525938/150‑1 (MQ=255)
                                                                                              tAGATTTCCATCCCGATAAATCCGGCAACAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaacg                                                         <  1:1439791/150‑1 (MQ=255)
                                                                                                gATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaacgaa                                                       >  1:957165/1‑150 (MQ=255)
                                                                                                gATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTCGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaacgaa                                                       >  2:278188/1‑150 (MQ=255)
                                                                                                 aTTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAgaacgaac                                                      >  1:452354/1‑150 (MQ=255)
                                                                                                         cccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTc                                              <  1:1233655/150‑1 (MQ=255)
                                                                                                          ccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCg                                             >  1:735009/1‑150 (MQ=255)
                                                                                                          ccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCg                                             >  2:1438495/1‑150 (MQ=255)
                                                                                                          ccGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCg                                             >  1:181795/1‑150 (MQ=255)
                                                                                                           ctaaaaaaCCGCCACCAAACACCCAGGTCAACGCCAGCCAGGACATACCCTGGCTTTTGTTTTTTTTGTCCATGGCGAACGCCGCCATGCCTTAGGTGATGGCGCTCAACAACACCAAGAAAGTTCCACCCACAACTAACGGCAGTTCGa                                            <  2:1382019/146‑1 (MQ=255)
                                                                                                           cGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGa                                            <  2:519614/150‑1 (MQ=255)
                                                                                                            gATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaa                                           >  1:1656117/1‑150 (MQ=255)
                                                                                                             aTAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGaaa                                          <  2:1593593/150‑1 (MQ=255)
                                                                                                               aaaTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACTGCCGTTCgttctg                                        <  1:606360/150‑6 (MQ=255)
                                                                                                                aaTCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATg                                       <  2:763083/150‑1 (MQ=255)
                                                                                                                   ccGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTcc                                    <  1:590360/150‑1 (MQ=255)
                                                                                                                       cACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAc                                <  2:724120/150‑1 (MQ=255)
                                                                                                                       cACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAc                                <  2:450279/150‑1 (MQ=255)
                                                                                                                       cACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAc                                <  2:987120/150‑1 (MQ=255)
                                                                                                                        aCCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAcc                               <  2:35947/150‑1 (MQ=255)
                                                                                                                        aCCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTAcc                               <  1:705957/150‑1 (MQ=255)
                                                                                                                          cAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGCTAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTg                             >  1:256137/1‑150 (MQ=255)
                                                                                                                           aaaCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGt                            <  2:1464435/150‑1 (MQ=255)
                                                                                                                           aaaCAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGCAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGt                            >  2:1108964/1‑150 (MQ=255)
                                                                                                                               aacCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCggg                        <  2:840917/150‑1 (MQ=255)
                                                                                                                                acCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGgc                       <  1:1462503/150‑1 (MQ=255)
                                                                                                                                  cAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGATCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGgccg                     <  2:1098850/150‑1 (MQ=255)
                                                                                                                                  cAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGgccg                     >  2:413319/1‑150 (MQ=255)
                                                                                                                                  cAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGATCGAAAATGTCCTTACCTCTCGGgccg                     >  2:1044389/1‑150 (MQ=255)
                                                                                                                                        aaCGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCg               >  1:700681/1‑150 (MQ=255)
                                                                                                                                          caaaagccagGTTATCTCCTGGCTGGCGTTGTTTTTGTACATGGCGATAGACGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGt             >  2:498658/5‑150 (MQ=255)
                                                                                                                                          cGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGt             <  2:302498/150‑1 (MQ=255)
                                                                                                                                            aaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGc           >  1:539877/3‑150 (MQ=255)
                                                                                                                                            aaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGc           >  1:1237896/3‑150 (MQ=255)
                                                                                                                                            aaagccagGTTATCTCCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGcc          >  1:1569874/3‑150 (MQ=255)
                                                                                                                                             cagccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGcc          <  2:1213185/150‑1 (MQ=255)
                                                                                                                                              agccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCg         >  1:1586110/1‑150 (MQ=255)
                                                                                                                                              agccagGAGATAACCTGGCTTTTGTGGTGTTTGTACATGGCTATAGCCGCCATGCCGTAGGTGAGGGAGCTGAACAACAGCAAGAAAGTATCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCg         <  2:445115/150‑1 (MQ=255)
                                                                                                                                               gccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGt        >  2:1185733/1‑150 (MQ=255)
                                                                                                                                                ccagGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCACCTGCGGTGCCGtt       <  2:44074/150‑1 (MQ=255)
                                                                                                                                                  agGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCa     <  2:1024518/150‑1 (MQ=255)
                                                                                                                                                    gagaTAACCTGGCTTTTGTTGTTATTGTACAGGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTGCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCAcc   <  1:840952/150‑1 (MQ=255)
                                                                                                                                                     agaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCa  >  1:1476513/1‑150 (MQ=255)
                                                                                                                                                     agaTAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCa  >  1:1670234/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCAGACCGTGCGTGCCGACCAGCGCAAAGAACGCTGACAGGAAGCCGCTGCGATCCGGACCCATGCCGTTAACAATCAGGTGATGGAATTCATAGATTTCCATCCCGATAAATCCGGCACCAAACAACCAGGTCAACGCCAGCCAGGAGATAACCTGGCTTTTGTTGTTTTTGTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCA  >  NC_000913/448214‑448512

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