Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 531,584 C→T 8.4% R213Q (CGA→CAA)  allS ← DNA‑binding transcriptional activator AllS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,5840CT 8.4% 288.3 / 10.1 95R213Q (CGA→CAA) allSDNA‑binding transcriptional activator AllS
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base T (3/5);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAG  >  NC_000913/531436‑531733
                                                                                                                                                    |                                                                                                                                                     
cGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCg                                                                                                                                                      <  2:843623/150‑1 (MQ=255)
cGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCg                                                                                                                                                      <  2:285469/150‑1 (MQ=255)
   tGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCa                                                                                                                                                   <  2:788914/150‑1 (MQ=255)
   tGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCa                                                                                                                                                   >  1:993594/1‑150 (MQ=255)
   tGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCa                                                                                                                                                   <  1:145607/150‑1 (MQ=255)
   tGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCa                                                                                                                                                   <  1:337598/150‑1 (MQ=255)
      gATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGc                                                                                                                                                <  1:233317/150‑1 (MQ=255)
       aTTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggc                                                                                                                                               <  1:119010/150‑3 (MQ=255)
       aTTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggc                                                                                                                                               <  2:559166/150‑3 (MQ=255)
        ttACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggca                                                                                                                                              <  2:1018218/150‑4 (MQ=255)
         tACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAATCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAc                                                                                                                                             <  1:330415/150‑1 (MQ=255)
         tACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAc                                                                                                                                             <  1:835367/150‑1 (MQ=255)
          aCCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcg                                                                                                                                            >  1:451599/1‑150 (MQ=255)
           cccGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggcaaaa                                                                                                                                           <  1:1409637/150‑7 (MQ=255)
            ccGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTCTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggcaaaaa                                                                                                                                          <  1:518070/150‑8 (MQ=255)
                 tGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTgg                                                                                                                                     <  1:276112/150‑1 (MQ=255)
                     cAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTaa                                                                                                                                 >  2:120485/1‑150 (MQ=255)
                       tttgttgATTATCGATCATTGACTGGCAAATCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAgg                                                                                                                               <  1:721661/149‑1 (MQ=255)
                       gttgttGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAgg                                                                                                                               >  2:1177660/1‑150 (MQ=255)
                         tgttgATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTg                                                                                                                             <  1:1444404/150‑1 (MQ=255)
                             gATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCggg                                                                                                                         <  1:971032/150‑1 (MQ=255)
                               ttATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGgcg                                                                                                                       <  2:1317996/150‑1 (MQ=255)
                                 aTCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCt                                                                                                                     <  2:7549/150‑1 (MQ=255)
                                 aTCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCt                                                                                                                     <  2:1648/150‑1 (MQ=255)
                                  tCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTg                                                                                                                    <  2:1284664/150‑1 (MQ=255)
                                     aTCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCt                                                                                                                 <  2:451599/150‑1 (MQ=255)
                                      tCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCtt                                                                                                                >  1:1320972/1‑150 (MQ=255)
                                       cATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTc                                                                                                               <  1:515425/150‑1 (MQ=255)
                                         ttGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCaa                                                                                                             >  1:390830/1‑150 (MQ=255)
                                           gACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAata                                                                                                           >  2:1175211/1‑150 (MQ=255)
                                             cTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATAtt                                                                                                         >  2:5325/1‑150 (MQ=255)
                                               ggCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGa                                                                                                       >  1:938929/1‑150 (MQ=255)
                                               ggCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATAACGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGa                                                                                                       >  1:6112/1‑150 (MQ=255)
                                                gCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCAAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGAc                                                                                                      >  2:1312214/1‑150 (MQ=255)
                                                      cGATTTTGGCAAAAACCCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCgg                                                                                                >  1:1040033/1‑150 (MQ=255)
                                                        aTTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTg                                                                                                                                      >  2:552350/1‑110 (MQ=255)
                                                        aTTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTg                                                                                                                                      <  1:552350/110‑1 (MQ=255)
                                                              gCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgcc                                                                                        >  2:1227124/1‑150 (MQ=255)
                                                               cAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgccg                                                                                       >  2:1173569/1‑150 (MQ=255)
                                                                          tGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCtt                                                                            <  2:530393/150‑1 (MQ=255)
                                                                          tGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCtt                                                                            >  2:368658/1‑150 (MQ=255)
                                                                           gCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTc                                                                           <  2:1320972/150‑1 (MQ=255)
                                                                           gCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTc                                                                           <  1:80756/150‑1 (MQ=255)
                                                                            ccAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTCACCGCCGGAAAGCGCCGCAACTGAGCTTCt                                                                          >  2:1085550/1‑150 (MQ=255)
                                                                             cAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGTTGCGGGCGCTGTCTGCAATATTGACCGCCTGAAAGCGCCGCAACTGCTCTTCTg                                                                         <  2:1024799/150‑1 (MQ=255)
                                                                               aCGCCCGCCAGAGGGGCGGCGATTTTCGTTTCCATATCGGGAAAAATAATCTCTTTTTGCCCTTGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGtt                                                                       <  2:1292941/150‑1 (MQ=255)
                                                                                 gcccgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAg                                                                     <  1:1177660/150‑1 (MQ=255)
                                                                                  cccgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGc                                                                    >  1:117629/1‑150 (MQ=255)
                                                                                   ccgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCg                                                                   >  1:1401986/1‑150 (MQ=255)
                                                                                   ccgccAGATGGGCGGCGATTTTCGTTTCAATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCg                                                                   >  2:561378/1‑150 (MQ=255)
                                                                                     gccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGc                                                                 >  1:183302/1‑150 (MQ=255)
                                                                                     gccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGc                                                                 <  2:1264762/150‑1 (MQ=255)
                                                                                        aGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGctct                                                              >  1:681275/1‑150 (MQ=255)
                                                                                        aGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGctct                                                              >  2:654791/1‑150 (MQ=255)
                                                                                          aTGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTc                                                            <  1:705301/150‑1 (MQ=255)
                                                                                          aTGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTc                                                            >  2:1563318/1‑150 (MQ=255)
                                                                                           tGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCa                                                           <  2:32283/150‑1 (MQ=255)
                                                                                           tGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCa                                                           <  1:1612053/150‑1 (MQ=255)
                                                                                                  cgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCg                                                    >  2:849542/1‑150 (MQ=255)
                                                                                                  cgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCg                                                    <  1:529679/150‑1 (MQ=255)
                                                                                                  cgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCg                                                    >  1:704809/1‑150 (MQ=255)
                                                                                                    aTTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGcc                                                  <  2:183302/150‑1 (MQ=255)
                                                                                                    aTTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGcc                                                  >  1:632619/1‑150 (MQ=255)
                                                                                                     ttttCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCa                                                 >  1:474968/1‑150 (MQ=255)
                                                                                                         cGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCgg                                             >  1:363118/1‑150 (MQ=255)
                                                                                                               cATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATc                                       <  2:245186/150‑1 (MQ=255)
                                                                                                               cATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATc                                       >  2:893416/1‑150 (MQ=255)
                                                                                                                  atCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgc                                    >  1:1042022/1‑150 (MQ=255)
                                                                                                                   tCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgcc                                   <  2:457799/150‑1 (MQ=255)
                                                                                                                    cGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgccg                                  >  2:973544/1‑150 (MQ=255)
                                                                                                                    cGGGAACAATAATCTCTTTTTGCCCTGGCAATCGACAGGCGACCCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAAAATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgccg                                  >  2:86745/1‑150 (MQ=255)
                                                                                                                    cGGGAACAATAATCTCTTTTTGCCCTGACAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgccg                                  >  2:1458822/1‑150 (MQ=255)
                                                                                                                      ggAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCgccgcc                                <  2:305554/150‑1 (MQ=255)
                                                                                                                      ggAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCgccgcc                                <  2:363118/150‑1 (MQ=255)
                                                                                                                         aCAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATg                             >  1:621956/1‑150 (MQ=255)
                                                                                                                          cAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGa                            >  2:811898/1‑150 (MQ=255)
                                                                                                                          cAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGa                            >  1:875716/1‑150 (MQ=255)
                                                                                                                          cAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGa                            >  1:742500/1‑150 (MQ=255)
                                                                                                                            ataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACa                          >  2:201732/1‑150 (MQ=255)
                                                                                                                            ataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACa                          <  1:1227124/150‑1 (MQ=255)
                                                                                                                            ataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACa                          >  2:207562/1‑150 (MQ=255)
                                                                                                                                  tctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAACgcgc                    >  1:1437501/1‑150 (MQ=255)
                                                                                                                                    tttttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCa                  >  2:873852/1‑150 (MQ=255)
                                                                                                                                      tttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCAtt                <  1:1540556/150‑1 (MQ=255)
                                                                                                                                      cgtcgccTGGCGATTGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCAtt                >  1:304362/6‑150 (MQ=255)
                                                                                                                                       ttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTg               <  1:1075102/150‑1 (MQ=255)
                                                                                                                                        tcgccTGGCGATTGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGc              >  1:1138821/4‑150 (MQ=255)
                                                                                                                                          gccTGGCGATTGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCAc            >  1:200877/2‑150 (MQ=255)
                                                                                                                                          cccTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCAc            >  2:1190187/1‑150 (MQ=255)
                                                                                                                                           ccTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCAcc           <  2:938929/150‑1 (MQ=255)
                                                                                                                                             tcgCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGCATGATCCGCCGCCATGACAAAGCGCCATTGCACCGa         >  1:388260/3‑150 (MQ=255)
                                                                                                                                             tGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGa         <  2:704809/150‑1 (MQ=255)
                                                                                                                                                cAATCGCCAGGCGACGCTTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATcc      <  2:1084553/150‑1 (MQ=255)
                                                                                                                                                 aaTCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCt     <  1:134433/150‑1 (MQ=255)
                                                                                                                                                 aaTCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCt     >  2:1242132/1‑150 (MQ=255)
                                                                                                                                                    cGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGTACCGATCCTAAg  <  2:744535/150‑1 (MQ=255)
                                                                                                                                                    cGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAg  <  1:487691/150‑1 (MQ=255)
                                                                                                                                                    cGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAg  <  2:993594/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAG  >  NC_000913/531436‑531733

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