Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,659,603 A→T 6.9% Y512F (TAT→TTT)  ynfE → putative selenate reductase YnfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,659,6030AT 6.9% 406.9 / 14.3 131Y512F (TAT→TTT) ynfEputative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base A (62/60);  new base T (4/5);  total (66/65)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 138 total aligned reads displayed.

TGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTG  >  NC_000913/1659454‑1659752
                                                                                                                                                     |                                                                                                                                                     
tgtgCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAata                                                                                                                                                       >  2:765194/1‑150 (MQ=255)
  tgCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTAGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCACGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCCGCAAAATATg                                                                                                                                                     >  1:1779496/1‑150 (MQ=255)
   gCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCTACAACACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGc                                                                                                                                                    >  2:1273200/1‑150 (MQ=255)
   gCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGc                                                                                                                                                    >  2:525371/1‑150 (MQ=255)
      cATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGa                                                                                                                                                 <  2:1738650/150‑1 (MQ=255)
       aTTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGAc                                                                                                                                                >  2:1642898/1‑150 (MQ=255)
        ttAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCAACGAAAACTTTATGACCTCATCGGCATATTTTGCCGatg                                                                                                                                               <  2:1310199/150‑3 (MQ=255)
         tAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCATATTTTGCCGatga                                                                                                                                              <  1:1950051/150‑4 (MQ=255)
          aaGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATCGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACAtt                                                                                                                                             >  2:666289/1‑150 (MQ=255)
           aGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCATATTTTGCCGatgagg                                                                                                                                            <  2:11500/150‑6 (MQ=255)
           aGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACAtta                                                                                                                                            >  2:1332820/1‑149 (MQ=255)
           aGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTc                                                                                                                                            >  1:1591418/1‑150 (MQ=255)
            gTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTct                                                                                                                                           <  2:66029/150‑1 (MQ=255)
               cATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTctgct                                                                                                                                        >  2:1194965/1‑150 (MQ=255)
                 tCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCtgctgc                                                                                                                                      >  2:2349185/1‑150 (MQ=255)
                  cTGGAACTACGCGGGTAATACCCTCGTTAATCAGCATTCTCACATCAACAAAACGCATGAAATTCTGCAGGACGAATGGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGcc                                                                                                                                     <  2:468984/150‑1 (MQ=255)
                  cTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGcc                                                                                                                                     <  2:737038/150‑1 (MQ=255)
                  cTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGcc                                                                                                                                     <  1:2273587/150‑1 (MQ=255)
                   tGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCa                                                                                                                                    >  1:2236199/1‑150 (MQ=255)
                    ggAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAg                                                                                                                                   >  1:2188889/1‑150 (MQ=255)
                    ggAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAg                                                                                                                                   >  1:723431/1‑150 (MQ=255)
                     gAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGa                                                                                                                                  <  1:966646/150‑1 (MQ=255)
                            gcgGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCtgatg                                                                                                                           >  2:294270/1‑150 (MQ=255)
                            gcgGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCtgatg                                                                                                                           >  1:1180207/1‑150 (MQ=255)
                             cgGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCtgatga                                                                                                                          >  1:1815681/1‑150 (MQ=255)
                              gggAAATCCCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGAc                                                                                                                         >  1:1326887/1‑150 (MQ=255)
                                    tACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGa                                                                                                                   <  2:1025953/150‑1 (MQ=255)
                                    tACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGa                                                                                                                   <  1:1039704/150‑1 (MQ=255)
                                    tACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGa                                                                                                                   <  1:1674356/150‑1 (MQ=255)
                                     aCCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAg                                                                                                                  <  1:2287259/150‑1 (MQ=255)
                                     aCCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAg                                                                                                                  <  1:1005017/150‑1 (MQ=255)
                                      cccTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGc                                                                                                                 <  2:2384608/150‑1 (MQ=255)
                                      cccTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGc                                                                                                                 >  2:738307/1‑150 (MQ=255)
                                       ccTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCa                                                                                                                <  1:2162319/150‑1 (MQ=255)
                                        cTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAg                                                                                                               <  1:1311902/150‑1 (MQ=255)
                                        cTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAg                                                                                                               <  1:1194965/150‑1 (MQ=255)
                                          cGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGa                                                                                                             >  2:550080/1‑150 (MQ=255)
                                                 cAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtatt                                                                                                      >  1:2234222/1‑150 (MQ=255)
                                                   gCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTGATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtattat                                                                                                    <  2:1368871/150‑1 (MQ=255)
                                                      ttCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTcc                                                                                                 <  2:1858849/150‑1 (MQ=255)
                                                       tCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCt                                                                                                <  1:44126/150‑1 (MQ=255)
                                                        cTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTa                                                                                               >  1:659985/1‑150 (MQ=255)
                                                         tGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTaa                                                                                              >  1:2104982/1‑150 (MQ=255)
                                                          taCATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAAc                                                                                             <  2:78769/149‑1 (MQ=255)
                                                             aTCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGAc                                                                                          <  1:1332820/150‑1 (MQ=255)
                                                             aTCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGAc                                                                                          <  2:2002665/150‑1 (MQ=255)
                                                                aacaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTAc                                                                                       >  1:354815/1‑150 (MQ=255)
                                                                aacaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTAc                                                                                       >  2:1133196/1‑150 (MQ=255)
                                                                aacaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGAAGGAAGATATTATTCCTAACGACTAc                                                                                       >  1:1655020/1‑150 (MQ=255)
                                                                 acaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACg                                                                                      >  2:1979568/1‑150 (MQ=255)
                                                                  caaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGc                                                                                     >  1:1446853/1‑150 (MQ=255)
                                                                   aaaaCGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGcc                                                                                    >  2:2253912/1‑150 (MQ=255)
                                                                       cGCATGAAATTCTGCAGGACTAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCAGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTa                                                                                <  1:142190/150‑1 (MQ=255)
                                                                         cATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAAc                                                                              >  2:579290/1‑150 (MQ=255)
                                                                         cATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAAc                                                                              >  2:1391866/1‑150 (MQ=255)
                                                                          aTGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACa                                                                             >  2:207486/1‑150 (MQ=255)
                                                                           tGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGCGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACAt                                                                            >  1:1632212/1‑150 (MQ=255)
                                                                           tGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACAt                                                                            >  2:1859882/1‑150 (MQ=255)
                                                                              aaTTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATggg                                                                         <  2:1591418/150‑1 (MQ=255)
                                                                                          cGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCAtttt                                                             <  2:2179018/150‑1 (MQ=255)
                                                                                            aaTCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTtcc                                                           >  2:2380759/1‑150 (MQ=255)
                                                                                            aaTCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTtcc                                                           >  1:1163194/1‑150 (MQ=255)
                                                                                                 aaatgcgaaatgATCGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAg                                                      <  1:666289/150‑1 (MQ=255)
                                                                                                  aatgcgaaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGc                                                     <  2:177687/150‑1 (MQ=255)
                                                                                                   atgcgaaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGcc                                                    <  2:1293560/150‑1 (MQ=255)
                                                                                                     gcgaaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTg                                                  <  1:2068705/150‑1 (MQ=255)
                                                                                                      cgaaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGt                                                 >  2:47240/1‑150 (MQ=255)
                                                                                                        aaatgATTGTGGTCATCGAATACTTTATGTCCTCATAGGCAAAATATGCCGACATTCTGCTGCCAGACGTGATGACCGTTGTGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCa                                               <  2:272483/150‑1 (MQ=255)
                                                                                                        aaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCa                                               <  1:742839/150‑1 (MQ=255)
                                                                                                        aaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCa                                               >  2:1125340/1‑150 (MQ=255)
                                                                                                         aatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCAc                                              >  2:169113/1‑150 (MQ=255)
                                                                                                             aTTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGc                                          >  1:233133/1‑150 (MQ=255)
                                                                                                                gtgGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGaa                                       <  1:2093400/150‑1 (MQ=255)
                                                                                                                  ggTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGaaaa                                     >  1:382280/1‑150 (MQ=255)
                                                                                                                  ggTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGaaaa                                     >  1:391046/1‑150 (MQ=255)
                                                                                                                       tCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCg                                <  1:1979568/150‑1 (MQ=255)
                                                                                                                         gAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGaa                              >  1:255061/1‑150 (MQ=255)
                                                                                                                         gAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCCCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGaa                              >  2:1767649/1‑150 (MQ=255)
                                                                                                                          aaaaCTTTATGACCTCATCGGCAAAATATGCCGAGATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAAc                             <  1:1216394/150‑1 (MQ=255)
                                                                                                                          aaaaCTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATATAATTTTACTCCAGACTGTCACCAGCGAAAAATTCGAAc                             >  1:950506/1‑150 (MQ=255)
                                                                                                                           aaaCTTTATGACCTCATCGACAAAATATGCCGACATTCTGCTGCCAGAACTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAAAATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACg                            >  1:1424455/1‑150 (MQ=255)
                                                                                                                              cTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCaa                         <  2:1180207/150‑1 (MQ=255)
                                                                                                                                 tATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAAcc                      <  1:620253/150‑1 (MQ=255)
                                                                                                                                     agaatgtCGGCATATTTTGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGAtt                  <  2:1384244/144‑1 (MQ=255)
                                                                                                                                      ccTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGTCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGAttt                 <  1:738307/150‑1 (MQ=255)
                                                                                                                                        tCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTAc               <  1:78372/150‑1 (MQ=255)
                                                                                                                                        atgtCGGCATATTTTGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTAc               >  2:2234221/4‑150 (MQ=255)
                                                                                                                                        atgtCGGCATATTTTGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTAc               >  2:371243/4‑150 (MQ=255)
                                                                                                                                        atgtCGGCATATTTTGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTAc               >  1:853542/4‑150 (MQ=255)
                                                                                                                                          gtCGGCATATTTTGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTAATCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTg             >  1:295883/2‑150 (MQ=255)
                                                                                                                                          aTCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTg             <  1:1237592/150‑1 (MQ=255)
                                                                                                                                            cGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGa           >  2:1405603/1‑150 (MQ=255)
                                                                                                                                            cGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGa           <  1:2146539/150‑1 (MQ=255)
                                                                                                                                             ggCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGAt          <  1:1928221/150‑1 (MQ=255)
                                                                                                                                              gCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATc         >  1:1536532/1‑150 (MQ=255)
                                                                                                                                                    tatGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGt   <  2:2236199/150‑1 (MQ=255)
                                                                                                                                                    tatGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGt   <  2:2206388/150‑1 (MQ=255)
                                                                                                                                                     atGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTg  <  2:723431/150‑1 (MQ=255)
                                                                                                                                                     atGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTg  <  1:423913/150‑1 (MQ=255)
                                                                                                                                                     atGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTg  >  1:1084499/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTG  >  NC_000913/1659454‑1659752

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: