Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,907,241 G→C 5.4% N65K (AAC→AAG cspC ← stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,907,2410GC 5.4% 357.2 / 9.4 112N65K (AAC→AAGcspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base G (55/51);  new base C (4/2);  total (59/53)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

Only 100 of 115 total aligned reads displayed.

GTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACC  >  NC_000913/1907094‑1907390
                                                                                                                                                   |                                                                                                                                                     
gtAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGtt                                                                                                                                                     >  1:1461418/1‑150 (MQ=255)
gtAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGtt                                                                                                                                                     >  1:739207/1‑150 (MQ=255)
   aCGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAAc                                                                                                                                                  >  2:1290330/1‑150 (MQ=255)
    cgcgAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACa                                                                                                                                                 <  2:1634447/150‑1 (MQ=255)
    cgcgAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTAGTCACACTTCAGATCAGTGGATTCGATCTGATAGCTGTTACGTTaaga                                                                                                                                                 >  2:83637/1‑148 (MQ=255)
     gcgAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAg                                                                                                                                                >  2:77605/1‑150 (MQ=255)
          tCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGACAGCTGTTAACGTAACAGCtatc                                                                                                                                           <  2:564626/150‑4 (MQ=255)
            aCGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGACAGCTGTTAACGTAACAGCtatctg                                                                                                                                         <  2:614476/150‑6 (MQ=255)
             cgcgGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGa                                                                                                                                        >  1:511200/1‑150 (MQ=255)
             cgcgGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGa                                                                                                                                        >  1:560467/1‑150 (MQ=255)
                 gcggGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCtt                                                                                                                                    >  1:1013181/1‑150 (MQ=255)
                     gCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTg                                                                                                                                >  1:362028/1‑150 (MQ=255)
                      cAAATATCCGACCGCCAGGAGCGAAGATAAAGCATAGTACACAGCCTCGCAATCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACATCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTgg                                                                                                                               <  2:1059391/150‑1 (MQ=255)
                           atCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTc                                                                                                                          <  2:1151394/150‑1 (MQ=255)
                            tCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTcc                                                                                                                         <  1:1526613/150‑1 (MQ=255)
                             ccGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCt                                                                                                                        <  1:1260333/150‑1 (MQ=255)
                               gAACGCCAGGAGCGAAGATTAAGCATAGTAAAAAGCCTCGCATTCGCCGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCATGa                                                                                                                      >  1:924870/1‑150 (MQ=255)
                                aaCGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGaa                                                                                                                     <  1:1121993/150‑1 (MQ=255)
                                 aCGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTACTGAAt                                                                                                                    <  1:1710904/150‑1 (MQ=255)
                                  cGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTAACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAAtt                                                                                                                   >  2:1275631/1‑150 (MQ=255)
                                   gCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAAttt                                                                                                                  >  1:843143/1‑150 (MQ=255)
                                    ccAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTc                                                                                                                 >  2:885322/1‑150 (MQ=255)
                                     cAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCg                                                                                                                <  1:885322/150‑1 (MQ=255)
                                         aGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACt                                                                                                            >  1:1138873/1‑150 (MQ=255)
                                         aGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACt                                                                                                            <  2:444614/150‑1 (MQ=255)
                                         aGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACt                                                                                                            >  1:497408/1‑150 (MQ=255)
                                         aGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACt                                                                                                            >  1:1338931/1‑150 (MQ=255)
                                         aGCGAAGATAAAGAATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGCCCGTCCTGAATTTCGAACt                                                                                                            >  1:503998/1‑150 (MQ=255)
                                            gAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCaa                                                                                                         >  1:1285851/1‑150 (MQ=255)
                                              aGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACg                                                                                                       <  1:54435/150‑1 (MQ=255)
                                              aGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACg                                                                                                       <  1:19487/150‑1 (MQ=255)
                                                aTAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGtt                                                                                                     >  1:835424/1‑150 (MQ=255)
                                                 tAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTc                                                                                                    >  1:1156510/1‑150 (MQ=255)
                                                   aaGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTg                                                                                                  <  2:425360/150‑1 (MQ=255)
                                                   aaGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTg                                                                                                  <  2:667420/150‑1 (MQ=255)
                                                     gCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGAc                                                                                                <  1:537836/150‑1 (MQ=255)
                                                     gCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGAc                                                                                                >  1:788369/1‑150 (MQ=255)
                                                      cATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGAcc                                                                                               <  1:1459489/150‑1 (MQ=255)
                                                       aTAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCt                                                                                              >  1:1614076/1‑150 (MQ=255)
                                                         aGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGTACTCAACGTTCTGACCTTc                                                                                            <  2:369116/150‑1 (MQ=255)
                                                         aGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTc                                                                                            >  1:822291/1‑150 (MQ=255)
                                                         aGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTc                                                                                            <  1:1226603/150‑1 (MQ=255)
                                                          gTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCa                                                                                           >  1:795221/1‑150 (MQ=255)
                                                           tAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAg                                                                                          <  1:1353323/150‑1 (MQ=255)
                                                           tAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAg                                                                                          <  2:1056877/150‑1 (MQ=255)
                                                            aaaaaGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGc                                                                                         <  2:1561865/150‑1 (MQ=255)
                                                             aaaaGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGcc                                                                                        >  1:771817/1‑150 (MQ=255)
                                                                       cATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGa                                                                              >  2:194424/1‑150 (MQ=255)
                                                                        cgaatgcgaGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGaa                                                                             <  2:896588/144‑1 (MQ=255)
                                                                        aTTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGaa                                                                             <  2:23326/150‑1 (MQ=255)
                                                                        aTTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGaa                                                                             >  1:559901/1‑150 (MQ=255)
                                                                         ttCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAg                                                                            <  1:77605/150‑1 (MQ=255)
                                                                         ttCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAg                                                                            <  2:1587303/150‑1 (MQ=255)
                                                                            gacgaGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCa                                                                         >  2:640651/1‑150 (MQ=255)
                                                                                aGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTGTTGAAGCCATTAc                                                                     >  1:1585740/1‑150 (MQ=255)
                                                                                 ggCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTAcc                                                                    >  1:1067016/1‑150 (MQ=255)
                                                                                   cTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGATCTCAACGTTCTGACCTTCAGCCAGCGTTTTGAAGCCAGTACCCt                                                                  >  1:104169/1‑150 (MQ=255)
                                                                                    tttATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTg                                                                 <  2:244603/150‑1 (MQ=255)
                                                                                      tataGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGCCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGa                                                               >  1:284019/1‑150 (MQ=255)
                                                                                       ataGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCGGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCACGCTTGTGACCTTCAGCCCGAGTTTTGTAGCGATTACCCTgggt                                                              >  2:1069657/1‑148 (MQ=255)
                                                                                         aGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAg                                                            <  1:1264866/150‑1 (MQ=255)
                                                                                           cGAGATACAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGTCGTCCTGAATTTCGAACTCCCCGTTCTGACCTTCAGCCAGAGTGTTGAAGCCATTACCCTGGATAGCg                                                          >  2:31959/1‑150 (MQ=255)
                                                                                            gagaTTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCgg                                                         >  1:724804/1‑150 (MQ=255)
                                                                                              gaTTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGgag                                                       <  2:1013181/150‑1 (MQ=255)
                                                                                                  gAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGGTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGACAGAGTGTAGAAGCCATTTCACTGGATAGCGGAGAAgt                                                   >  1:296681/1‑150 (MQ=255)
                                                                                                   aaGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                  <  1:1712519/150‑1 (MQ=255)
                                                                                                   aaGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                  <  2:1461418/150‑1 (MQ=255)
                                                                                                   aaGCGTAATCAGACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtg                                                  <  2:87971/150‑1 (MQ=255)
                                                                                                    aGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAgtgt                                                 <  2:1463247/150‑1 (MQ=255)
                                                                                                    aGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCATTCTGGCCCTCCTGAAATTCGAAATCAACGATCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGCATAGCGGAGAAgtgt                                                 >  2:333314/1‑150 (MQ=255)
                                                                                                    aGCGTATTCACACTTCAGATCAGTGGATACGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCATGAATTTCGAACTCAACGTTCTGACCTTCCGCCAGAGTTTTGCAGCCATTACCCTGGATAGCGGAGAAgagt                                                 >  1:1548907/1‑150 (MQ=255)
                                                                                                      cGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTAc                                               <  2:446767/150‑1 (MQ=255)
                                                                                                       gTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACg                                              >  1:883238/1‑150 (MQ=255)
                                                                                                         aTTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGaa                                            >  2:1087913/1‑150 (MQ=255)
                                                                                                           tCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAaca                                          <  2:395399/150‑1 (MQ=255)
                                                                                                             acacTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGCCCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca                                        >  1:139470/1‑150 (MQ=255)
                                                                                                             acacTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAAcaca                                        <  2:511200/150‑1 (MQ=255)
                                                                                                                     gATCAGTGGATGCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTtgct                                <  1:83637/150‑1 (MQ=255)
                                                                                                                         aGTGGATTCGAACAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCa                            >  2:432747/1‑150 (MQ=255)
                                                                                                                            ggATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGCATTTCGAACTCAACGATCTGACCTTCAGCCAGAGTATTGAAGACATTACACTCGATAGCTGAGAAGAGTACGAACACATCGTTGCTGCCATCa                         >  2:1047894/1‑150 (MQ=255)
                                                                                                                                 cGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGCACACATCTTTGCTGCCATCAGCCgg                    >  1:736315/1‑150 (MQ=255)
                                                                                                                                 cGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCgg                    <  2:835424/150‑1 (MQ=255)
                                                                                                                                 cGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCgg                    >  2:209518/1‑150 (MQ=255)
                                                                                                                                  gATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGa                   >  1:378371/1‑150 (MQ=255)
                                                                                                                                   aTCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAg                  >  2:693377/1‑150 (MQ=255)
                                                                                                                                    tCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACGTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGATGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                 <  1:1275631/150‑1 (MQ=255)
                                                                                                                                     ccggcaGCTGTTAACGTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTTGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACAGCTTTGCTGCCATCAGCCGGAGtt                >  2:1261113/6‑149 (MQ=255)
                                                                                                                                      cggcaGCTGTTAACGTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa               >  1:656886/5‑150 (MQ=255)
                                                                                                                                      cggcaGCTGTTAACGTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa               >  2:645859/5‑150 (MQ=255)
                                                                                                                                      aGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa               >  1:1424079/1‑150 (MQ=255)
                                                                                                                                       ggcaGCTGTTAACGTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAAATCAACCTTCTGACCTTCAACCAGAGTTTTGAAGCCATTACACTGGATAGCGGAGAAGTGTACGAACACAACTTAGCTGCCATCAGCCGGAGTAAt              >  2:630024/4‑150 (MQ=255)
                                                                                                                                           gCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGCCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg          <  1:64241/150‑1 (MQ=255)
                                                                                                                                           gCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGATAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg          <  1:438767/150‑1 (MQ=255)
                                                                                                                                            cTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGc         <  1:1099583/150‑1 (MQ=255)
                                                                                                                                              gTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCa       <  2:338732/150‑1 (MQ=255)
                                                                                                                                                 aCGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg          <  1:350574/144‑1 (MQ=255)
                                                                                                                                                 aCGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg          >  2:350574/1‑144 (MQ=255)
                                                                                                                                                  cGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc   >  2:525251/1‑150 (MQ=255)
                                                                                                                                                  cGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAc   >  2:518702/1‑150 (MQ=255)
                                                                                                                                                   gTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAAcc  >  2:1294641/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAACGCGAATCACGCGGCGGGCAAATATCCGAACGCCAGGAGCGAAGATAAAGCATAGTAAAAAGCCTCGCATTCGACGAGGCTTTATAGCGAGATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACC  >  NC_000913/1907094‑1907390

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