BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,318,228	C→A	64.9%	S251* (TCA→TAA)	yciV →	PHP domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,318,228	0	C	A	64.9%	128.6 / 242.8	202	S251* (TCA→TAA)	yciV	PHP domain protein
Reads supporting (aligned to +/- strand): ref base C (37/34); new base A (72/59); total (109/93)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 209 total aligned reads displayed.

CAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTT  >  CP009273/1318082‑1318370

cAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATca                                                                                                                                               <  2:2064820/148‑1 (MQ=255)

  aTCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATcaca                                                                                                                                             <  1:90435/148‑1 (MQ=255)

  aTCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATcaca                                                                                                                                             <  1:636023/148‑1 (MQ=255)

  aTCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATaaca                                                                                                                                             <  2:2488824/148‑1 (MQ=255)

     ttAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCCTTTATGGGCATAACAAgg                                                                                                                                          >  2:260256/1‑148 (MQ=255)

     ttAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAgg                                                                                                                                          <  1:2253184/148‑1 (MQ=255)

      tAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGa                                                                                                                                         <  2:2673140/148‑1 (MQ=255)

         tGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCt                                                                                                                                      >  1:2008593/1‑148 (MQ=255)

         tGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCt                                                                                                                                      >  1:2008610/1‑148 (MQ=255)

         tGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCt                                                                                                                                      >  1:1661286/1‑148 (MQ=255)

          gCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTg                                                                                                                                     >  1:770504/1‑148 (MQ=255)

          gCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTg                                                                                                                                     <  2:364036/148‑1 (MQ=255)

          gCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTg                                                                                                                                     <  1:1621548/148‑1 (MQ=255)

           cTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGa                                                                                                                                    <  2:1835219/148‑1 (MQ=255)

            tAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGAt                                                                                                                                   <  1:99628/148‑1 (MQ=255)

             aaaTGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGAtt                                                                                                                                  >  2:2094777/1‑148 (MQ=255)

             aaaTGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGAtt                                                                                                                                  <  2:3533241/148‑1 (MQ=255)

             aaaTGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGAtt                                                                                                                                  <  1:412888/148‑1 (MQ=255)

             aaaTGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGAtt                                                                                                                                  <  2:552358/148‑1 (MQ=255)

                      aaaaGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGc                                                                                                                         <  1:3270722/148‑1 (MQ=255)

                      aaaaGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGc                                                                                                                         >  1:1201765/1‑148 (MQ=255)

                       aaaGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGcc                                                                                                                        >  2:1901101/1‑148 (MQ=255)

                        aaGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCa                                                                                                                       >  1:1431729/1‑148 (MQ=255)

                          gATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATg                                                                                                                     >  1:514635/1‑148 (MQ=255)

                           aTTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGc                                                                                                                    <  2:2651318/148‑1 (MQ=255)

                            ttGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGcc                                                                                                                   >  1:2855956/1‑148 (MQ=255)

                             tGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGccc                                                                                                                  <  2:1349831/148‑1 (MQ=255)

                                 tGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTgg                                                                                                              >  1:1913961/1‑148 (MQ=255)

                                 tGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTgg                                                                                                              <  1:1146466/148‑1 (MQ=255)

                                 tGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCAGTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTgg                                                                                                              >  1:1327489/1‑148 (MQ=255)

                                  gcgcATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGa                                                                                                             <  2:2745079/148‑1 (MQ=255)

                                      aTTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGa                                                                                                         <  2:2827156/148‑1 (MQ=255)

                                       ttttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGaa                                                                                                        >  1:2720522/1‑148 (MQ=255)

                                       ttttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGaa                                                                                                        >  1:935629/1‑148 (MQ=255)

                                        tttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAAc                                                                                                       <  2:3081252/148‑1 (MQ=255)

                                        tttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAAc                                                                                                       >  1:1731884/1‑148 (MQ=255)

                                         ttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACt                                                                                                      <  1:3317741/148‑1 (MQ=255)

                                           gCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTgg                                                                                                    <  2:1095643/148‑1 (MQ=255)

                                            ccGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTggg                                                                                                   <  1:1073101/148‑1 (MQ=255)

                                            ccGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTggg                                                                                                   <  2:474421/148‑1 (MQ=255)

                                               aaCACCACGGTGACGCGATGGAAGTCGAGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCg                                                                                                <  1:685786/148‑1 (MQ=255)

                                                acacCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGt                                                                                               >  1:1784461/1‑148 (MQ=255)

                                                 caccacGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTa                                                                                              >  2:1173571/1‑148 (MQ=255)

                                                         tGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTgg                                                                                      >  2:1962675/1‑148 (MQ=255)

                                                                   gAAGTCGCGCAGTGCCAGCAATCGCCCAATTAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAg                                                                            >  1:1510889/1‑148 (MQ=255)

                                                                    taGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAgg                                                                           >  1:269577/2‑148 (MQ=255)

                                                                      gTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCg                                                                         >  1:1086407/1‑148 (MQ=255)

                                                                       tCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGt                                                                        <  1:2108634/148‑1 (MQ=255)

                                                                           gcAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGaa                                                                    <  2:1079143/148‑1 (MQ=255)

                                                                                 ccAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTc                                                              >  1:851144/1‑148 (MQ=255)

                                                                                     cAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGc                                                          >  1:983211/1‑148 (MQ=255)

                                                                                           cccAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTAt                                                    <  1:2209769/148‑1 (MQ=255)

                                                                                           cccAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTAt                                                    >  2:4269/1‑148 (MQ=255)

                                                                                             cAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATgg                                                  <  2:1728258/148‑1 (MQ=255)

                                                                                             cAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATgg                                                  >  1:1959946/1‑148 (MQ=255)

                                                                                             cAATGAACGCACCCAGTTGGCGACCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATgg                                                  >  2:1599588/1‑148 (MQ=255)

                                                                                              aaTGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATggg                                                 >  2:3373778/1‑148 (MQ=255)

                                                                                              aaTGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATggg                                                 >  2:2195892/1‑148 (MQ=255)

                                                                                               aTGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGa                                                >  2:814945/1‑148 (MQ=255)

                                                                                               aTGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGa                                                >  2:2630270/1‑148 (MQ=255)

                                                                                                 gAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAAc                                              >  2:1801962/1‑148 (MQ=255)

                                                                                                      tACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCg                                         >  2:2469754/1‑148 (MQ=255)

                                                                                                        cccAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCa                                       >  1:3374304/1‑148 (MQ=255)

                                                                                                            gCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAac                                   >  2:1288803/1‑148 (MQ=255)

                                                                                                             cTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAaca                                  >  2:2546515/1‑148 (MQ=255)

                                                                                                                gcggcCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcacca                               >  1:3339254/1‑148 (MQ=255)

                                                                                                                 cggcCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcaccac                              <  1:3205021/148‑1 (MQ=255)

                                                                                                                   gcCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACag                            >  2:1336937/1‑148 (MQ=255)

                                                                                                                    cccTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACaga                           >  1:3061494/1‑148 (MQ=255)

                                                                                                                    cccTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACaga                           >  2:2904969/1‑148 (MQ=255)

                                                                                                                    cccTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACaga                           >  2:1616182/1‑148 (MQ=255)

                                                                                                                     ccTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagag                          <  2:1104844/148‑1 (MQ=255)

                                                                                                                     ccTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagag                          <  1:1343079/148‑1 (MQ=255)

                                                                                                                      cTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagaga                         <  1:106616/148‑1 (MQ=255)

                                                                                                                      cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCGGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTCTGGGAACAGCCGCAGAACACCACagaga                         >  2:2382196/1‑148 (MQ=255)

                                                                                                                      cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagaga                         >  1:108297/1‑148 (MQ=255)

                                                                                                                         gcgcGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGACGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAggg                      <  1:2779443/148‑1 (MQ=255)

                                                                                                                          cgcgTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGa                     >  2:2935054/1‑148 (MQ=255)

                                                                                                                           gcgTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGaa                    <  1:2040733/148‑1 (MQ=255)

                                                                                                                             gTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAAtt                  <  1:704135/148‑1 (MQ=255)

                                                                                                                              tCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                 >  1:1937039/1‑148 (MQ=255)

                                                                                                                              tCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                 >  1:1078313/1‑148 (MQ=255)

                                                                                                                              tCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                 <  1:2971654/148‑1 (MQ=255)

                                                                                                                              tCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                 <  2:1882585/148‑1 (MQ=255)

                                                                                                                              tCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                 >  1:3324911/1‑148 (MQ=255)

                                                                                                                              tCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                 <  2:514635/148‑1 (MQ=255)

                                                                                                                               cAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTAt                >  1:2057791/1‑148 (MQ=255)

                                                                                                                                 gCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGa              <  1:3491313/148‑1 (MQ=255)

                                                                                                                                  catcatTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAg             >  2:1335726/1‑148 (MQ=255)

                                                                                                                                   atcatTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGc            <  2:1398410/148‑1 (MQ=255)

                                                                                                                                     catTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGCCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCa          >  2:232372/1‑148 (MQ=255)

                                                                                                                                     catTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGCCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCa          <  1:232372/148‑1 (MQ=255)

                                                                                                                                     catTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCa          >  1:1188082/1‑148 (MQ=255)

                                                                                                                                      atTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAg         <  1:1962675/148‑1 (MQ=255)

                                                                                                                                      atTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAg         >  2:1510929/1‑148 (MQ=255)

                                                                                                                                       tttATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGt        >  2:740958/1‑148 (MQ=255)

                                                                                                                                          aTGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGtttt     <  2:3469027/148‑1 (MQ=255)

                                                                                                                                          aTGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGtttt     <  1:3244508/148‑1 (MQ=255)

                                                                                                                                           tGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGttttt    >  1:2567501/1‑148 (MQ=255)

                                                                                                                                            gggCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGtttttt   <  1:1549882/148‑1 (MQ=255)

CAATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTT  >  CP009273/1318082‑1318370

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