Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,507 (C)7→8 intergenic (+152/‑203) rnlB → / → yfjP CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,5001.C95.6% 391.0 / 12.0 117intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base C (52/56);  total (57/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 131 total aligned reads displayed.

TGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  NC_000913/2767402‑2767594
                                                                                                    |                                                                                                
tGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑c                                                                                                <  2:2324301/100‑1 (MQ=255)
  atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑ccc                                                                                              <  2:2252083/100‑1 (MQ=255)
  atatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑ccc                                                                                              <  2:3008624/100‑1 (MQ=255)
   tatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTGCCGGCTCAGGAAAACGTCTTCCCCAGTACACCCCCGACCGTATTT‑‑cccc                                                                                             >  2:1260529/1‑100 (MQ=255)
    atTAACTTGTCCAATATGATCAAATAGCATTAACTCCCCCTCACAATATACTGCATAAGGAACACATATTCCCCTGTGCACCCACGACTAAATTT‑‑ccc                                                                                              >  1:2223171/1‑98 (MQ=255)
     ttAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccccc                                                                                           <  2:971509/100‑1 (MQ=255)
     ttAACTTGTCCAATATGATCAAATAGCATTAACACCCGCTCAAAGCGTCGTGCAGAGGAGCCACGTTATCGACAGAGCACCAAAGACTAAAGTT‑‑ccacca                                                                                           >  2:1167476/1‑99 (MQ=255)
       aaCTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑cccccccc                                                                                         <  2:2995909/100‑2 (MQ=255)
           tGTCCAATATGATCAAATAGCATTAACCCCCCCTCCCAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCctt                                                                                      <  1:70545/98‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                  <  2:2514851/100‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAAc                                                                                  <  2:417735/100‑1 (MQ=255)
              ccAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACATTTTCCCCTGTGCACAAACGACTACATAT‑‑cccaaatcccaa                                                                                     >  2:2095921/1‑88 (MQ=255)
              ccAATATGATCAAATAGCATTAACACCCCCTCCCAAGTCACTGCATAGGACTCCCCCTTTCCCCTGTGCATACACAACTAACTTTACCTTCCCTTTaact                                                                                  >  1:2262983/1‑98 (MQ=255)
                   atGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTtctc                                                                             <  1:1013660/100‑1 (MQ=255)
                    tGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAAGCAACACGTTTCCCCCTGTGCATCCACGCCCAAATTTCCCCCCCCTTTAAACTtctc                                                                             >  1:2822159/1‑99 (MQ=255)
                     gATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaa                                                                           <  2:2634571/100‑1 (MQ=255)
                     gATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCaa                                                                           <  1:3058832/100‑1 (MQ=255)
                           aTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGt                                                                     <  2:296831/100‑1 (MQ=255)
                            tAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTc                                                                    <  2:3039757/100‑1 (MQ=255)
                              gCATTAACCCCCCCTTACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCAc                                                                  <  2:3009636/100‑1 (MQ=255)
                                aTTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCa                                                                <  2:1155251/100‑1 (MQ=255)
                                     cccccgcTGACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTTAACTTCTCAATTACGGCACCAAtttt                                                           >  1:310509/1‑100 (MQ=255)
                                        ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  2:1870040/100‑1 (MQ=255)
                                        ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTg                                                        <  2:495383/100‑1 (MQ=255)
                                          ccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                      <  2:458724/100‑1 (MQ=255)
                                           cTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaa                                                     >  1:1616923/1‑100 (MQ=255)
                                            tCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGaaa                                                    >  2:87715/1‑100 (MQ=255)
                                             cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGGACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                   >  1:1640276/1‑100 (MQ=255)
                                              acaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATc                                                  <  2:650520/100‑1 (MQ=255)
                                               caACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATcc                                                 <  2:551181/100‑1 (MQ=255)
                                                aaCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTTAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCt                                                <  1:18277/100‑1 (MQ=255)
                                                 aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTTAATTTCCCCCCCCTTTAAACTTCTCAATTATGTCACCAATTTTCTGGAAATCCTa                                               <  1:1193962/100‑1 (MQ=255)
                                                  cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                              >  1:2500297/1‑100 (MQ=255)
                                                  cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAACTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTa                                               <  2:2635767/99‑1 (MQ=255)
                                                   gatcttGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                              <  2:1796090/95‑1 (MQ=255)
                                                   gTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                             >  1:3011688/1‑100 (MQ=255)
                                                    tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatga                                            <  2:88025/100‑1 (MQ=255)
                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                          >  1:1683206/1‑100 (MQ=255)
                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                          >  1:541166/1‑100 (MQ=255)
                                                       cTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg                                             <  1:1381274/97‑1 (MQ=255)
                                                       cTGCATAGCGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg                                          <  1:2808878/100‑1 (MQ=255)
                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         >  1:391664/1‑100 (MQ=255)
                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         <  2:392909/100‑1 (MQ=255)
                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         >  2:36802/1‑100 (MQ=255)
                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         <  1:2951075/100‑1 (MQ=255)
                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         <  2:333752/100‑1 (MQ=255)
                                                        tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTACCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGc                                         <  1:2951846/100‑1 (MQ=255)
                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        <  2:268694/100‑1 (MQ=255)
                                                         gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                        <  2:268972/100‑1 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  1:457907/1‑100 (MQ=255)
                                                            tAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCttt                                     >  2:929004/1‑100 (MQ=255)
                                                             aGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                    >  1:1060917/1‑100 (MQ=255)
                                                             aGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                    >  2:1564978/1‑100 (MQ=255)
                                                             aGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                    >  1:1802376/1‑100 (MQ=255)
                                                              gggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                   >  2:1978522/1‑100 (MQ=255)
                                                              gggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                   >  2:2579551/1‑100 (MQ=255)
                                                              gggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                   >  2:137128/1‑100 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  >  1:128763/1‑100 (MQ=255)
                                                               ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  <  1:927559/100‑1 (MQ=255)
                                                                gAACACGTTTTCCCCTGTGCACCCACGACTAAATTGCCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTccc                                  <  1:517955/100‑1 (MQ=255)
                                                                gAACACGTTTGCCCCTGTGCACCCACGACTAATTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                 <  2:778763/100‑1 (MQ=255)
                                                                 aaCACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                                <  2:2670760/100‑1 (MQ=255)
                                                                  acacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                               >  2:1777962/1‑100 (MQ=255)
                                                                   cacGTTTTCCCCTGTGCACCCACTACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                              <  1:2920367/100‑1 (MQ=255)
                                                                   cacGTTTTCCCCTGTGCACCCACGACTAGATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                              >  2:466556/1‑100 (MQ=255)
                                                                   cacGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCgg                                >  1:2674262/1‑98 (MQ=255)
                                                                     cattttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                            >  2:2063977/3‑100 (MQ=255)
                                                                     cGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaag                            >  1:174510/1‑99 (MQ=255)
                                                                     cGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaaa                            >  1:1964439/1‑100 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCTCGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           <  2:2095573/100‑1 (MQ=255)
                                                                      gTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                           >  1:689133/1‑100 (MQ=255)
                                                                       ttttCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCATTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                          <  2:2151438/100‑1 (MQ=255)
                                                                         ttCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGACGCCTTTCCCGGTAAATGGt                        <  2:78021/100‑1 (MQ=255)
                                                                          tCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTg                       <  2:1693773/100‑1 (MQ=255)
                                                                           ccccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGc                      <  1:2551468/100‑1 (MQ=255)
                                                                             ccTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                    >  2:2235934/1‑100 (MQ=255)
                                                                                gtgCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGa                 >  1:1101470/1‑100 (MQ=255)
                                                                                  gCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGtt              >  2:2629368/1‑100 (MQ=255)
                                                                                   cACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGa                 >  1:265540/1‑97 (MQ=255)
                                                                                   cACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAgn               <  2:1385477/99‑2 (MQ=255)
                                                                                    aCCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTc             <  2:170353/100‑1 (MQ=255)
                                                                                     cccATGATTATATTTCCCCCCCCTCTTAACTTCTCAATAAGGTCACCATTTTTCTGGAAATACTATGATGCCTTTCCCGGTAAATGGTACTGTGAGTTCa            <  1:325904/100‑1 (MQ=255)
                                                                                     cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa            >  2:2524730/1‑100 (MQ=255)
                                                                                     cccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCa            >  2:1270257/1‑100 (MQ=255)
                                                                                      ccACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAg           >  2:1237888/1‑100 (MQ=255)
                                                                                       cACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa          >  2:2997850/1‑100 (MQ=255)
                                                                                       cACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGa          >  2:1984148/1‑100 (MQ=255)
                                                                                         cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        >  1:90570/1‑100 (MQ=255)
                                                                                         cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        >  1:1251363/1‑100 (MQ=255)
                                                                                         cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        >  1:2791410/1‑100 (MQ=255)
                                                                                         cGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa        >  1:2264792/1‑100 (MQ=255)
                                                                                          gACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg       >  2:2842362/1‑100 (MQ=255)
                                                                                           aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAg       <  1:1739489/99‑1 (MQ=255)
                                                                                           aCTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa      <  2:312848/100‑1 (MQ=255)
                                                                                            cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     <  1:2063969/100‑1 (MQ=255)
                                                                                            cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     >  1:2475646/1‑100 (MQ=255)
                                                                                            cTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt     >  1:2593711/1‑100 (MQ=255)
                                                                                              aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc   >  2:1818523/1‑100 (MQ=255)
                                                                                              aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc   >  1:1732280/1‑100 (MQ=255)
                                                                                              aaaTTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTc   >  2:1415428/1‑100 (MQ=255)
                                                                                                    |                                                                                                
TGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTT‑‑CCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCT  >  NC_000913/2767402‑2767594

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