Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,751,558 G→T 100% T134T (ACC→ACA yfjH ← CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,751,5580GT97.8% 883.9 / 8.0 278T134T (ACC→ACAyfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base T (131/141);  total (137/141)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 278 total aligned reads displayed.

ATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGT  >  CP009273/2751414‑2751701
                                                                                                                                                |                                                                                                                                               
 tGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTaa                                                                                                                                             <  1:2286873/148‑1 (MQ=255)
  gAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAg                                                                                                                                            <  2:209715/148‑1 (MQ=255)
  gAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAg                                                                                                                                            >  1:1519740/1‑148 (MQ=255)
  gAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAg                                                                                                                                            <  1:3147505/148‑1 (MQ=255)
   aaaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAgg                                                                                                                                           <  1:1380498/148‑1 (MQ=255)
   aaaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAgg                                                                                                                                           >  2:2334799/1‑148 (MQ=255)
   aaaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAgg                                                                                                                                           >  2:1331730/1‑148 (MQ=255)
    aaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGt                                                                                                                                          <  1:3328233/148‑1 (MQ=255)
    aaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGt                                                                                                                                          <  1:3407029/148‑1 (MQ=255)
    aaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGt                                                                                                                                          <  1:3585864/148‑1 (MQ=255)
    aaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGt                                                                                                                                          <  1:3473104/148‑1 (MQ=255)
     aTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTACAAACGTTTCCTTTGTAAGGtt                                                                                                                                         >  2:2927426/1‑148 (MQ=255)
         tttAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGaa                                                                                                                                     >  2:1669633/1‑148 (MQ=255)
            aaaGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACtt                                                                                                                                  <  1:2499719/148‑1 (MQ=255)
              aGGCAACGTTTTCATTTTCAGACATTATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACACGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAg                                                                                                                                <  1:1816741/148‑1 (MQ=255)
              aGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAg                                                                                                                                >  2:2946437/1‑148 (MQ=255)
               ggCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGa                                                                                                                               >  1:1412758/1‑148 (MQ=255)
                gCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAt                                                                                                                              >  2:3749545/1‑148 (MQ=255)
                  aaCGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATACGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAttt                                                                                                                            <  2:24512/148‑1 (MQ=255)
                  aaCGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAttt                                                                                                                            <  2:858555/148‑1 (MQ=255)
                    cGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAttttt                                                                                                                          <  1:3558436/148‑1 (MQ=255)
                             tttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGAGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTc                                                                                                                 <  2:1998911/148‑1 (MQ=255)
                                caGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTg                                                                                                              <  1:2750090/148‑1 (MQ=255)
                                 aGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGa                                                                                                             >  2:3785061/1‑148 (MQ=255)
                                    cATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCt                                                                                                          >  2:1667894/1‑148 (MQ=255)
                                     aTGATAGATTTTCCACACAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGTTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTg                                                                                                         <  2:3771724/148‑1 (MQ=255)
                                      tGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATACGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTgg                                                                                                        <  1:2012949/148‑1 (MQ=255)
                                        aTAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGAt                                                                                                      >  2:3092768/1‑148 (MQ=255)
                                           gATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCat                                                                                                   <  2:1854620/148‑1 (MQ=255)
                                           gATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCat                                                                                                   <  1:1769177/148‑1 (MQ=255)
                                            aTTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCata                                                                                                  <  2:587328/148‑1 (MQ=255)
                                              tttCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCAtata                                                                                                <  1:760697/148‑1 (MQ=255)
                                              tttCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCAtata                                                                                                >  1:3311783/1‑148 (MQ=255)
                                               ttCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATaa                                                                                               >  1:955379/1‑148 (MQ=255)
                                                   ccAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTGTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCa                                                                                           <  2:2879329/148‑1 (MQ=255)
                                                     aGAAAAAAGAGAGTGAGGCGAACATACGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAgt                                                                                         <  2:3046088/148‑1 (MQ=255)
                                                       aaaaaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAgtgt                                                                                       >  1:912156/1‑148 (MQ=255)
                                                        aaaaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTa                                                                                      <  2:1758484/148‑1 (MQ=255)
                                                         aaaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAt                                                                                     >  2:3616320/1‑148 (MQ=255)
                                                         aaaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAt                                                                                     >  1:3689827/1‑148 (MQ=255)
                                                          aaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAtt                                                                                    >  1:1029626/1‑148 (MQ=255)
                                                             gagaTTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttctt                                                                                 >  1:3211920/1‑148 (MQ=255)
                                                              agaTTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttcttc                                                                                <  1:2503516/148‑1 (MQ=255)
                                                               gaTTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTcc                                                                               <  1:203419/148‑1 (MQ=255)
                                                                  tGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCttt                                                                            <  2:1827160/148‑1 (MQ=255)
                                                                  tGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCttt                                                                            <  2:3311736/148‑1 (MQ=255)
                                                                    aGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTAt                                                                          <  1:655276/148‑1 (MQ=255)
                                                                    aGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTAt                                                                          <  1:2946437/148‑1 (MQ=255)
                                                                    aGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTGTAt                                                                          >  2:1971647/1‑148 (MQ=255)
                                                                      gCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTAtct                                                                        <  2:1373156/148‑1 (MQ=255)
                                                                         aaCATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCa                                                                     >  2:1488076/1‑148 (MQ=255)
                                                                           cATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTACGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTAAGTGTATTCTTCCTTTATCTCCAGt                                                                   >  1:1447247/1‑148 (MQ=255)
                                                                              aaGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGc                                                                >  1:3718476/1‑148 (MQ=255)
                                                                              aaGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGc                                                                >  1:1449241/1‑148 (MQ=255)
                                                                              aaGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGc                                                                >  2:435518/1‑148 (MQ=255)
                                                                                ggTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGccc                                                              >  2:261092/1‑148 (MQ=255)
                                                                                ggTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGccc                                                              >  2:2113952/1‑148 (MQ=255)
                                                                                ggTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGccc                                                              >  2:712004/1‑148 (MQ=255)
                                                                                     gCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCa                                                         <  2:3815689/148‑1 (MQ=255)
                                                                                           aaaaTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTAt                                                   >  2:1421066/1‑148 (MQ=255)
                                                                                            aaaTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATg                                                  <  1:3622772/148‑1 (MQ=255)
                                                                                            aaaTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATg                                                  <  2:3631575/148‑1 (MQ=255)
                                                                                             aaTGTCATTAATGGAAAATGTCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATgg                                                 <  2:3498978/148‑1 (MQ=255)
                                                                                              aTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGt                                                >  1:1662343/1‑148 (MQ=255)
                                                                                              aTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGt                                                >  2:2052/1‑148 (MQ=255)
                                                                                               tGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCAATCACTGTATGGTg                                               >  1:2343268/1‑148 (MQ=255)
                                                                                                gTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGc                                              >  1:1153684/1‑148 (MQ=255)
                                                                                                  cATTAATGGAAAATTGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCt                                            <  1:2198148/148‑1 (MQ=255)
                                                                                                  cATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCt                                            >  2:1238125/1‑148 (MQ=255)
                                                                                                  cATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCt                                            <  1:830404/148‑1 (MQ=255)
                                                                                                    ttAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTg                                          <  2:1521359/148‑1 (MQ=255)
                                                                                                    ttAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTg                                          <  2:3420350/148‑1 (MQ=255)
                                                                                                     tAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGt                                         <  2:948622/148‑1 (MQ=255)
                                                                                                      aaTGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGtt                                        <  2:2765303/148‑1 (MQ=255)
                                                                                                       aTGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGttt                                       >  1:2015596/1‑148 (MQ=255)
                                                                                                               tGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCTTATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGt                               >  1:1849336/1‑148 (MQ=255)
                                                                                                               tGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGt                               >  1:3074988/1‑148 (MQ=255)
                                                                                                                  cTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACa                            >  2:1015441/1‑148 (MQ=255)
                                                                                                                    cGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATc                          >  1:480947/1‑148 (MQ=255)
                                                                                                                      aTAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGc                        >  1:3415974/1‑148 (MQ=255)
                                                                                                                        aaCAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAg                      >  2:2142139/1‑148 (MQ=255)
                                                                                                                         aCAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAgg                     >  1:3787666/1‑148 (MQ=255)
                                                                                                                          cAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGa                    >  1:2995323/1‑148 (MQ=255)
                                                                                                                          cAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTAGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGa                    <  2:1433682/148‑1 (MQ=255)
                                                                                                                           aaGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGaa                   >  2:3342232/1‑148 (MQ=255)
                                                                                                                           aaGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGaa                   >  2:1607956/1‑148 (MQ=255)
                                                                                                                           aaGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGaa                   <  1:206433/148‑1 (MQ=255)
                                                                                                                            agagGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAat                  <  2:677409/148‑1 (MQ=255)
                                                                                                                            agagGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAat                  <  2:560308/148‑1 (MQ=255)
                                                                                                                                gTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATaa              <  2:1895741/148‑1 (MQ=255)
                                                                                                                                gTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATaa              >  2:785334/1‑148 (MQ=255)
                                                                                                                                 tAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAAc             >  1:2870517/1‑148 (MQ=255)
                                                                                                                                  aaaaaCGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAAcc            >  2:1004675/1‑148 (MQ=255)
                                                                                                                                  aaaaaCGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAAcc            >  1:2539052/1‑148 (MQ=255)
                                                                                                                                     aaCGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTca         >  1:1681486/1‑148 (MQ=255)
                                                                                                                                       cGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTcaca       <  2:1153684/148‑1 (MQ=255)
                                                                                                                                         tttCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACg     >  1:1873835/1‑148 (MQ=255)
                                                                                                                                          ttCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGa    >  1:2335224/1‑148 (MQ=255)
                                                                                                                                           tCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAg   <  2:3439695/148‑1 (MQ=255)
                                                                                                                                            ccTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGt  <  2:2938861/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGT  >  CP009273/2751414‑2751701

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