Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,724,587 C→T G80G (GGG→GGA)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,5870CT95.1% 668.7 / NA 226G80G (GGG→GGA) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base C (0/11);  new base T (123/92);  total (123/103)
Fisher's exact test for biased strand distribution p-value = 1.30e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

GATCTCGTATTCTTTTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTT  >  CP006881/3724440‑3724732
                                                                                                                                                   |                                                                                                                                                 
 aTCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGGAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGGGCTCTCCt                                                                                                                                                >  1:2705637/1‑150 (MQ=255)
 aTCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCt                                                                                                                                                >  1:1695470/1‑150 (MQ=255)
  tctcATATTCTTTTTTAAAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCtt                                                                                                                                               <  2:1072241/150‑1 (MQ=255)
   ctcATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTc                                                                                                                                              <  1:136469/150‑1 (MQ=255)
    tcatATTCTTTTTTATTGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTcc                                                                                                                                             <  1:3049713/147‑1 (MQ=255)
     catATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAGTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCTGCGGCGCCGCGAGCTCTCCTTCCg                                                                                                                                            >  1:919505/3‑150 (MQ=255)
     catATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCg                                                                                                                                            <  2:3174484/148‑1 (MQ=255)
       tATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                                          >  2:957460/1‑150 (MQ=255)
        aTTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGGGCTCTCCTTCCGCAt                                                                                                                                         >  1:31206/1‑150 (MQ=255)
          tCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTCTTCTGCACAAGCACCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAg                                                                                                                                       <  1:406056/150‑1 (MQ=255)
          tCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGGGCCGCGAGCTCTCCTTCCGCATAg                                                                                                                                       >  1:2583939/1‑150 (MQ=255)
          tCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAg                                                                                                                                       >  1:1606542/1‑150 (MQ=255)
           cTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                                      <  2:1784251/150‑1 (MQ=255)
               tttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCtgat                                                                                                                                  >  1:2848330/1‑150 (MQ=255)
               tttATAGAGCGCCTTTACCTCATGCAGCAATTTATTCTGCACAAGCAGCCAATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCACTCCTTCCGCATAGCtgat                                                                                                                                  <  1:2329900/150‑1 (MQ=255)
                ttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCtgatg                                                                                                                                 >  1:1916557/1‑150 (MQ=255)
                ttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAAGGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGGtgatg                                                                                                                                 >  1:3773088/1‑150 (MQ=255)
                  ataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGGGCTCTCCTTCCGCATAGCTggtgat                                                                                                                               >  2:891307/1‑150 (MQ=255)
                  ataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCtgatggt                                                                                                                               >  1:2116708/1‑148 (MQ=255)
                  ataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTgatgat                                                                                                                               >  2:3752096/1‑150 (MQ=255)
                  ataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTgatgat                                                                                                                               <  1:1112136/150‑1 (MQ=255)
                    agagCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGGTGatat                                                                                                                             >  2:2809735/1‑150 (MQ=255)
                     gagCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCCCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATc                                                                                                                            <  2:1528025/150‑1 (MQ=255)
                       gcgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCt                                                                                                                          <  2:2652952/150‑1 (MQ=255)
                        cgtctttATCTCATGCAGCAATTTATTCTGCCCAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTc                                                                                                                         <  1:552214/146‑1 (MQ=255)
                        cgcCTTTATCTCCTGCAGCATTTTCTTCTGCCCAAGCAGCCCATTTTGAACCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTc                                                                                                                         <  2:838832/150‑1 (MQ=255)
                        cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGGAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTc                                                                                                                         >  1:323486/1‑150 (MQ=255)
                         gcctttATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCa                                                                                                                        >  2:2553092/4‑150 (MQ=255)
                          cctttATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                                                                                                       <  2:3316934/148‑1 (MQ=255)
                             ttATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                                                                                                    >  1:3180467/1‑150 (MQ=255)
                              tATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTAt                                                                                                                   >  2:350533/1‑150 (MQ=255)
                                    aTGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTat                                                                                                             <  2:808460/150‑1 (MQ=255)
                                     tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGGGCTCTCCTTCCGCATGGCTGTTGCAATCCTCAGCTATATCtata                                                                                                            >  2:2168525/1‑150 (MQ=255)
                                     tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCtata                                                                                                            <  2:2026450/150‑1 (MQ=255)
                                      gcagcaATTTATTCTGCACAAGCTGCCCATTTTGAATCCTTTCGATCGCAAACGACAAATGGTCGGAAAGCGCGATATGAATGTGGTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATAACCTCAGCTATTTctgtat                                                                                                           >  2:695649/1‑145 (MQ=255)
                                      gcagcaATTTATTCTGCACAAGCAGCCCATTTTGGATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTatat                                                                                                           >  1:3350804/1‑150 (MQ=255)
                                         gcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGACCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGct                                                                                                        <  2:2517965/150‑1 (MQ=255)
                                          caATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCCCTCCTTCCGCATAGCTGATGATATCCTCAGCTCTTTCTCTATGctc                                                                                                       >  1:2738813/1‑150 (MQ=255)
                                            aTTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGGGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGctcct                                                                                                     >  2:625558/1‑150 (MQ=255)
                                            aTTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATAGCCTTAGCTATTTCTATATGctcct                                                                                                     >  1:3512358/1‑150 (MQ=255)
                                             tttATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGGGCTCTCCTTCCGCATAGCCGATGATATCCTCACCCATTTCTATATGctcctc                                                                                                    >  1:1315348/1‑150 (MQ=255)
                                                aTTCTGCACAAGCAGACCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGGGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCACATAGACGATGATATCCTCACCTCTTTCTAGATGTTCCTCCgg                                                                                                 >  1:3310132/1‑150 (MQ=255)
                                                 ttCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGc                                                                                                >  1:3636512/1‑150 (MQ=255)
                                                    tGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTGGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAgt                                                                                             <  2:2793535/150‑1 (MQ=255)
                                                     gCACAAGCAGCCCATTTTGATTCCTTTCGATCGCAAAGGACAAAAGGTCGGAAAGCGCGATTTGAATGTGGTCGCTCAGCGGCGCCGCGAGCGCTCCTTCCGCATCGCGGGCGATGTCCTCAGCTATTTCTAAATGCTCCTCCGCCCgtg                                                                                            >  1:3384460/1‑150 (MQ=255)
                                                      cacaAGCAGCCCTTTTTGAATCCTTTGGATGGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGACCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAgtgt                                                                                           <  1:1151464/150‑1 (MQ=255)
                                                        caAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGttt                                                                                         <  2:1399504/150‑1 (MQ=255)
                                                             agcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATACCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGGTTGCAgg                                                                                    >  2:2464949/1‑150 (MQ=255)
                                                             agcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCCGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCACTGTTTGCAgg                                                                                    >  1:355418/1‑150 (MQ=255)
                                                               cccATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAggtt                                                                                  >  2:1545146/1‑148 (MQ=255)
                                                                ccATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGGTATTTCTATATGCTCCTCCCGCAGGGTTTGCAGGaat                                                                                 >  2:219217/1‑148 (MQ=255)
                                                                ccATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCTCATATCTGATGATATCTTAAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGAtt                                                                                 <  1:997965/150‑1 (MQ=255)
                                                                  aTTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCCCGAGCTCTCCTTCCGCATAGCTGGTGATATCCTCCGCTATTTCTATATGCTCCCCCGGCAGTGTTTGCCGGATTTg                                                                               >  1:3536232/1‑150 (MQ=255)
                                                                  aTTTTGAATCCTTCCGATCGCAAAGGCCAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTg                                                                               <  2:534532/150‑1 (MQ=255)
                                                                   ttttGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGt                                                                              >  2:1992072/1‑150 (MQ=255)
                                                                   ttttGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGACCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGt                                                                              <  2:2093976/150‑1 (MQ=255)
                                                                       gAATCCTTTCGACCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCCCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGtttga                                                                          <  1:438860/150‑3 (MQ=255)
                                                                         aTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGAATTGtttggac                                                                        >  2:1659415/1‑146 (MQ=255)
                                                                             tttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCCCCGGCAGTGTTTGCAGGATTTGTTTGAACtttt                                                                    <  1:2414807/150‑1 (MQ=255)
                                                                                cGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTTAt                                                                 >  1:3792494/1‑150 (MQ=255)
                                                                                cGATCGCAAACGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCAt                                                                 <  2:3162212/150‑1 (MQ=255)
                                                                                 gatgccAAAGGACAAATGGCCGGGAAGGGCGATATGAATGTGACCGCTCAGCGCCGCTGCGAGCTCTCCTTCCACATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCAtt                                                                <  2:2685901/145‑1 (MQ=255)
                                                                                  aTCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCCATTTCTATATGCTCCTCCGGCAGTGTGTGCAGGATGTGTTTGAACTTTTCATtc                                                               >  2:2908890/1‑150 (MQ=255)
                                                                                    cGCACAGGACAAATGGCCGGAAAGCGCGATCTGAATGTGACGGCCCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATtctc                                                             <  2:2452968/150‑1 (MQ=255)
                                                                                    cGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGGACTTTTCATtctc                                                             >  2:1722654/1‑150 (MQ=255)
                                                                                      cAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGCTGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGGGTTTGCAGGATTTGTTTGGACTTTTCATTCtcgt                                                           >  2:2744153/1‑150 (MQ=255)
                                                                                      cAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgt                                                           >  2:3943619/1‑150 (MQ=255)
                                                                                      cAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgt                                                           <  1:3782830/150‑1 (MQ=255)
                                                                                       aaaGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgtc                                                          >  1:3594030/1‑150 (MQ=255)
                                                                                           gACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCgcgc                                                      >  2:3775649/1‑150 (MQ=255)
                                                                                             cAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCAc                                                    <  2:1916557/150‑1 (MQ=255)
                                                                                              aaaTGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGGTATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGGTTGAACTTTTCATTCTCGTCGCGCACg                                                   >  2:415241/1‑150 (MQ=255)
                                                                                              aaaTGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGGTTGAACTTTTCATTCTCGTCGCGCACg                                                   >  2:3517098/1‑150 (MQ=255)
                                                                                               aaTGGTCGGAAAGCGCGATATGAATGTGATCGCTCTGCGGCGCCGCGAGCTCTCCTTGCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGa                                                  >  2:3075633/1‑150 (MQ=255)
                                                                                                 tGGTCGGAAAGCGCGATATGAATGTGACCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACa                                                <  1:3404546/150‑1 (MQ=255)
                                                                                                  ggTCTGAAAGAGCGATATGCTTTTGATCGCTCAGCGGCGCCGCATGCTCTCCTTCCGCATAGCTGATGATTTCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACaa                                               <  1:1995691/150‑1 (MQ=255)
                                                                                                      ggAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATTTCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGGACTTTTCATTCTCGTCGCGCCCGAcccagtt                                           >  2:3875546/2‑144 (MQ=255)
                                                                                                      ggAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATCCCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGAt                                           <  2:3699734/149‑1 (MQ=255)
                                                                                                            cgcgATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTc                                     >  2:665453/1‑150 (MQ=255)
                                                                                                             gcgATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTcc                                    >  2:224375/1‑150 (MQ=255)
                                                                                                                atatGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                 <  2:2848330/150‑1 (MQ=255)
                                                                                                                  atGATTGTGATCGCCCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACttt                               <  2:1177042/150‑1 (MQ=255)
                                                                                                                     aaTGTGATCGCCCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTAtt                            <  1:350533/150‑1 (MQ=255)
                                                                                                                       tgtgATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCa                          <  1:665453/150‑1 (MQ=255)
                                                                                                                         tgATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATTTGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGGCAAAGGTCTTTTCCACTTTATTCATa                        >  1:3704180/1‑150 (MQ=255)
                                                                                                                           aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGTACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAgg                      >  2:1588857/1‑150 (MQ=255)
                                                                                                                              gcccAGCGCCGCCCCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGAt                   <  2:2426780/147‑1 (MQ=255)
                                                                                                                               cccAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATg                  <  2:3137273/148‑1 (MQ=255)
                                                                                                                                  aGCGGCGCCGCGAGCTCTCCTTCCGCATAACTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAga               >  2:1042852/1‑149 (MQ=255)
                                                                                                                                   gcggcgCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAgat              <  2:2636046/150‑3 (MQ=255)
                                                                                                                                    cggcgCCGCGGGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATc             <  1:1087689/150‑1 (MQ=255)
                                                                                                                                     ggcgCCGCGGGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCCGCAGTGTTTGTAGGATTTGGTTGAACTTTTCATTCTCGTCCCGCACGACAAAGATCTTTTCCACTTTATTCATAGGgatcagatca            >  1:1507549/1‑150 (MQ=255)
                                                                                                                                        gccgcGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGACCTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCAttt         <  1:3943619/150‑1 (MQ=255)
                                                                                                                                        gccgcGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCAttt         >  2:3193980/1‑150 (MQ=255)
                                                                                                                                           gcgAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTc      <  1:2354787/150‑1 (MQ=255)
                                                                                                                                            cgAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGGGTTTGCAGGATTTGGTTGAACTTTTCATTCTCGTCGCGCACGGCAAAGATCTTTTCCACTTTATTCATAGGGATGAGAGCATTTTTCt     >  1:596471/1‑150 (MQ=255)
                                                                                                                                             gAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCtt    >  2:1447420/1‑150 (MQ=255)
                                                                                                                                               gCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCtttt  <  2:336955/150‑1 (MQ=255)
                                                                                                                                               gCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCtttt  <  2:1613827/150‑1 (MQ=255)
                                                                                                                                               gCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCtttt  <  2:294153/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATCTCGTATTCTTTTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTT  >  CP006881/3724440‑3724732

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑