Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,802,411 C→A 5.6% A355E (GCA→GAA)  nrdE → ribonucleoside‑diphosphate reductase 2, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,802,4110CA 5.6% 316.6 / 6.3 108A355E (GCA→GAA) nrdEribonucleoside‑diphosphate reductase 2, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (60/41);  new base A (3/3);  total (63/45)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 108 total aligned reads displayed.

ATGGCGCTGTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCT  >  NC_000913/2802263‑2802558
                                                                                                                                                    |                                                                                                                                                   
aTGGCGCTGTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCa                                                                                                                                                    <  2:145778/150‑1 (MQ=255)
   gcgcTGTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGaa                                                                                                                                                 >  2:780850/1‑150 (MQ=255)
     gcTGTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTAGCCGATGTGGCCATCAGCCAACACTATGACGATCTGGTTGACGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAAt                                                                                                                                               >  2:205460/1‑150 (MQ=255)
      cTGTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATc                                                                                                                                              >  1:1047025/1‑150 (MQ=255)
       tGTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGACATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACATCAACGCCCGTGATTTCTTCCAGCGAATGGCAGAAATcc                                                                                                                                             >  1:1306723/1‑150 (MQ=255)
       tGTTTTAGGCTTCTGACGTAGAGAGAGTTTATGGGAAGACGTTTGCCGATGGGGCCATCAGCCAACACTATGACGAACTGGGTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATcc                                                                                                                                             <  2:312065/150‑1 (MQ=255)
        gTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTGCCAGTCGCTGGAAGAAATcac                                                                                                                                            <  1:800715/150‑3 (MQ=255)
           ttCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTGCCAGTCGCTGGAAGAAATcacggg                                                                                                                                         <  1:327677/150‑6 (MQ=255)
            tCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGttt                                                                                                                                        <  1:523428/150‑1 (MQ=255)
             cGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTg                                                                                                                                       <  1:547919/150‑1 (MQ=255)
              gCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCTATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGa                                                                                                                                      <  1:883250/150‑1 (MQ=255)
              gCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAACTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGa                                                                                                                                      >  1:292547/1‑150 (MQ=255)
                 ttATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGACCTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCAAGTTTGAGTc                                                                                                                                   >  1:1061851/1‑150 (MQ=255)
                 ttATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATACAGTTTGAGTc                                                                                                                                   >  1:1056498/1‑150 (MQ=255)
                  tATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTcc                                                                                                                                  >  1:442086/1‑150 (MQ=255)
                     gACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGc                                                                                                                               <  1:605884/150‑1 (MQ=255)
                       cGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTa                                                                                                                             >  1:1373111/1‑150 (MQ=255)
                       cGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTa                                                                                                                             >  1:141143/1‑150 (MQ=255)
                        gTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTAt                                                                                                                            >  1:905583/1‑150 (MQ=255)
                            agCGAGTTTTTGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCGGTGATTTCTTCCAGCGTCTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCt                                                                                                                        <  2:267122/150‑1 (MQ=255)
                             gCGAGTTTATGGCTAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTa                                                                                                                       <  1:193602/150‑1 (MQ=255)
                                 gTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatc                                                                                                                   >  1:1213368/1‑150 (MQ=255)
                                   ttATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatcat                                                                                                                 >  2:1027817/1‑150 (MQ=255)
                                   ttATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTAcatcat                                                                                                                 >  1:861874/1‑150 (MQ=255)
                                     aTGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGt                                                                                                               <  1:1228675/150‑1 (MQ=255)
                                        gCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCACCGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATg                                                                                                            >  1:931095/1‑150 (MQ=255)
                                         cAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGa                                                                                                           >  2:339414/1‑150 (MQ=255)
                                         cAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGa                                                                                                           >  2:1587193/1‑150 (MQ=255)
                                         cAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCAGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGa                                                                                                           >  1:1341391/1‑150 (MQ=255)
                                          aaGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGaa                                                                                                          >  1:1571297/1‑150 (MQ=255)
                                          aaGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGaa                                                                                                          <  1:12494/150‑1 (MQ=255)
                                          aaGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGaa                                                                                                          <  1:1324850/150‑1 (MQ=255)
                                           aGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAg                                                                                                         <  2:323064/150‑1 (MQ=255)
                                               gTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGacac                                                                                                     >  2:1692914/1‑150 (MQ=255)
                                                tttGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACg                                                                                                    <  1:1478395/150‑1 (MQ=255)
                                                   gCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTa                                                                                                 >  1:1391643/1‑150 (MQ=255)
                                                      gATGTGGCAATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATAACTCAACGCCCGTGATTTCTTCATGCGACTGGCAGAAATACATTTTGAGTCCGGGTATCACTTCATCATTTATTAAGAAACTGTaaag                                                                                              >  2:149541/1‑149 (MQ=255)
                                                          tgGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTg                                                                                          >  1:1442304/1‑150 (MQ=255)
                                                                tCAGCCAACAATATGACGACCTGGTTGCCGATGAACGCATGCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTGTGAGTCCGGCTCTCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAAcc                                                                                    <  1:867845/150‑1 (MQ=255)
                                                                  aGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCt                                                                                  >  1:437522/1‑150 (MQ=255)
                                                                  aGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCt                                                                                  <  1:1730309/150‑1 (MQ=255)
                                                                   gCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTa                                                                                 >  2:1591497/1‑150 (MQ=255)
                                                                   gCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTa                                                                                 >  1:1387978/1‑150 (MQ=255)
                                                                          cTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCgg                                                                          >  1:1176682/1‑150 (MQ=255)
                                                                           tATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGGCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCggg                                                                         <  2:1047025/150‑1 (MQ=255)
                                                                            aTGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgc                                                                        >  2:1323792/1‑150 (MQ=255)
                                                                             tGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgcg                                                                       <  1:1692914/150‑1 (MQ=255)
                                                                              taCGAACTGGTTGCCGATGAACGCATTCGCAAAAAATCCCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCAACATCAAGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGgcgc                                                                      <  2:1523215/149‑1 (MQ=255)
                                                                               aCGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCa                                                                     >  2:674985/1‑150 (MQ=255)
                                                                                  aaCTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGCAGCCACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATaa                                                                  >  2:1214650/1‑150 (MQ=255)
                                                                                          gCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGt                                                          >  1:3885/1‑150 (MQ=255)
                                                                                          gCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGt                                                          <  1:1418959/150‑1 (MQ=255)
                                                                                            cGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTaa                                                        >  1:650871/1‑150 (MQ=255)
                                                                                            cGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTaa                                                        >  1:1525845/1‑150 (MQ=255)
                                                                                                 aaCGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGACGGGCGCATAAATATGAGTAATctct                                                   >  2:105193/1‑150 (MQ=255)
                                                                                                  aCGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTACCCCTATCGCCGGGCGCATAAATATGAGTAATCTCTg                                                  >  2:745235/1‑150 (MQ=255)
                                                                                                  aCGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTg                                                  >  2:188833/1‑150 (MQ=255)
                                                                                                   cGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGc                                                 <  1:859751/150‑1 (MQ=255)
                                                                                                    gCATTCGCAAAAAATACCTCAACGCCCTTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCt                                                <  1:636984/150‑1 (MQ=255)
                                                                                                     cATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCTACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTc                                               >  2:1222234/1‑150 (MQ=255)
                                                                                                     cATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTc                                               >  1:523532/1‑150 (MQ=255)
                                                                                                      aTTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCACTACATCATGTATGAAGACACGGTAAACCTTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCtaa                                              >  1:1418865/1‑148 (MQ=255)
                                                                                                         cGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGaa                                           >  1:1370084/1‑150 (MQ=255)
                                                                                                          gCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGaaa                                          >  1:555721/1‑150 (MQ=255)
                                                                                                             aaaaaTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAAttt                                       >  1:181155/1‑150 (MQ=255)
                                                                                                              aaaaTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAAtttt                                      <  1:294913/150‑1 (MQ=255)
                                                                                                                aaTACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGTGTAATCTCTGCTCAGAAATTTTGc                                    <  2:1052634/150‑1 (MQ=255)
                                                                                                                  tACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAg                                  <  1:1528370/150‑1 (MQ=255)
                                                                                                                  tACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGAAACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTACTCTCTGCTCAGAAATTTTGCAg                                  >  1:1028338/1‑150 (MQ=255)
                                                                                                                   aCCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAgg                                 <  1:1587193/150‑1 (MQ=255)
                                                                                                                   aCCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAgg                                 <  1:339414/150‑1 (MQ=255)
                                                                                                                     cTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATTAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGtt                               >  2:351364/1‑150 (MQ=255)
                                                                                                                     cTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGtt                               >  1:914036/1‑150 (MQ=255)
                                                                                                                     cTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATAGGCGGGCGCATTAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGtt                               >  2:320629/1‑150 (MQ=255)
                                                                                                                      tacaCGCCCGTTATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCTTAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTa                              <  2:254490/147‑1 (MQ=255)
                                                                                                                       cAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGATTCCGGCTATCCCTACAGCATGTATGCAGCCACGGTAAACCGTGCTCCCCCGCTCGCCGGGCGCAGAAATCGGAGGACGCGCTGCTCAGAAAGGTTGCAGGTGaa                             >  2:836258/1‑150 (MQ=255)
                                                                                                                       cAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTaa                             >  2:1209554/1‑150 (MQ=255)
                                                                                                                        aaCGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAAc                            <  2:460596/150‑1 (MQ=255)
                                                                                                                        aaCGCCCGTGATTTATTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAAc                            <  1:209999/150‑1 (MQ=255)
                                                                                                                          cGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAg                          <  2:442086/150‑1 (MQ=255)
                                                                                                                             ccGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATCAAGACACGGTAAACCGTTCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAgcgc                       >  2:640991/1‑150 (MQ=255)
                                                                                                                             actggATTTCTGCCAGTCGCTGGAAGATATCCAGTTTGAGTCCGGCTATCACTACATCATGTATGAAGACACGGTATACAGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTAAGAAATTTTGCAGGTTATCAgcgc                       >  1:53189/5‑150 (MQ=255)
                                                                                                                             actggATTTCTGCCAGTCGCTGGAAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAgcgc                       >  1:313432/5‑150 (MQ=255)
                                                                                                                                  aTTTCTGCCAGTCGCTGGAAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCag                  >  1:1389673/1‑150 (MQ=255)
                                                                                                                                   tttCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCaga                 >  2:963263/1‑150 (MQ=255)
                                                                                                                                    ttcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCagag                <  2:1387978/150‑1 (MQ=255)
                                                                                                                                    ttcttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCagag                <  2:722142/150‑1 (MQ=255)
                                                                                                                                      cttcCAGCGACTGGGAGAAATCCAGGTGGAGTCCGGCTATCCCTACATCATGTATGCAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTa              <  2:1476192/150‑1 (MQ=255)
                                                                                                                                      cttcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTa              >  2:686396/1‑150 (MQ=255)
                                                                                                                                        tcCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATTAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGATAACAGCGCCTCAGAGTATg            >  1:1176925/1‑150 (MQ=255)
                                                                                                                                          cAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGAc          <  2:905583/150‑1 (MQ=255)
                                                                                                                                           cgCGAGTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGACGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACg         <  2:292547/149‑1 (MQ=255)
                                                                                                                                           cgCGACTGGCAGAAATCCAGTTTGAGTCCGGATATCCCTACATCATGGATGTAGAAACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACg         <  2:269878/149‑1 (MQ=255)
                                                                                                                                           aGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACg         <  2:519343/150‑1 (MQ=255)
                                                                                                                                            gtcgcTGGAAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACga        <  2:265997/146‑1 (MQ=255)
                                                                                                                                            gCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCAAGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACga        >  1:108804/1‑150 (MQ=255)
                                                                                                                                              gACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACgaga      <  2:1358796/150‑1 (MQ=255)
                                                                                                                                               aCTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACgagat     >  2:1436925/1‑149 (MQ=255)
                                                                                                                                                cTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTATGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAAt    >  2:1287982/1‑150 (MQ=255)
                                                                                                                                                  ggCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAAtct  >  1:737999/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGCGCTGTTTTCGCCTTATGACGTAGAGCGAGTTTATGGCAAGCCGTTTGCCGATGTGGCCATCAGCCAACACTATGACGAACTGGTTGCCGATGAACGCATTCGCAAAAAATACCTCAACGCCCGTGATTTCTTCCAGCGACTGGCAGAAATCCAGTTTGAGTCCGGCTATCCCTACATCATGTATGAAGACACGGTAAACCGTGCTAACCCTATCGCCGGGCGCATAAATATGAGTAATCTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCT  >  NC_000913/2802263‑2802558

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: