Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 325,853 A→G R465R (CGT→CGC)  betA ← choline dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913325,8530AG99.1% 401.6 / ‑4.6 112R465R (CGT→CGC) betAcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (56/55);  total (56/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

Only 100 of 113 total aligned reads displayed.

ATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCC  >  NC_000913/325716‑325986
                                                                                                                                         |                                                                                                                                     
aTCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCg                                                                                                                                     <  1:30035/140‑1 (MQ=255)
aTCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCg                                                                                                                                     <  2:504919/140‑1 (MQ=255)
   tGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCAc                                                                                                                                  <  1:1068631/140‑1 (MQ=255)
       gCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAAc                                                                                                                              <  1:365120/140‑1 (MQ=255)
       gCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAAc                                                                                                                              >  2:289521/1‑140 (MQ=255)
             tCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCg                                                                                                                        >  1:687170/1‑140 (MQ=255)
               gACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAg                                                                                                                      <  2:1007753/140‑1 (MQ=255)
               gACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAg                                                                                                                      <  2:704312/140‑1 (MQ=255)
               gACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAg                                                                                                                      <  2:147060/140‑1 (MQ=255)
                aCGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACCGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGc                                                                                                                     <  2:1071946/140‑1 (MQ=255)
                            cacGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTc                                                                                                         <  2:322597/140‑1 (MQ=255)
                               gCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTgg                                                                                                      >  1:330503/1‑140 (MQ=255)
                                 aGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCa                                                                                                    <  2:57284/140‑1 (MQ=255)
                                 aGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCa                                                                                                    >  1:57284/1‑140 (MQ=255)
                                      ttCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGa                                                                                               >  1:432650/1‑140 (MQ=255)
                                       tCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGac                                                                                              >  2:824163/1‑140 (MQ=255)
                                        cTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTg                                                                                                       <  2:774730/130‑1 (MQ=255)
                                        cTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTg                                                                                                       >  1:774730/1‑130 (MQ=255)
                                        cTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGaca                                                                                             >  1:1096862/1‑140 (MQ=255)
                                         tAACCCGTGTACGCGGTCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGATTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCAt                                                                                                              >  2:466896/1‑122 (MQ=255)
                                         tAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCAt                                                                                                              >  2:842136/1‑122 (MQ=255)
                                         tAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCAt                                                                                                              <  1:842136/122‑1 (MQ=255)
                                         tAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGacac                                                                                            <  2:345310/140‑1 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTg                                                                                                       >  2:533037/1‑122 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTg                                                                                                       <  1:533037/122‑1 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCAtt                                                                                                  <  1:991330/127‑1 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCt                                                                                     >  1:1191825/1‑140 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCt                                                                                     >  2:1176498/1‑140 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCt                                                                                     >  2:14168/1‑140 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCt                                                                                     >  2:593825/1‑140 (MQ=255)
                                                tgtACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTCGCATTCGACACCGGGGCt                                                                                     >  2:310717/1‑140 (MQ=255)
                                                       ggCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTcg                                                                              >  2:969061/1‑140 (MQ=255)
                                                        gCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTcgc                                                                             <  2:1096862/140‑1 (MQ=255)
                                                        gCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTcgc                                                                             <  1:729344/140‑1 (MQ=255)
                                                         ccTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTcgcg                                                                            >  1:919915/1‑140 (MQ=255)
                                                         ccTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGAGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTcgcg                                                                            >  2:435198/1‑140 (MQ=255)
                                                          cTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCgg                                                                           >  1:464015/1‑140 (MQ=255)
                                                          cTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCgg                                                                           >  2:569951/1‑140 (MQ=255)
                                                          cTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCgg                                                                           >  1:412618/1‑140 (MQ=255)
                                                           ttCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGc                                                                          <  2:330503/140‑1 (MQ=255)
                                                             cGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCa                                                                        <  1:1176498/140‑1 (MQ=255)
                                                             cGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCa                                                                        <  1:574133/140‑1 (MQ=255)
                                                             cGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCa                                                                        <  2:303585/140‑1 (MQ=255)
                                                              gCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCAc                                                                       >  2:160533/1‑140 (MQ=255)
                                                                cGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGa                                                                     <  1:824163/140‑1 (MQ=255)
                                                                   cAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCAc                                                                       <  1:85460/135‑1 (MQ=255)
                                                                   cAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCAc                                                                       <  1:769043/135‑1 (MQ=255)
                                                                   cAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCAc                                                                       >  2:85460/1‑135 (MQ=255)
                                                                   cAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCAc                                                                       >  2:769043/1‑135 (MQ=255)
                                                                   cAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATAc                                                                  >  2:1207127/1‑140 (MQ=255)
                                                                        aCGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATc                                                             <  2:822213/140‑1 (MQ=255)
                                                                             cATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAgcg                                                        <  1:478295/140‑1 (MQ=255)
                                                                               tctcGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGcgcg                                                      >  2:50971/1‑140 (MQ=255)
                                                                                ctcGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCgg                                                     >  2:416538/1‑140 (MQ=255)
                                                                                ctcGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCgg                                                     >  1:915835/1‑140 (MQ=255)
                                                                                 tcgtcgTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCggg                                                    >  2:720737/1‑140 (MQ=255)
                                                                                  cgtcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGt                                                   <  2:464015/140‑1 (MQ=255)
                                                                                  cgtcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGt                                                   <  2:120817/140‑1 (MQ=255)
                                                                                  cgtcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGt                                                   <  1:91943/140‑1 (MQ=255)
                                                                                  cgtcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGt                                                   <  1:569951/140‑1 (MQ=255)
                                                                                  cgtcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGt                                                   <  1:720737/140‑1 (MQ=255)
                                                                                    tcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTg                                                 <  2:1220728/140‑1 (MQ=255)
                                                                                    tcgtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTg                                                 >  2:1193647/1‑140 (MQ=255)
                                                                                      gtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGAt                                               >  2:794565/1‑140 (MQ=255)
                                                                                      gtAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGAt                                               >  1:318097/1‑140 (MQ=255)
                                                                                        aaCCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGc                                             <  1:870181/140‑1 (MQ=255)
                                                                                          cccATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCAt                                           >  2:591421/1‑140 (MQ=255)
                                                                                          cccATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCAt                                           >  1:43152/1‑140 (MQ=255)
                                                                                          cccATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCAt                                           <  1:14168/140‑1 (MQ=255)
                                                                                              ttttGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGt                                                   <  1:236207/128‑1 (MQ=255)
                                                                                               tttGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGAtct                                      >  2:1173767/1‑140 (MQ=255)
                                                                                                     ggTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCggg                                >  2:73803/1‑140 (MQ=255)
                                                                                                     ggTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCggg                                <  1:160533/140‑1 (MQ=255)
                                                                                                     ggTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCggg                                >  1:1010271/1‑140 (MQ=255)
                                                                                                     ggTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCgg                                 <  1:1212611/139‑1 (MQ=255)
                                                                                                     ggTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCgg                                 >  2:1212611/1‑139 (MQ=255)
                                                                                                      gTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGt                               <  1:416538/140‑1 (MQ=255)
                                                                                                                gATGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCTCGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAAtt                     >  2:478221/1‑140 (MQ=255)
                                                                                                                 aTGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGGTGTTCATCCGTCTGGCAATCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGGGGGTTGATGCATGATCTCGCGGGt                               <  2:837075/129‑1 (MQ=255)
                                                                                                                 aTGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGt                               >  1:598777/1‑129 (MQ=255)
                                                                                                                 aTGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGt                               >  1:837075/1‑129 (MQ=255)
                                                                                                                 aTGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGt                               <  2:598777/129‑1 (MQ=255)
                                                                                                                 aTGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTg                    >  2:256446/1‑140 (MQ=255)
                                                                                                                 aTGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTg                    >  2:1187812/1‑140 (MQ=255)
                                                                                                                  tGGAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGc                   <  2:687170/140‑1 (MQ=255)
                                                                                                                   ggAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTcg                                    <  2:436436/122‑1 (MQ=255)
                                                                                                                   ggAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTcg                                    >  1:436436/1‑122 (MQ=255)
                                                                                                                    gAAGGCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGt                 <  1:256446/140‑1 (MQ=255)
                                                                                                                        gCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTcgcg             <  1:1193647/140‑1 (MQ=255)
                                                                                                                        gCGGTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTcgcg             <  1:50971/140‑1 (MQ=255)
                                                                                                                           gTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTg                    >  2:661789/1‑130 (MQ=255)
                                                                                                                           gTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTg                    <  1:661789/130‑1 (MQ=255)
                                                                                                                           gTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGaa          >  1:452016/1‑140 (MQ=255)
                                                                                                                           gTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGaa          <  1:1187812/140‑1 (MQ=255)
                                                                                                                           gTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGaa          >  1:545997/1‑140 (MQ=255)
                                                                                                                           gTTTCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGaa          <  1:1173767/140‑1 (MQ=255)
                                                                                                                             ttCGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACt        <  1:607024/140‑1 (MQ=255)
                                                                                                                               cGGCGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTcc      >  2:646128/1‑140 (MQ=255)
                                                                                                                                  cGTGGTTGCGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGc   <  2:432650/140‑1 (MQ=255)
                                                                                                                                   gTGGTTGCGCACGAACTCAACGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGcc  >  2:255226/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                     
ATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCATCGAGCTGTTCATCCGTCTGGCATTCGACACCGGGGCTGATTTCGCGGCCACGATACTGATCCAGCGCGGGTTGATGCATGATCTCGCGGGTGATGCGAATTGCGTCGCGGAACTCCTGCC  >  NC_000913/325716‑325986

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