BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,941	A→T	44.2%	intergenic (‑1024/‑144)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,941	0	A	T	44.2%	28.2 / 158.5	127	intergenic (‑1024/‑144)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (40/29); new base T (29/25); total (72/55)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-02

Only 100 of 128 total aligned reads displayed.

TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTT  >  CP000730/1929802‑1930077

tgatATAATAGTAGAGTCGCCTATCTCTCAGGAGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa                                                                                                                                         >  1:653723/1‑141 (MQ=255)

 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                       >  1:104876/3‑141 (MQ=255)

 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                       >  2:404668/3‑141 (MQ=255)

 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTACTTGCTTAGt                                                                                                                                       >  2:491458/3‑141 (MQ=255)

 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATGATACTTATTGGCTTAGt                                                                                                                                       <  1:219423/139‑1 (MQ=255)

 ggtacaataatAGAGTCGCCTATCTCTCAGCCGTCAATT‑‑GACGAAGAGAGGAGATGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAga                                                                                                                                       >  2:554136/6‑140 (MQ=255)

   tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                     >  2:133234/1‑140 (MQ=255)

   tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTaa                                                                                                                                     >  2:698714/1‑140 (MQ=255)

   tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                     >  1:255172/1‑140 (MQ=255)

    ataataataGAGTCGCCTATCTCTCAGGCGTCTATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATATCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                    <  1:609189/141‑3 (MQ=255)

    ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                    <  1:517830/141‑3 (MQ=255)

    ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                    >  1:777417/1‑139 (MQ=255)

    ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGGTAAGtca                                                                                                                                     >  2:578842/1‑139 (MQ=255)

    ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAg                                                                                                                                     >  2:202846/1‑141 (MQ=255)

     taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                                                    <  1:230257/141‑1 (MQ=255)

     taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                                                    >  2:65728/1‑141 (MQ=255)

     taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                                                    >  1:289535/1‑141 (MQ=255)

       ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                 >  1:291904/5‑141 (MQ=255)

       ataataGGGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TCGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                 <  2:601457/141‑1 (MQ=255)

        taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAAGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                <  2:852428/138‑1 (MQ=255)

         aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                               <  1:870882/139‑1 (MQ=255)

         aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑GTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTTATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                               <  2:671873/139‑1 (MQ=255)

           aagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                             <  2:342829/140‑3 (MQ=255)

            tgagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                                             <  1:846571/140‑1 (MQ=255)

            agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                                             <  1:846460/141‑1 (MQ=255)

            agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                            >  2:630967/1‑138 (MQ=255)

            agagTCGCCTACCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                            <  2:242657/141‑4 (MQ=255)

             gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                                            >  1:316406/1‑141 (MQ=255)

              agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                                                                                                           >  1:521234/1‑141 (MQ=255)

              agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                          <  1:677652/141‑6 (MQ=255)

               gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGATTAGTAAACGGaacaaata                                                                                                                         >  2:590923/1‑135 (MQ=255)

               gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGTTTGCGTTTTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                         >  1:189030/1‑135 (MQ=255)

               gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGAATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACACaacaaata                                                                                                                         >  2:597656/1‑135 (MQ=255)

                tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                                        <  2:73604/141‑8 (MQ=255)

                ccGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                                        <  1:737207/140‑8 (MQ=255)

                 cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                                                                                        <  1:239109/141‑1 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                                                                                       >  2:432488/1‑141 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGAA‑TAA‑G‑G‑T‑GATGGTTATTGTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                                                                                       <  2:277341/141‑1 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                      >  1:889378/1‑132 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                      <  2:623208/141‑10 (MQ=255)

                  gCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                      <  1:878105/141‑10 (MQ=255)

                   ccTATCTCTCAGGCATCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGccacaaataatta                                                                                                                     >  1:595418/1‑129 (MQ=255)

                     tATCTCTCAGGCTTCAATT‑‑GGCGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGTTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                   <  2:365251/141‑13 (MQ=255)

                     tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                   >  2:530364/1‑129 (MQ=255)

                     tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                   <  2:59483/141‑13 (MQ=255)

                     tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGGTTGCGTATTATACTTATTGGCTGAGTAGACACaacaaataattaaa                                                                                                                   >  1:141676/1‑129 (MQ=255)

                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                  <  1:819367/141‑14 (MQ=255)

                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                  >  1:364646/1‑128 (MQ=255)

                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                  >  1:606306/1‑128 (MQ=255)

                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatcattaatg                                                                                                                  >  1:59342/1‑128 (MQ=255)

                       tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCTTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTg                                                                                                                  >  1:802774/1‑141 (MQ=255)

                       tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattattag                                                                                                                 <  2:790341/141‑15 (MQ=255)

                       tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                                 <  1:728110/141‑15 (MQ=255)

                       tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                                 <  2:633684/141‑15 (MQ=255)

                       tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                                 >  2:808231/1‑127 (MQ=255)

                       tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATAATTATTGGCTGAGTAAACCCaacaactaattaatag                                                                                                                 >  1:790687/1‑127 (MQ=255)

                       tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGAACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                                 <  2:781193/141‑15 (MQ=255)

                        ctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACACAATAGCTGa                                                                                                                 >  1:149792/1‑141 (MQ=255)

                        ctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGATCACGCGGGTTCACGTAATAGTCCCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTGCTGATTTGCTAAGTAGACGCAATACAAAATAGGTGa                                                                                                                 >  1:901647/1‑141 (MQ=255)

                             aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                                                                                            >  1:787245/1‑141 (MQ=255)

                                   cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                                                                                                      >  1:111358/1‑141 (MQ=255)

                                       ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                                                                                  >  1:610001/1‑141 (MQ=255)

                                         tgACTCAGAGAGGAGGTGTA‑TAA‑G‑A‑T‑GATGCTGATCTTCGTGCACACGATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                                                <  1:802507/140‑1 (MQ=255)

                                         aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                                                <  1:832211/141‑1 (MQ=255)

                                              cAGAGAGGAGGTGTC‑TAG‑G‑G‑T‑GAGGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGAAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                                           <  2:100939/141‑1 (MQ=255)

                                              cAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                                           <  2:714853/141‑1 (MQ=255)

                                              cAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                                           <  2:714883/141‑1 (MQ=255)

                                              cAGAGAGGAGGGGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                                           <  1:368091/141‑1 (MQ=255)

                                                 agagGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcc                                                                                        <  2:653723/141‑1 (MQ=255)

                                                 agagGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTAACATATACGCCCACCAATAAAAAAcc                                                                                        >  2:82701/1‑141 (MQ=255)

                                                   aggaggTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCAACAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAAAAAATAGGTGACATATAGCCGCACGAATAAAAATaccc                                                                                      >  1:581602/1‑141 (MQ=255)

                                                   aggaggAGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcccc                                                                                      >  2:597906/1‑141 (MQ=255)

                                                         tgtA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAATTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa                                                                                >  2:736959/1‑141 (MQ=255)

                                                                a‑G‑G‑T‑GATGCTTATTTACGTTCACATCTTAGCACCTGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGTCGCAATACAAAATGGGTGACATATAGCCGCACCAATATAAATCCCCTCACTACCGCaa                                                                          <  1:551202/141‑1 (MQ=255)

                                                                a‑G‑G‑T‑CATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCCACTGCGTTTTTTTCTTATTGGCTACGTAGACGCAATACAACATCGGTCACATATAGCCGCACCAATACAAACCCCCTCACTACCGCaa                                                                          <  2:477213/141‑1 (MQ=255)

                                                                a‑G‑G‑G‑GATGCTTATTTTCGTTCATATCATAGCACCAGTCAGCAGTGGCTGTGCCATTGCGTTGTTTTCTTATTGGCTGAGTAGACGCAATACAAAATAGGCGACATATAGCCGCACCAATAACAATCACGTCGTTACCGCaa                                                                          <  2:106536/141‑1 (MQ=255)

                                                                         aTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAATAgg                                                                     >  1:229113/1‑140 (MQ=255)

                                                                           gCTTATTTTCGTTCACATCATAGCACCAGTCATAAGTGGCTGGGGCATTGCGTGTTTTTCTTATTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCACATCACTACCGCAAATAGTGa                                                                   >  1:352959/1‑141 (MQ=255)

                                                                            cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                  >  2:161154/1‑141 (MQ=255)

                                                                              tATTTTCGTTTACATCATAGCACCAGTCATCAGTGGCTTTGCCATTGCGTTTTTTCCTTATTGTCTAAGTAGACGCAATACAAAATAAGTCACATATAGCCGCACCAATAAAAAACCCCTCACTACCGCAAATAGTGAggg                                                                <  1:777971/141‑1 (MQ=255)

                                                                              tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATACAAATCCCCTCACTCCCCCCAACAGCGAggg                                                                >  2:106029/1‑141 (MQ=255)

                                                                              tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCGATACAAAATAGGTGACATATAGCCGCACCAATACAAATCCCCTCACTACAGCAAACAGCGCggg                                                                >  2:658861/1‑141 (MQ=255)

                                                                               tttttCGTTCCCATCATAGCACCGGTCGTCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAACTCCCCTCACTACCGCAAATAGTGAgggg                                                               <  2:229402/140‑1 (MQ=255)

                                                                                ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                              <  1:459763/141‑1 (MQ=255)

                                                                                ttttCGTTCACATCAAAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                              <  2:140705/141‑1 (MQ=255)

                                                                                ttttCGGTCACATCATAGGACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                              <  2:208213/141‑1 (MQ=255)

                                                                                  ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                                            >  1:183539/1‑141 (MQ=255)

                                                                                      ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGTTTGgtg                                                        >  1:248434/1‑141 (MQ=255)

                                                                                        cacaTCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGta                                                      >  1:401170/1‑141 (MQ=255)

                                                                                              atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATGAGTa                                                >  2:115433/1‑141 (MQ=255)

                                                                                                       gTCATCAGTGGCTGTGGCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt                                       <  1:688181/141‑1 (MQ=255)

                                                                                                       gTCAGCAGTGGCGGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTACACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCCCCATACCGGAACTAGTGAGGGGATTGGTGTATAAGTAAATCCTTAt                                       >  1:642473/1‑141 (MQ=255)

                                                                                                        tcatcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                                      <  2:610001/141‑1 (MQ=255)

                                                                                                         catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTTAATACTTGttt                                     >  1:242302/1‑141 (MQ=255)

                                                                                                            caGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATCGGATACAGACAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCg                                  <  1:432488/141‑1 (MQ=255)

                                                                                                                    gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAGACAAAATAGGGGACATATAGCCGCACCAATAAACATCCCCTCACTACCGAAAGTAGAGAGGGGGTTGGTGTATAAGTAAATACTTATTTTACTTGGCTTa                          >  1:423806/1‑141 (MQ=255)

                                                                                                                                   tttCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt           <  2:316406/141‑1 (MQ=255)

                                                                                                                                   tttCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt           <  2:264965/141‑1 (MQ=255)

                                                                                                                                   tttCTGATTGGCTAAGTAGACGCAATACAAAAGAGGTGACATATAGCCGCACCAATAAACATCCCCTCCCTACCGCAAATAGTGAGGGGATTGGTGTATAATGAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt           <  1:180522/141‑1 (MQ=255)

                                                                                                                                            ggCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTGTCGTTGTCTTAATTATACTGCTAATTTTTCttttt  >  1:168411/1‑141 (MQ=255)

TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTT  >  CP000730/1929802‑1930077

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