Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,281,666 A→T 6.5% intergenic (+37/+117) yjcE → / ← yjcF putative transporter YjcE/pentapeptide repeat‑containing protein YjcF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,6660AT 6.5% 313.4 / 9.8 107intergenic (+37/+117)yjcE/yjcFputative transporter YjcE/pentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base A (52/48);  new base T (5/2);  total (57/50)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 110 total aligned reads displayed.

TGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAAT  >  NC_000913/4281518‑4281815
                                                                                                                                                    |                                                                                                                                                     
tGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGaaa                                                                                                                                                      >  1:440384/1‑150 (MQ=255)
 gCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGaaaa                                                                                                                                                     <  1:1422079/150‑1 (MQ=255)
   tGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAg                                                                                                                                                   <  1:912740/150‑1 (MQ=255)
    gAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAgg                                                                                                                                                  >  2:403078/1‑150 (MQ=255)
       cTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGc                                                                                                                                               <  1:377094/150‑1 (MQ=255)
       cTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGc                                                                                                                                               >  2:760041/1‑150 (MQ=255)
             caccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtt                                                                                                                                         >  1:509223/1‑150 (MQ=255)
             caccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtt                                                                                                                                         >  1:834519/1‑150 (MQ=255)
               ccTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCtttt                                                                                                                                       <  1:243856/150‑1 (MQ=255)
                cTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTg                                                                                                                                      <  1:221552/150‑1 (MQ=255)
                 tGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGc                                                                                                                                     <  1:1446570/150‑1 (MQ=255)
                  gcgcgcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCa                                                                                                                                    >  2:768861/1‑150 (MQ=255)
                      gcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTTCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGACGAAAGCGACGTAACCAGGCCCGAGCTATCCTGTTAAAGGTGCAAAAGATTGCACCtt                                                                                                                                <  2:445958/150‑1 (MQ=38)
                      gcCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCtt                                                                                                                                >  2:1390661/1‑150 (MQ=255)
                        cacaCGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTa                                                                                                                              >  2:929969/1‑150 (MQ=255)
                        cacaCGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTa                                                                                                                              <  1:1429501/150‑1 (MQ=255)
                        cacaCGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCtttt                                                                                                                              <  2:104355/150‑2 (MQ=255)
                        cacaCGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAAAAGATTGCACCtttt                                                                                                                              <  2:109906/150‑2 (MQ=255)
                          cacGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAAc                                                                                                                            >  1:436502/1‑150 (MQ=255)
                           acGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAACGGTGCAATCTTTTGCACCTTTAACa                                                                                                                           <  1:74911/150‑1 (MQ=255)
                               gagaTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGtt                                                                                                                       <  1:358041/150‑1 (MQ=255)
                                 gaTCAGTAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAg                                                                                                                     <  2:529569/150‑1 (MQ=255)
                                 gaTCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAg                                                                                                                     >  1:434718/1‑150 (MQ=255)
                                   tCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCt                                                                                                                   >  2:1188253/1‑150 (MQ=255)
                                   tCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGAGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTATTGCACCTTTAACAGGTTAGCt                                                                                                                   >  1:165796/1‑150 (MQ=255)
                                    cAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTc                                                                                                                  <  2:436502/150‑1 (MQ=255)
                                    cAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTc                                                                                                                  >  2:960414/1‑150 (MQ=255)
                                     aGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCg                                                                                                                 >  2:1064840/1‑150 (MQ=255)
                                     aGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCg                                                                                                                 <  2:685104/150‑1 (MQ=25)
                                      gCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCgg                                                                                                                >  2:1208344/1‑150 (MQ=255)
                                                   gCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATg                                                                                                   <  1:729185/150‑1 (MQ=255)
                                                    cAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGa                                                                                                  >  1:35159/1‑150 (MQ=255)
                                                       aaaTTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTc                                                                                               <  2:5715/150‑1 (MQ=255)
                                                           tACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGACGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGAGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaaa                                                                                           >  1:259090/1‑150 (MQ=255)
                                                            aCTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTaaaa                                                                                          >  1:346617/1‑150 (MQ=255)
                                                             cTGCTCGATCTCGATTTGCTTGAAGGGTTGCTAATTGAGTAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCTAGGTAATCTGAAAAATGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATTATTCTAAAAc                                                                                         <  2:259090/150‑1 (MQ=255)
                                                             cTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAAc                                                                                         >  1:636143/1‑150 (MQ=255)
                                                             cTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCGCGGGCTAACAATCATGCTGCTAAAAc                                                                                         >  1:687706/1‑150 (MQ=255)
                                                               gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                       >  1:1293280/1‑150 (MQ=255)
                                                               gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                       <  1:835867/150‑1 (MQ=255)
                                                               gCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaa                                                                                       >  2:260042/1‑150 (MQ=255)
                                                                cACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaa                                                                                      >  1:591017/1‑150 (MQ=255)
                                                                   gATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataa                                                                                   <  1:823031/150‑1 (MQ=255)
                                                                    aTCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaa                                                                                  >  1:366011/1‑150 (MQ=255)
                                                                    aTCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaa                                                                                  <  1:1250748/150‑1 (MQ=255)
                                                                    aTCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACATTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaa                                                                                  >  2:48532/1‑150 (MQ=255)
                                                                     tctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGATGTAATCTGAAAAAGGTGCAATCTTTTTCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaat                                                                                 <  1:93265/150‑1 (MQ=255)
                                                                     tctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaat                                                                                 <  2:1099032/150‑1 (MQ=255)
                                                                     tctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACaaataaat                                                                                 >  1:566710/1‑150 (MQ=255)
                                                                      ctcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAAtt                                                                                >  2:695464/1‑150 (MQ=255)
                                                                       tcGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTa                                                                               <  2:1207643/150‑1 (MQ=255)
                                                                              gCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATctctct                                                                        >  1:380503/1‑150 (MQ=255)
                                                                                 tGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtttt                                                                     >  2:366422/1‑150 (MQ=255)
                                                                                 tGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCtttt                                                                     >  1:460826/1‑150 (MQ=255)
                                                                                     gCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAg                                                                 <  1:865189/150‑1 (MQ=255)
                                                                                        ttGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCa                                                              <  2:1100640/150‑1 (MQ=255)
                                                                                          gCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAAt                                                            <  2:834519/150‑1 (MQ=255)
                                                                                          gCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAAt                                                            <  2:829029/150‑1 (MQ=255)
                                                                                          gCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAAtt                                                           <  1:978068/150‑1 (MQ=255)
                                                                                                 gAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTgg                                                     <  2:1070678/150‑1 (MQ=255)
                                                                                                   ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTa                                                                               <  1:1016612/122‑1 (MQ=255)
                                                                                                   ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTa                                                                               >  2:1016612/1‑122 (MQ=255)
                                                                                                   ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                   >  2:319476/1‑150 (MQ=255)
                                                                                                   ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                   >  2:113583/1‑150 (MQ=255)
                                                                                                   ggAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCt                                                   >  1:443122/1‑150 (MQ=255)
                                                                                                    gAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTg                                                  >  2:942683/1‑150 (MQ=255)
                                                                                                        aTCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGttct                                              <  1:1078713/150‑1 (MQ=255)
                                                                                                           aGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCg                                           >  1:896036/1‑150 (MQ=255)
                                                                                                            gTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGa                                          <  2:591017/150‑1 (MQ=255)
                                                                                                             tAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGAc                                         <  1:1064840/150‑1 (MQ=255)
                                                                                                             tAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGAc                                         <  1:768861/150‑1 (MQ=255)
                                                                                                              aaGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTTAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACt                                        <  1:45842/150‑1 (MQ=255)
                                                                                                                  cGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAc                                    <  1:1321775/150‑1 (MQ=255)
                                                                                                                  cGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAc                                    >  2:1308449/1‑150 (MQ=255)
                                                                                                                  cGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAc                                    <  1:946231/150‑1 (MQ=255)
                                                                                                                   gAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAAcc                                   >  2:609707/1‑150 (MQ=255)
                                                                                                                             tAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCa                         >  1:171114/1‑150 (MQ=255)
                                                                                                                             tAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCa                         <  1:609707/150‑1 (MQ=255)
                                                                                                                               aCCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATc                       >  2:532225/1‑150 (MQ=255)
                                                                                                                                ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCt                      <  2:165796/150‑1 (MQ=255)
                                                                                                                                ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCt                      >  1:1374459/1‑150 (MQ=255)
                                                                                                                                ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCt                      >  2:758702/1‑150 (MQ=255)
                                                                                                                                ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTATCTTTTaaa                                                                  >  2:1534805/1‑106 (MQ=255)
                                                                                                                                ccAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTATCTTTTaaa                                                                  <  1:1534805/106‑1 (MQ=255)
                                                                                                                                  tagCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTg                    >  1:749849/3‑150 (MQ=21)
                                                                                                                                  tagCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTg                    >  1:1082613/3‑150 (MQ=21)
                                                                                                                                   ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                   >  1:1172920/1‑150 (MQ=255)
                                                                                                                                   ggCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                   <  1:642296/150‑1 (MQ=255)
                                                                                                                                   agCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTGTTAAAGTCAATTAATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGc                   >  1:980240/2‑150 (MQ=21)
                                                                                                                                    gCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCa                  >  1:283154/1‑150 (MQ=255)
                                                                                                                                    gCCCGAGCTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCa                  >  2:155001/1‑150 (MQ=21)
                                                                                                                                      ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                <  2:1374459/150‑1 (MQ=255)
                                                                                                                                      ccGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCata                <  1:1390661/150‑1 (MQ=255)
                                                                                                                                         aGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACGTCAATCTGCATCTTGCataata             >  1:1043227/1‑150 (MQ=255)
                                                                                                                                            tAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCt          >  1:1481396/1‑150 (MQ=255)
                                                                                                                                             aaTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTg         <  1:532225/150‑1 (MQ=255)
                                                                                                                                              aTCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGc        <  1:760041/150‑1 (MQ=255)
                                                                                                                                               tCTGAAAAAGGTGCAATCTTTTGCACGTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTTTCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCa       <  2:509788/150‑1 (MQ=255)
                                                                                                                                                 tGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGt     >  1:1136667/1‑150 (MQ=255)
                                                                                                                                                    aaaaGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAAt  >  1:1369893/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCTGAACTTTACCACCTGCGCGCCACACGGGAGATCAGCAACGAAACGCTGCAAAAATTACTGCACGATCTCGATTTGCTTGAAGCGTTGCTAATTGAGGAAAATCAGTAAGCCGAAAGCGACGTAACCAGGCCCGAGGTAATCTGAAAAAGGTGCAATCTTTTGCACCTTTAACAGGTTAGCTCGGGCTAACATTCATGATTCTAAAACAAATAAATTATCTCTCTTTTAAAGTCAATTCATTGGCTGTTCTTCGACTTAACCTCAATCTGCATCTTGCATAATATCTGCAGTAAT  >  NC_000913/4281518‑4281815

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