Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,038,281 G→C 9.3% A181G (GCC→GGC)  dsbC ← protein disulfide isomerase DsbC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,038,2810GC 9.3% 308.2 / 18.1 97A181G (GCC→GGC) dsbCprotein disulfide isomerase DsbC
Reads supporting (aligned to +/- strand):  ref base G (37/51);  new base C (4/5);  total (41/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 = 8.74e-01

Only 100 of 104 total aligned reads displayed.

TTTGGTGTTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATCCCCAGCGCGTT  >  NC_000913/3038132‑3038429
                                                                                                                                                     |                                                                                                                                                    
tttGGTGTTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAActgg                                                                                                                                                      <  1:498459/150‑1 (MQ=255)
 ttGGTGTTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAActggc                                                                                                                                                     <  1:406855/150‑1 (MQ=255)
 ttGGTGTTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAActggc                                                                                                                                                     <  2:1153/150‑1 (MQ=255)
    gtgtTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAActggctgg                                                                                                                                                  <  1:1299523/150‑1 (MQ=255)
    gtgtTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAActggctgg                                                                                                                                                  >  1:1058530/1‑150 (MQ=255)
     tgtTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGt                                                                                                                                                 >  1:108072/1‑150 (MQ=255)
     tgtTAGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGt                                                                                                                                                 >  1:1075462/1‑150 (MQ=255)
      gtTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTg                                                                                                                                                <  2:1045263/150‑1 (MQ=255)
        tcgtcgAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCCCGTCGCAACTGGCTGGTGCg                                                                                                                                              >  2:1053137/1‑150 (MQ=255)
         cgtcgAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGa                                                                                                                                             >  2:588717/1‑150 (MQ=255)
           tcgAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCACCAGCCAGTTGCGACg                                                                                                                                           <  1:141302/150‑1 (MQ=255)
            cgAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCACCAGCCAGTTGCGACgt                                                                                                                                          <  1:1526744/150‑2 (MQ=255)
             gAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCt                                                                                                                                         <  2:1316857/150‑1 (MQ=255)
              aGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCtt                                                                                                                                        <  2:483250/150‑1 (MQ=255)
                gAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCACCAGCCAGTTGCGACgtggat                                                                                                                                      <  2:796864/150‑6 (MQ=255)
                gAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCtttt                                                                                                                                      <  1:1243289/150‑1 (MQ=255)
                 aaTTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTa                                                                                                                                     >  2:695693/1‑150 (MQ=255)
                 aaTTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTa                                                                                                                                     >  1:1366250/1‑150 (MQ=255)
                   ttCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTAcc                                                                                                                                   <  1:291899/150‑1 (MQ=255)
                        tCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTTGCAACTGGCTGGTGCGACGCTTTTACCTGCca                                                                                                                              <  1:137702/150‑1 (MQ=255)
                        tCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCca                                                                                                                              <  1:164142/150‑1 (MQ=255)
                         cATCTCTTTCGGCGGCTGGTAAACCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCcat                                                                                                                             >  1:1233304/1‑150 (MQ=255)
                          aTCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCGAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCcatc                                                                                                                            >  1:1403005/1‑150 (MQ=255)
                           tctcTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCcatca                                                                                                                           >  1:1305372/1‑150 (MQ=255)
                            ctctTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCcatcac                                                                                                                          >  2:1191988/1‑150 (MQ=255)
                             tctTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCcatcaca                                                                                                                         <  1:325164/150‑1 (MQ=255)
                                   ggcggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                   <  1:885916/150‑1 (MQ=255)
                                   ggcggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                   <  1:107184/150‑1 (MQ=255)
                                   ggcggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                   >  1:24170/1‑150 (MQ=255)
                                     cggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                   <  2:107989/148‑1 (MQ=255)
                                     cggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                   >  1:107989/1‑148 (MQ=255)
                                      ggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAAc                                                                                                                <  2:934434/150‑1 (MQ=255)
                                       gcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACATGCCATCACATCATCAAACg                                                                                                               >  1:323680/1‑150 (MQ=255)
                                        cTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGc                                                                                                              >  2:583615/1‑150 (MQ=255)
                                            tAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCACCAGCCAGTTGCGACgtgg                                                                                                                                        >  1:154369/1‑117 (MQ=255)
                                            tAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTg                                                                                                          <  2:1014177/150‑1 (MQ=255)
                                              aCCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGtt                                                                                                        >  1:513746/1‑150 (MQ=255)
                                               cccGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGttt                                                                                                       >  1:86646/1‑150 (MQ=255)
                                               cccGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCGTTGttt                                                                                                       >  1:356808/1‑150 (MQ=255)
                                                ccGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGtttt                                                                                                      <  2:1075462/150‑1 (MQ=255)
                                                 cGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGttttt                                                                                                     <  1:588717/150‑1 (MQ=255)
                                                  ggAACAAGTGTGCCATTGCTCAGCACAAGTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTa                                                                                                    <  1:495680/150‑1 (MQ=255)
                                                  ggAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTa                                                                                                    <  2:543759/150‑1 (MQ=255)
                                                  ggAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTa                                                                                                    <  1:684923/150‑1 (MQ=255)
                                                   gAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTAt                                                                                                   <  1:683832/150‑1 (MQ=255)
                                                   gAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTAt                                                                                                   <  2:323680/150‑1 (MQ=255)
                                                    aaCAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATc                                                                                                  <  2:1319243/150‑1 (MQ=255)
                                                    aaCAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATc                                                                                                  >  2:1298491/1‑150 (MQ=255)
                                                       aaGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCttt                                                                                               <  1:1013535/150‑1 (MQ=255)
                                                        aGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTAGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTc                                                                                              <  1:66840/150‑1 (MQ=255)
                                                          tgtgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGc                                                                                            <  1:1297888/150‑1 (MQ=255)
                                                          tgtgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGc                                                                                            >  1:385507/1‑150 (MQ=255)
                                                          tgtgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGc                                                                                            <  2:677049/150‑1 (MQ=255)
                                                            tgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGcac                                                                                          <  2:1305372/150‑1 (MQ=255)
                                                            tgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGcac                                                                                          <  2:1305163/150‑1 (MQ=255)
                                                                  tGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGa                                                                                    >  2:575440/1‑150 (MQ=255)
                                                                   gCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGAt                                                                                   <  1:1439111/150‑1 (MQ=255)
                                                                   gCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGAt                                                                                   >  2:272694/1‑150 (MQ=255)
                                                                    cTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATa                                                                                  <  1:1475650/150‑1 (MQ=255)
                                                                       aGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCt                                                                               >  2:1484314/1‑150 (MQ=255)
                                                                       aGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCt                                                                               <  2:922867/150‑1 (MQ=255)
                                                                             aCTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtt                                                                         >  1:113599/1‑150 (MQ=255)
                                                                               tGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtttc                                                                       <  2:925575/150‑1 (MQ=255)
                                                                               tGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtttc                                                                       <  1:1030781/150‑1 (MQ=255)
                                                                                gCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTTCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtttct                                                                      <  2:1371781/150‑1 (MQ=255)
                                                                                 ccGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtttctt                                                                     <  2:108072/150‑1 (MQ=255)
                                                                                 ccGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtttctt                                                                     >  2:1321058/1‑150 (MQ=255)
                                                                                   ggAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtttctttc                                                                   <  1:98753/150‑1 (MQ=255)
                                                                                   ggAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCAtttctttc                                                                   <  1:204615/150‑1 (MQ=255)
                                                                                      gTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTctct                                                                <  1:159227/150‑1 (MQ=255)
                                                                                         ccGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCa                                                             >  1:176427/1‑150 (MQ=255)
                                                                                                cGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGt                                                      >  1:1046604/1‑150 (MQ=255)
                                                                                                  ccAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTcc                                                    >  2:34534/1‑150 (MQ=255)
                                                                                                  ccAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTcc                                                    >  2:983105/1‑150 (MQ=255)
                                                                                                  ccAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTcc                                                    <  1:1191988/150‑1 (MQ=255)
                                                                                                        tGGACGCCGAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGcccc                                              <  2:1403005/150‑1 (MQ=255)
                                                                                                          gACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTg                                            <  2:24170/150‑1 (MQ=255)
                                                                                                            cGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCACCCGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGgc                                          >  2:1144367/1‑150 (MQ=255)
                                                                                                                   tGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGaa                                   <  1:2228/150‑1 (MQ=255)
                                                                                                                   tGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGaa                                   <  1:1388574/150‑1 (MQ=255)
                                                                                                                     cGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAg                                 <  1:1484314/150‑1 (MQ=255)
                                                                                                                       tAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCa                               <  1:272694/150‑1 (MQ=255)
                                                                                                                         aTGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCGATCACATCATCTAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAg                             <  2:1366250/150‑1 (MQ=255)
                                                                                                                         aTGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAg                             >  1:812856/1‑150 (MQ=255)
                                                                                                                                gCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGc                      >  2:599803/1‑150 (MQ=255)
                                                                                                                                gCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGc                      >  1:1049728/1‑150 (MQ=255)
                                                                                                                                 cAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCa                     <  1:361642/150‑1 (MQ=255)
                                                                                                                                 cAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTACCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCa                     <  1:76434/150‑1 (MQ=255)
                                                                                                                                     aaaagcGTCGCACCAGCCAGTTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGt                 >  2:830875/6‑150 (MQ=255)
                                                                                                                                      tCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTg                >  2:97514/1‑150 (MQ=255)
                                                                                                                                      aaagcGTCGCACCAGCCAGTTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTg                >  2:829841/5‑150 (MQ=255)
                                                                                                                                       ccACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGa               >  1:201540/1‑150 (MQ=255)
                                                                                                                                       ccACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATAATCAAACGCTTTGTTTTTATCTTTAGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGa               >  2:784496/1‑150 (MQ=255)
                                                                                                                                       aagcGTCGCACCAGCCAGTTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGa               >  1:381965/4‑150 (MQ=255)
                                                                                                                                          cGTCGCACCAGCCAGTTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATcc            <  2:381965/150‑1 (MQ=255)
                                                                                                                                           gTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATccc           >  2:828949/1‑150 (MQ=255)
                                                                                                                                            tCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATcccc          <  2:86646/150‑1 (MQ=255)
                                                                                                                                            tCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATAGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCTCGGTGATaccc          <  2:356808/150‑1 (MQ=255)
                                                                                                                                                    ggctggTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATCCCCAGCGCGtt  >  1:1423606/1‑150 (MQ=255)
                                                                                                                                                    ggctggTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATCCCCAGCGCGtt  >  2:1286340/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTGGTGTTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATCCCCAGCGCGTT  >  NC_000913/3038132‑3038429

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: