BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,783	0	C	A	68.0%	945.6 / inf	895	I69I (ATC→ATA)	ybcX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (165/121); new base A (275/334); total (440/455)
Fisher's exact test for biased strand distribution p-value = 5.68e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 921 total aligned reads displayed.

TATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAA  >  NC_000913/580691‑580875

   ctcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGt                                                                                           <  2:680499/93‑1 (MQ=255)

    tcCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTg                                                                                          >  1:3994923/1‑93 (MQ=255)

     ccccAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGa                                                                                         <  1:41936/93‑1 (MQ=255)

     ccccAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGa                                                                                         <  1:725818/93‑1 (MQ=255)

      cccaaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAg                                                                                        <  1:4113895/90‑1 (MQ=255)

      ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAg                                                                                        >  2:2661487/1‑93 (MQ=255)

      ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAg                                                                                        >  1:3560345/1‑93 (MQ=255)

         aaaGAGTCGGCCTACCAGGAAGGGCGCTGGGAACCGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtc                                                                                     <  1:5784044/93‑2 (MQ=255)

         aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTg                                                                                     >  1:2514694/1‑93 (MQ=255)

          aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtca                                                                                    >  1:2603421/1‑91 (MQ=255)

          aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtca                                                                                    <  2:371488/93‑3 (MQ=255)

          aaGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGtca                                                                                    <  2:3863122/93‑3 (MQ=255)

             aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtg                                                                                 <  1:4977268/93‑1 (MQ=255)

             aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtg                                                                                 >  1:1474366/1‑93 (MQ=255)

                 ggcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgt                                                                             <  2:5967268/89‑1 (MQ=255)

                  gcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgtt                                                                            >  2:1484801/4‑92 (MQ=255)

                  gcctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACgttgtt                                                                            >  2:5144013/4‑92 (MQ=255)

                   cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                                                           >  2:4669552/1‑93 (MQ=255)

                   cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                                                           >  1:1389489/1‑93 (MQ=255)

                    ctACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTaa                                                                          >  1:29382/2‑93 (MQ=255)

                     tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                         >  2:3973782/1‑93 (MQ=255)

                     tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                         >  2:5042905/1‑93 (MQ=255)

                     tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                         >  2:3080921/1‑93 (MQ=255)

                     tACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAg                                                                         >  2:5229064/1‑93 (MQ=255)

                        cAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCt                                                                      >  1:564259/1‑93 (MQ=255)

                        cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                      >  2:3784289/1‑93 (MQ=255)

                        cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                      >  1:911611/1‑93 (MQ=255)

                         aGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTg                                                                     >  2:4080592/1‑93 (MQ=255)

                         aGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTg                                                                     >  2:421975/1‑93 (MQ=255)

                          ggAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGc                                                                    >  2:2992849/1‑93 (MQ=255)

                          ggAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                                                    >  1:1204904/1‑93 (MQ=255)

                           gAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGcc                                                                   >  2:78288/1‑93 (MQ=255)

                             aGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCg                                                                 <  1:191396/93‑1 (MQ=255)

                               ggCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtg                                                               <  1:2049669/93‑1 (MQ=255)

                                gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                                              <  2:1620366/93‑1 (MQ=255)

                                gcgcTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgt                                                              <  1:5077444/93‑1 (MQ=255)

                                gcgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                                              <  2:915693/93‑1 (MQ=255)

                                 cgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTc                                                             <  1:4258949/93‑1 (MQ=255)

                                   cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                                                           <  2:4807059/93‑1 (MQ=255)

                                   cTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTgg                                                           <  2:1273341/93‑1 (MQ=255)

                                    tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                          <  1:1081075/93‑2 (MQ=255)

                                    tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                          >  1:2200337/1‑92 (MQ=255)

                                    tGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTggc                                                          >  1:4423866/1‑92 (MQ=255)

                                     gggAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTcagct                                                         >  1:5187925/1‑89 (MQ=255)

                                      ggtaaCACTGCCCTTTCAGCGGGCCAGCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAACGTGGTGAAGTCTGCCCGTGTCGGCTa                                                        >  1:5500491/4‑93 (MQ=255)

                                      ggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                                                        <  2:6044473/93‑1 (MQ=255)

                                        aaaCACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAtt                                                      >  2:1366912/1‑93 (MQ=255)

                                              tGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                                <  1:1820950/93‑1 (MQ=255)

                                              tGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAaaaa                                                <  1:991844/93‑1 (MQ=255)

                                               gCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAt                                               >  2:2561961/1‑93 (MQ=255)

                                                ccgtttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtg                                              <  1:5382829/90‑1 (MQ=255)

                                                cccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTATGCCCGTGTCGGTTATTCCAAAAtg                                              >  1:2753450/1‑93 (MQ=255)

                                                 ccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGGCCGTGTCGGTTATTCCAAAAtgc                                             >  2:2430490/1‑93 (MQ=255)

                                                 ccTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgc                                             <  2:1138984/93‑1 (MQ=255)

                                                   tttCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                                           >  1:3053289/1‑93 (MQ=255)

                                                   tttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttg                                           <  1:5967634/93‑1 (MQ=255)

                                                    ttCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAAtgttgc                                          <  1:3855858/93‑1 (MQ=255)

                                                     tCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATgttgct                                         <  1:1710890/93‑1 (MQ=255)

                                                        gCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg                                      <  2:531463/93‑1 (MQ=255)

                                                           ggCCACCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt                                   >  1:4689676/1‑93 (MQ=255)

                                                           cgcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGgtgt                                   <  1:4858610/89‑1 (MQ=255)

                                                            gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGtt                                  >  1:3948452/4‑93 (MQ=255)

                                                            gcgatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGtt                                  <  1:825418/90‑1 (MQ=255)

                                                              gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                <  1:5876504/92‑1 (MQ=255)

                                                              gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                <  2:748022/92‑1 (MQ=255)

                                                              gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                <  2:1861053/92‑1 (MQ=255)

                                                              gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                <  1:2697952/92‑1 (MQ=255)

                                                              gatcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTa                                <  2:5384994/92‑1 (MQ=255)

                                                               atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                               >  2:4422660/1‑93 (MQ=255)

                                                               atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                               >  2:435177/1‑93 (MQ=255)

                                                               atcatGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                               >  2:5499905/1‑93 (MQ=255)

                                                               atcatGAATGCGATGGGTAACGACTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTAt                               <  1:4488412/93‑1 (MQ=255)

                                                                tcatGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATg                              >  1:2828952/1‑93 (MQ=255)

                                                                 catGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGc                             >  1:324836/1‑93 (MQ=255)

                                                                 catGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGc                             >  2:2208126/1‑93 (MQ=255)

                                                                   tGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGcgt                           >  2:3515160/1‑91 (MQ=255)

                                                                     aaTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                         >  1:3170759/1‑93 (MQ=255)

                                                                     aaTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                         <  2:1789574/93‑1 (MQ=255)

                                                                      aTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                        <  2:5577478/93‑1 (MQ=255)

                                                                      aTGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                        >  1:1584997/1‑93 (MQ=255)

                                                                        gCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTAtttt                      <  2:1931423/93‑1 (MQ=255)

                                                                         cGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTta                     >  2:5868149/1‑93 (MQ=255)

                                                                          gATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtat                    >  2:4892401/1‑93 (MQ=255)

                                                                           aTGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtatt                   >  2:3222396/1‑92 (MQ=255)

                                                                            tGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattc                  <  1:3737322/93‑3 (MQ=255)

                                                                            tGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTtattc                  <  1:2411029/93‑3 (MQ=255)

                                                                            tGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATag                  <  1:2308628/93‑1 (MQ=255)

                                                                             gggCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATaga                 >  1:523772/1‑93 (MQ=255)

                                                                                   cGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACaa           <  2:1690561/93‑1 (MQ=255)

                                                                                    gattACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAag          <  1:334278/90‑1 (MQ=255)

                                                                                     attACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc         <  2:2732109/91‑1 (MQ=255)

                                                                                     attACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc         <  1:663946/91‑1 (MQ=255)

                                                                                     attACATACGTGAGGTCAACGGTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagc         <  1:606042/91‑1 (MQ=255)

                                                                                      ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca        >  2:283773/2‑93 (MQ=255)

                                                                                      ttACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagca        >  2:5799989/2‑93 (MQ=255)

                                                                                      cTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca        >  2:3586016/1‑93 (MQ=255)

                                                                                       tACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagcag       >  2:3394201/1‑93 (MQ=255)

                                                                                       tACATACGTGAGGTCAACGTTGTTAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTATGCGTATTTTATTCAGCACAagcag       >  2:5670425/1‑93 (MQ=255)

                                                                                          aTCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGATGTTTATGCCTACTTTATAGAGCATAAGCAgcgc    <  1:457389/93‑1 (MQ=255)

                                                                                           tCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCa   <  2:2067145/93‑1 (MQ=255)

TATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAA  >  NC_000913/580691‑580875

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