Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,154 T→A 17.8% intergenic (‑124/+22) dsbB ← / ← nhaB oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,233,1540TA17.8% 363.0 / 74.1 160intergenic (‑124/+22)dsbB/nhaBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (56/74);  new base A (12/16);  total (68/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-02

Only 100 of 164 total aligned reads displayed.

AGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCC  >  NC_000913/1233056‑1233253
                                                                                                  |                                                                                                   
agTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAttt                                                                                                    <  1:907615/100‑1 (MQ=255)
 gTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTta                                                                                                   <  2:3290917/100‑1 (MQ=255)
 gTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTta                                                                                                   <  2:555491/100‑1 (MQ=255)
  tGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtat                                                                                                  >  2:3221241/1‑100 (MQ=255)
       aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                             <  1:548197/100‑1 (MQ=255)
          ttAAACCAATTCATTCGCTGCGCCACCAGTAGGCAGCAAATATAGGACGCAATTTGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCgg                                                                                          <  1:766609/100‑1 (MQ=255)
          ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCgg                                                                                          <  1:2758740/100‑1 (MQ=255)
          ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                          <  2:2452827/100‑1 (MQ=255)
          ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                          <  2:1197784/100‑1 (MQ=255)
          ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                          <  2:1144266/100‑1 (MQ=255)
          ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                          <  2:1136818/100‑1 (MQ=255)
            aaaCCAATTCATTCGCTGCGCCACTAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAACGCAAACCGGGCATTTATATGCCCGGTa                                                                                        <  1:747975/100‑1 (MQ=255)
            aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtt                                                                                        <  2:2104666/100‑2 (MQ=255)
            aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTAGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                             >  1:581286/1‑95 (MQ=255)
            aaaCCAATTCACCCGCTGCTCCACAAGAAAAAAAAAAAAATAGAACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGt                                                                                         <  1:2848384/99‑1 (MQ=255)
             agccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACGGAGCATTGATATGCCCGGTa                                                                                        <  1:3640900/97‑1 (MQ=255)
             aaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTaa                                                                                       >  2:2302360/1‑100 (MQ=255)
             aaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTaa                                                                                       >  1:3759790/1‑100 (MQ=255)
             aaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                        >  1:1457462/1‑99 (MQ=255)
              aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAg                                                                                      >  2:52486/1‑100 (MQ=255)
              aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttg                                                                                      <  2:4401509/100‑4 (MQ=255)
              aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttg                                                                                      <  2:1537925/100‑4 (MQ=255)
               ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgc                                                                                     <  2:688979/100‑5 (MQ=39)
                cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcc                                                                                    <  1:1238459/100‑6 (MQ=39)
                 aaTTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAAATAAAAAAAAAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttg                                                                                   <  1:1682467/100‑1 (MQ=255)
                  attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGTACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgt                                                                                  <  1:793044/100‑1 (MQ=255)
                  attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgt                                                                                  <  1:377748/100‑1 (MQ=255)
                      attcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATc                                                                              >  1:3222543/1‑100 (MQ=255)
                         cGCTGCGCCACCATGAGGCAGCAAAAATAGGACGCTATTGGGTGAGTATCGGGCGAAATGGCAAACCGTTCATTTATATGCTCGGTAAGTTGTTATCaaa                                                                           <  1:4383163/100‑1 (MQ=255)
                         cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                           <  2:2910936/100‑1 (MQ=255)
                         cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                           >  2:1970338/1‑100 (MQ=255)
                         cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                           >  2:3708991/1‑100 (MQ=255)
                          gCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAg                                                                          <  1:2211169/100‑1 (MQ=255)
                            tGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCg                                                                        >  1:2275485/1‑100 (MQ=255)
                             gcgcCACCAGGAGGCAGCAAAAATAGGACGCGATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGt                                                                       <  1:2635219/100‑1 (MQ=255)
                               gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                     >  1:1876733/1‑100 (MQ=255)
                                ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                    <  1:311597/100‑1 (MQ=255)
                                  accaGGAGGCAGCAAAAATAGGACGCAATTGTACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGATAAGTTGTTATCAAAGCGTTg                                                                     <  1:1571421/97‑1 (MQ=255)
                                  accaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTACCAAAGCGTTGCTa                                                                  <  2:2986195/100‑1 (MQ=255)
                                     aggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                               <  1:3337276/100‑1 (MQ=255)
                                         ggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGcc                                                           <  2:2275495/100‑1 (MQ=255)
                                         ggCAGCAAAAATAGGACGCAATTGGACAATTAGCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGTGTTGCTATCCAGcc                                                           <  1:3684493/100‑1 (MQ=255)
                                            agcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                        >  2:1366479/1‑100 (MQ=255)
                                            agcaAAAAAACGACGCAATTGGACCATTACCGGGCGAAAAGTCAAACCGGGCATCTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGccc                                                          <  1:3668105/98‑1 (MQ=255)
                                             gcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                        >  1:2069400/1‑99 (MQ=255)
                                              caAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAttt                                                      >  2:4186042/1‑100 (MQ=255)
                                              caAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAttt                                                      >  1:3139881/1‑100 (MQ=255)
                                              caAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAttt                                                      >  2:3289832/1‑100 (MQ=255)
                                               aaaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                     <  2:4079649/100‑1 (MQ=255)
                                               aaaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                     <  1:968577/100‑1 (MQ=255)
                                                aaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGc                                                    <  1:1463726/100‑1 (MQ=255)
                                                aaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGc                                                    <  2:876765/100‑1 (MQ=255)
                                                 aaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                   <  1:985198/100‑1 (MQ=255)
                                                 aaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                   <  1:3556071/100‑1 (MQ=255)
                                                 aaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                   <  2:3139898/100‑1 (MQ=255)
                                                     aGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                               <  2:872495/100‑1 (MQ=255)
                                                       gACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAcc                                             >  2:2557444/1‑100 (MQ=255)
                                                           cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         >  2:3607421/1‑100 (MQ=255)
                                                           cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         >  2:169684/1‑100 (MQ=255)
                                                            aaTTGGACAATTATCGCGCGAAAAGGCAAACCGGGCATTTCTATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGTCCATTTGCATAAACCATACg                                        <  1:1977229/100‑1 (MQ=255)
                                                             aTTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCgg                                       <  2:924500/100‑1 (MQ=255)
                                                              ttGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGt                                      <  1:3279107/100‑1 (MQ=255)
                                                              ttGGACAATTATCGGGCGAAAAGGCAAACCGGGCATCTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGt                                      <  2:3088671/100‑1 (MQ=255)
                                                                ggACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTa                                    >  1:343028/1‑100 (MQ=255)
                                                                ggACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGGTa                                    >  1:4404676/1‑100 (MQ=255)
                                                                  aCAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACa                                  <  2:1050937/100‑1 (MQ=255)
                                                                    aaTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgg                                <  2:1989022/100‑1 (MQ=255)
                                                                     aTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                               <  1:3928159/100‑1 (MQ=255)
                                                                      ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                              <  2:1744307/100‑1 (MQ=255)
                                                                      ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                              >  2:3681140/1‑100 (MQ=255)
                                                                       tATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                              <  1:2902521/99‑1 (MQ=255)
                                                                       atctgggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGc                             <  2:4095427/96‑1 (MQ=255)
                                                                        aTCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                            <  2:2801781/100‑1 (MQ=255)
                                                                        aTCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                            <  1:2384942/100‑1 (MQ=255)
                                                                         tCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCaa                           <  2:3535563/100‑1 (MQ=255)
                                                                         tCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTACTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                            >  1:3094807/1‑99 (MQ=255)
                                                                           gggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCagat                         >  1:4095398/1‑97 (MQ=255)
                                                                             gCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGt                       <  2:2024196/100‑1 (MQ=255)
                                                                              cGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTa                      <  2:326203/100‑1 (MQ=255)
                                                                                 aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                   <  2:1986548/100‑1 (MQ=255)
                                                                                  acaacttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACt                  >  2:1897152/8‑100 (MQ=38)
                                                                                   aGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTaaattc                 >  2:426293/1‑97 (MQ=255)
                                                                                    ggCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCg                <  1:1658428/100‑1 (MQ=255)
                                                                                    ggCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCg                >  2:20370/1‑100 (MQ=255)
                                                                                    ggCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCg                >  2:610643/1‑100 (MQ=255)
                                                                                     gCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGAGCAAGCGTAAACTc                 >  2:4046045/1‑98 (MQ=255)
                                                                                     acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa               >  2:3182006/5‑100 (MQ=39)
                                                                                     acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa               >  2:1806525/5‑100 (MQ=39)
                                                                                        taCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGc            >  1:125149/2‑100 (MQ=255)
                                                                                         aCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa           >  1:4262854/1‑100 (MQ=255)
                                                                                         aCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa           >  1:1005381/1‑100 (MQ=255)
                                                                                          ccGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGGCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          <  1:1058551/100‑1 (MQ=255)
                                                                                          ccGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          <  1:1366474/100‑1 (MQ=255)
                                                                                          ccGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          >  2:1454967/1‑100 (MQ=255)
                                                                                          ccGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          >  2:2567666/1‑100 (MQ=255)
                                                                                           cGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCaga         <  1:3434707/100‑1 (MQ=255)
                                                                                               cATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTCCAGGGCCAAACGTAAACTCGACGCattgca      >  2:264986/1‑94 (MQ=255)
                                                                                               cATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa     <  1:599985/100‑1 (MQ=255)
                                                                                                 tttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGc   <  1:2586316/100‑1 (MQ=255)
                                                                                                  ttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGc   <  2:4379977/99‑1 (MQ=255)
                                                                                                  |                                                                                                   
AGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCC  >  NC_000913/1233056‑1233253

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