Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 531,581 A→G 8.6% L214S (TTG→TCG)  allS ← DNA‑binding transcriptional activator AllS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,5810AG 8.6% 268.3 / 11.8 93L214S (TTG→TCG) allSDNA‑binding transcriptional activator AllS
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base G (3/5);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCT  >  NC_000913/531433‑531730
                                                                                                                                                    |                                                                                                                                                     
cATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCaa                                                                                                                                                      >  1:530393/1‑150 (MQ=255)
   cGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCg                                                                                                                                                   <  2:285469/150‑1 (MQ=255)
   cGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCg                                                                                                                                                   <  2:843623/150‑1 (MQ=255)
      tGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCa                                                                                                                                                <  1:337598/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)
         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:1648/150‑1 (MQ=255)
                                    aTCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCt                                                                                                                  <  2:7549/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                                                                                                                                   <  1:552350/110‑1 (MQ=255)
                                                           aTTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTg                                                                                                                                   >  2:552350/1‑110 (MQ=255)
                                                                 gCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgcc                                                                                     >  2:1227124/1‑150 (MQ=255)
                                                                  cAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgccg                                                                                    >  2:1173569/1‑150 (MQ=255)
                                                                             tGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCtt                                                                         >  2:368658/1‑150 (MQ=255)
                                                                             tGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCtt                                                                         <  2:530393/150‑1 (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                                                              <  2:1264762/150‑1 (MQ=255)
                                                                                        gccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGc                                                              >  1:183302/1‑150 (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                                                 <  1:529679/150‑1 (MQ=255)
                                                                                                     cgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCg                                                 >  1:704809/1‑150 (MQ=255)
                                                                                                     cgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCg                                                 >  2:849542/1‑150 (MQ=255)
                                                                                                       aTTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGcc                                               >  1:632619/1‑150 (MQ=255)
                                                                                                       aTTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGcc                                               <  2:183302/150‑1 (MQ=255)
                                                                                                        ttttCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCa                                              >  1:474968/1‑150 (MQ=255)
                                                                                                            cGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCgg                                          >  1:363118/1‑150 (MQ=255)
                                                                                                                  cATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATc                                    >  2:893416/1‑150 (MQ=255)
                                                                                                                  cATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATc                                    <  2:245186/150‑1 (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:207562/1‑150 (MQ=255)
                                                                                                                               ataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACa                       <  1:1227124/150‑1 (MQ=255)
                                                                                                                               ataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACa                       >  2:201732/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)
                                                                                                                                                    |                                                                                                                                                     
CATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCT  >  NC_000913/531433‑531730

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: