Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,751,558 G→T 93.6% 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,5580GT93.6% 704.8 / 35.9 248T134T (ACC→ACAyfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (7/9);  new base T (113/119);  total (120/128)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 249 total aligned reads displayed.

CGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGT  >  CP009273/2751412‑2751701
                                                                                                                                                  |                                                                                                                                               
 gATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGt                                                                                                                                               >  1:2850478/1‑148 (MQ=255)
   tGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTaa                                                                                                                                             >  1:147972/1‑148 (MQ=255)
     aaaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAgg                                                                                                                                           <  2:1914040/148‑1 (MQ=255)
     aaaTGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAgg                                                                                                                                           >  1:747193/1‑148 (MQ=255)
         gCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTg                                                                                                                                       >  2:1515523/1‑148 (MQ=255)
          cTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGa                                                                                                                                      <  2:2882856/148‑1 (MQ=255)
           tttAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGaa                                                                                                                                     >  1:2279093/1‑148 (MQ=255)
               aaGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTa                                                                                                                                 >  1:610532/1‑148 (MQ=255)
                aGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATACGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAg                                                                                                                                <  1:457275/148‑1 (MQ=255)
                aGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAg                                                                                                                                >  2:2211479/1‑148 (MQ=255)
                 ggCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGACCATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGa                                                                                                                               <  1:2582936/148‑1 (MQ=255)
                  gCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAt                                                                                                                              >  2:526277/1‑148 (MQ=255)
                   cAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAtt                                                                                                                             >  1:1403722/1‑148 (MQ=255)
                   cAACGATTTCATTTTCAGACATGATAGATTTTCCCAATAAAAAAGAGATTGAGGAGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAtt                                                                                                                             <  2:3031117/148‑1 (MQ=255)
                     aCGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGAtttt                                                                                                                           >  1:1002554/1‑148 (MQ=255)
                     aCGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGAtttt                                                                                                                           <  2:1394814/148‑1 (MQ=255)
                       gTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTc                                                                                                                         <  1:3037282/148‑1 (MQ=255)
                        ttttcattttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTcc                                                                                                                        <  2:16154/148‑1 (MQ=255)
                           tcattttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGc                                                                                                                     <  1:638503/148‑1 (MQ=255)
                            cattttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCa                                                                                                                    <  2:1140854/148‑1 (MQ=255)
                              ttttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCAtt                                                                                                                  <  1:350904/148‑1 (MQ=255)
                               tttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTc                                                                                                                 >  1:2725574/1‑148 (MQ=255)
                                ttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCt                                                                                                                <  2:1724457/148‑1 (MQ=255)
                                ttcaGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCt                                                                                                                <  2:2606562/148‑1 (MQ=255)
                                    gACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAg                                                                                                            <  2:2100133/148‑1 (MQ=255)
                                     aCATGATAGATTTTCCCCAGAAAAAATATATTGAGTCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGATCGATAACAAGAGGTAAAAAAGTTTCCTTTGTAAGGTTTGAACTTATATTTTTCCAGCATTCTTGAGc                                                                                                           >  2:2963321/1‑148 (MQ=255)
                                     aCATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGc                                                                                                           >  1:853859/1‑148 (MQ=255)
                                     aCATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGc                                                                                                           <  1:2673533/148‑1 (MQ=255)
                                          aTAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGAt                                                                                                      >  2:1464691/1‑148 (MQ=255)
                                          aTAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGAt                                                                                                      >  2:2163192/1‑148 (MQ=255)
                                           tAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATc                                                                                                     <  1:1874754/148‑1 (MQ=255)
                                           tAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATc                                                                                                     <  1:2682388/148‑1 (MQ=255)
                                           tAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATc                                                                                                     >  1:1999714/1‑148 (MQ=255)
                                            aGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCa                                                                                                    <  2:1341077/148‑1 (MQ=255)
                                             gATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCat                                                                                                   >  2:2187432/1‑148 (MQ=255)
                                               ttttCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCatat                                                                                                 >  2:2415332/1‑148 (MQ=255)
                                                tttCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCAtata                                                                                                >  1:2634776/1‑148 (MQ=255)
                                                 ttCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATaa                                                                                               <  2:1195027/148‑1 (MQ=255)
                                                    cccAGAAAAAAGAGATTGAGGCGAACATAAGGTGTTGATTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTc                                                                                            <  2:362745/148‑1 (MQ=255)
                                                    cccAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTc                                                                                            >  1:1258444/1‑148 (MQ=255)
                                                     ccAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCa                                                                                           <  1:773709/148‑1 (MQ=255)
                                                      cAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAg                                                                                          <  1:839739/148‑1 (MQ=255)
                                                      cAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAg                                                                                          <  2:1900429/148‑1 (MQ=255)
                                                       aGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAgt                                                                                         >  2:2523859/1‑148 (MQ=255)
                                                         aaaaaaGAGATTGAGGCTAACATAAGGTGTGCCTTGAAAATGTCATTAATGTAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAgtgt                                                                                       <  2:1962661/148‑1 (MQ=255)
                                                         aaaaaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAgtgt                                                                                       <  1:1624180/148‑1 (MQ=255)
                                                           aaaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAt                                                                                     >  2:1067920/1‑148 (MQ=255)
                                                           aaaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAt                                                                                     <  2:2261114/148‑1 (MQ=255)
                                                            aaaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAtt                                                                                    >  1:2995219/1‑148 (MQ=255)
                                                             aaGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttc                                                                                   >  2:1809191/1‑148 (MQ=255)
                                                              agagATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttct                                                                                  <  1:526277/148‑1 (MQ=255)
                                                              agagATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttct                                                                                  >  2:826969/1‑148 (MQ=255)
                                                               gagaTTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttctt                                                                                 >  1:2572600/1‑148 (MQ=255)
                                                               gagaTTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttctt                                                                                 >  1:566193/1‑148 (MQ=255)
                                                                agaTTGAGGCTAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTAttcttc                                                                                <  2:1678139/148‑1 (MQ=255)
                                                                  aTTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCt                                                                              >  1:67018/1‑148 (MQ=255)
                                                                      aGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTAt                                                                          >  1:2523250/1‑148 (MQ=255)
                                                                       ggCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTAtc                                                                         >  1:2644304/1‑148 (MQ=255)
                                                                       ggCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTAtc                                                                         >  2:1145585/1‑148 (MQ=255)
                                                                           aaCATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCa                                                                     <  2:2834929/148‑1 (MQ=255)
                                                                            aCATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAg                                                                    >  1:2196976/1‑148 (MQ=255)
                                                                               tAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTg                                                                 >  1:1135600/1‑148 (MQ=255)
                                                                                 aGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGcc                                                               <  1:2681096/148‑1 (MQ=255)
                                                                                   gtgtGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGATGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCa                                                             <  2:1031276/148‑1 (MQ=255)
                                                                                   gtgtGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCa                                                             <  2:1031156/148‑1 (MQ=255)
                                                                                           tGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGt                                                     >  1:2945492/1‑148 (MQ=255)
                                                                                             aaaaTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTAt                                                   >  2:1046393/1‑148 (MQ=255)
                                                                                              aaaTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATg                                                  <  2:3158903/148‑1 (MQ=255)
                                                                                              aaaTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATg                                                  >  1:453347/1‑148 (MQ=255)
                                                                                               aaTGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATgg                                                 >  1:2154711/1‑148 (MQ=255)
                                                                                                   tCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGcc                                             <  2:1002380/148‑1 (MQ=255)
                                                                                                      ttAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTg                                          >  2:2295906/1‑148 (MQ=255)
                                                                                                       tAATTGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGt                                         <  1:915255/148‑1 (MQ=255)
                                                                                                       tAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGt                                         >  2:872946/1‑148 (MQ=255)
                                                                                                            gAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTccc                                    >  1:1636829/1‑148 (MQ=255)
                                                                                                            gAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTccc                                    <  1:944112/148‑1 (MQ=255)
                                                                                                             aaaaTGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCa                                   >  2:2857371/1‑148 (MQ=255)
                                                                                                               aaTGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATg                                 >  2:664719/1‑148 (MQ=255)
                                                                                                                aTGGCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATgg                                >  1:2950004/1‑148 (MQ=255)
                                                                                                                aTGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATgg                                >  1:3113175/1‑148 (MQ=255)
                                                                                                                   gCTCGATAACAATAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTAc                             <  2:2694411/148‑1 (MQ=255)
                                                                                                                   gCTCGATAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTAc                             <  2:1251308/148‑1 (MQ=255)
                                                                                                                   gCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTAc                             >  1:2324895/1‑148 (MQ=255)
                                                                                                                     tCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACAt                           <  2:1326876/148‑1 (MQ=255)
                                                                                                                        aTAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGc                        >  1:1037096/1‑148 (MQ=255)
                                                                                                                         tAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCa                       >  2:3250335/1‑148 (MQ=255)
                                                                                                                         tAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCa                       <  2:561346/148‑1 (MQ=255)
                                                                                                                         tAACAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCa                       <  2:1999714/148‑1 (MQ=255)
                                                                                                                          aaCAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAg                      >  2:776378/1‑148 (MQ=255)
                                                                                                                           aCAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAgg                     >  1:2438670/1‑148 (MQ=255)
                                                                                                                            cAAGAGGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGa                    <  1:1067749/148‑1 (MQ=255)
                                                                                                                              agagGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAat                  <  2:2725574/148‑1 (MQ=255)
                                                                                                                              agagGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAat                  >  2:814/1‑148 (MQ=255)
                                                                                                                              agagGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCCTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAat                  >  2:1639346/1‑148 (MQ=255)
                                                                                                                               gagGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAata                 >  1:68979/1‑148 (MQ=255)
                                                                                                                               gagGTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAata                 >  1:1522672/1‑148 (MQ=255)
                                                                                                                                 ggTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAAtata               <  2:424611/148‑1 (MQ=255)
                                                                                                                                  gTAAAAACGTTTCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATaa              >  1:1304763/1‑148 (MQ=255)
                                                                                                                                         cGTTTCCTTTGTAATGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTcaca       <  2:2634776/148‑1 (MQ=255)
                                                                                                                                           tttCCTTTGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACg     <  1:2163192/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGATGAAATGCTTTAAAGGCAACGTTTTCATTTTCAGACATGATAGATTTTCCCCAGAAAAAAGAGATTGAGGCGAACATAAGGTGTGCCTTGAAAATGTCATTAATGGAAAATGGCTCGATAACAAGAGGTAAAAACGTTTCCTTGGTAAGGTTTGAACTTAGATTTTTCCAGCATTCTTGAGCTGGATCATATAATTCAGTGTATTCTTCCTTTATCTCCAGTTGCCCATTCACTGTATGGTGCCTTGTTTCCCATGGTACATCGCAGGAATATAACCTCACACGAGT  >  CP009273/2751412‑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: