Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,038,283 T→A 9.3% P180P (CCA→CCT) ‡ dsbC ← protein disulfide isomerase DsbC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,038,2830TA 9.3% 277.6 / 15.9 97P180P (CCA→CCT) ‡dsbCprotein disulfide isomerase DsbC
Reads supporting (aligned to +/- strand):  ref base T (37/51);  new base A (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 = 2.97e-01

Only 100 of 102 total aligned reads displayed.

GTGTTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATCCCCAGCGCGTTGT  >  NC_000913/3038136‑3038431
                                                                                                                                                   |                                                                                                                                                    
gtgtTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAActggctgg                                                                                                                                                    >  1:1058530/1‑150 (MQ=255)
gtgtTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAActggctgg                                                                                                                                                    <  1:1299523/150‑1 (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                                                                                                                                       >  1:1366250/1‑150 (MQ=255)
             aaTTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTa                                                                                                                                       >  2:695693/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)
                           tttCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCcatcacatc                                                                                                                         <  1:1486179/150‑1 (MQ=255)
                               ggcggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                     <  1:885916/150‑1 (MQ=255)
                               ggcggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                     >  1:24170/1‑150 (MQ=255)
                               ggcggcTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACAtcatca                                                                                                                     <  1:107184/150‑1 (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                                                                                                                                          <  2:154369/120‑4 (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:385507/1‑150 (MQ=255)
                                                      tgtgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGc                                                                                              <  2:677049/150‑1 (MQ=255)
                                                      tgtgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGc                                                                                              <  1:1297888/150‑1 (MQ=255)
                                                        tgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGcac                                                                                            <  2:1305163/150‑1 (MQ=255)
                                                        tgCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGcac                                                                                            <  2:1305372/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:1388574/150‑1 (MQ=255)
                                                                                                               tGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGaa                                     <  1:2228/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                        >  1:1049728/1‑150 (MQ=255)
                                                                                                                            gCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGc                        >  2:599803/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)
                                                                                                                                                  ctggTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATCCCCAGCGCGTTGt  >  1:472433/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTGTTCGTCGAGGAATTCTTTCATCTCTTTCGGCGGCTGGTAACCCGGAACAAGTGTGCCATTGCTCAGCACAACTGCCGGAGTACCGCTAACGCCAAGCTGGACGCCAAGTGCGTAATGGTCGGCAATATCCACGTCGCAACTGGCTGGTGCGACGCTTTTACCTGCCATCACATCATCAAACGCTTTGTTTTTATCTTTCGCACACCAGATAGCTTTCATTTCTTTCTCTGCATCGCTGTCCAGCCCCTGGCGCGGGAAAGCAAGATAACGCACGGTGATCCCCAGCGCGTTGT  >  NC_000913/3038136‑3038431

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: