Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,157 T→A 8.3% intergenic (‑127/+19) dsbB ← / ← nhaB protein thiol:quinone oxidoreductase DsbB/Na(+):H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,233,1570TA 8.3% 201.1 / 4.6 67intergenic (‑127/+19)dsbB/nhaBprotein thiol:quinone oxidoreductase DsbB/Na(+):H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base A (4/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGC  >  NC_000913/1233009‑1233300
                                                                                                                                                    |                                                                                                                                               
aTTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                                                                <  2:1376409/150‑1 (MQ=255)
aTTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                                                                <  1:179672/150‑1 (MQ=255)
   ggTTCNAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGc                                                                                                                                             >  1:3/1‑150 (MQ=255)
     ttCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGccc                                                                                                                                           <  2:1242383/150‑1 (MQ=255)
       cAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                                                                         <  1:930306/150‑1 (MQ=255)
       cAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCg                                                                                                                                          <  1:1118333/149‑1 (MQ=255)
       cAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCg                                                                                                                                          >  2:1118333/1‑149 (MQ=255)
        aaaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGt                                                                                                                                        >  2:904242/1‑150 (MQ=255)
         aaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                                                                       <  1:947699/150‑1 (MQ=255)
             tCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgtt                                                                                                                                   >  1:863894/1‑150 (MQ=255)
               gCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgt                                                                                                                                 <  1:326080/150‑1 (MQ=255)
                cAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgtt                                                                                                                                <  1:1177050/150‑1 (MQ=255)
                 aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                                                               <  2:321578/150‑1 (MQ=255)
                    aTAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCa                                                                                                                            <  1:1368059/150‑1 (MQ=255)
                       aTCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAg                                                                                                                         <  1:793751/150‑1 (MQ=255)
                       aTCATTTCCCTGCAATATGCATAGAGTGCCCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAg                                                                                                                         >  2:1405173/1‑150 (MQ=255)
                        tCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGc                                                                                                                        >  1:443628/1‑150 (MQ=255)
                         cATTTCCCTGCAATATGCATAGAGTGCTCAGTTTAAACCAATTCATTCGCTGAGACACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATCATCGGGCGAAAAGGCAAACCGGGCATTTATAAGACCGGTAAGTTGTTATCAAAGCg                                                                                                                       >  2:542716/1‑150 (MQ=255)
                            ttCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                                                                    <  2:1281765/150‑1 (MQ=255)
                             tCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGATGCATCAAAAATAGGACGCAATTTGACATTTATCTTGCGATAAGGCAAAGCGGTCATTTATATTCGCGGTAAGTTGTTATAAATGAGTTGc                                                                                                                   <  2:69279/150‑1 (MQ=255)
                             tCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                                                                   <  2:38661/150‑1 (MQ=255)
                                cTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTAt                                                                                                                <  2:1182941/150‑1 (MQ=255)
                                 tGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATc                                                                                                               >  1:209358/1‑150 (MQ=255)
                                  gCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                                                                              <  2:1129007/150‑1 (MQ=255)
                                    aaTATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACAGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAg                                                                                                            >  1:1289768/1‑150 (MQ=255)
                                     atatGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGc                                                                                                           <  2:1201238/150‑1 (MQ=255)
                                          cATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAtt                                                                                                      >  1:458101/1‑150 (MQ=255)
                                          cATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAtt                                                                                                      >  1:1079339/1‑150 (MQ=255)
                                              gagTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTGTATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                                                                  <  1:937737/150‑1 (MQ=255)
                                              gagTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                                                                  >  2:631097/1‑150 (MQ=255)
                                              gagTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                                                                  <  1:978901/150‑1 (MQ=255)
                                                  gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                                                                              >  1:1174547/1‑150 (MQ=255)
                                                    gcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAcc                                                                                            <  1:1405173/150‑1 (MQ=255)
                                                       gTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAtt                                                                                         >  1:983148/1‑150 (MQ=255)
                                                         ttAAACCGATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCg                                                                                       <  1:784720/150‑1 (MQ=255)
                                                         ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCg                                                                                       <  1:670252/150‑1 (MQ=255)
                                                            aaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGtt                                                                                    >  2:1316200/1‑150 (MQ=255)
                                                              ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTAc                                                                                  >  1:563160/1‑150 (MQ=255)
                                                               cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATATGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACa                                                                                 <  1:813339/150‑1 (MQ=255)
                                                               cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACa                                                                                 <  2:1079339/150‑1 (MQ=255)
                                                                    cattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGc                                                                            <  2:1174547/150‑1 (MQ=255)
                                                                        cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGc                                                                        >  1:516240/1‑150 (MQ=255)
                                                                         gCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGTTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCg                                                                       <  2:68429/150‑1 (MQ=255)
                                                                          cTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                                                                              <  2:958336/142‑1 (MQ=255)
                                                                          cTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                                                                              >  1:958336/1‑142 (MQ=255)
                                                                                     ggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGc                                                           <  2:422765/150‑1 (MQ=255)
                                                                                        ggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCaga                                                        <  2:3/150‑1 (MQ=255)
                                                                                        ggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCaga                                                        <  1:205636/150‑1 (MQ=255)
                                                                                          cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGc                                                      <  2:863894/150‑1 (MQ=255)
                                                                                            gcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa                                                    >  1:712087/1‑150 (MQ=255)
                                                                                              aaaaaTAGGACGCAATTGGACAATTATCTGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGc                                                  <  2:443628/150‑1 (MQ=255)
                                                                                                          cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCa                                      >  1:1310387/1‑150 (MQ=255)
                                                                                                             ttGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGa                                   <  2:458101/150‑1 (MQ=255)
                                                                                                               ggACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAATCCGACGAGTGTCAGGACg                                 >  2:556332/1‑150 (MQ=255)
                                                                                                                 aCAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAg                               <  2:563160/150‑1 (MQ=255)
                                                                                                                  cAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAgg                              >  1:894749/1‑150 (MQ=255)
                                                                                                                       aTCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTa                         >  1:90186/1‑150 (MQ=255)
                                                                                                                         cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAg                       <  2:666647/150‑1 (MQ=255)
                                                                                                                         cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAg                       <  2:490210/150‑1 (MQ=255)
                                                                                                                         cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTCGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAg                       <  1:304081/150‑1 (MQ=255)
                                                                                                                          gggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAgg                      >  1:1202871/1‑150 (MQ=255)
                                                                                                                           ggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGc                     <  1:34712/150‑1 (MQ=255)
                                                                                                                                aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGc                <  1:904242/150‑1 (MQ=255)
                                                                                                                                    acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATc            >  2:1415841/5‑150 (MQ=255)
                                                                                                                                     cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATcc           >  1:176813/1‑150 (MQ=255)
                                                                                                                                      ttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCca          >  2:1208446/3‑150 (MQ=255)
                                                                                                                                       aaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCcac         >  2:147696/1‑150 (MQ=255)
                                                                                                                                        aCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCcaca        >  1:489090/1‑150 (MQ=255)
                                                                                                                                              cATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTATACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGc  <  1:13779/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGC  >  NC_000913/1233009‑1233300

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