Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,157 T→A 5.9% 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 5.9% 201.8 / 4.3 68intergenic (‑127/+19)dsbB/nhaBprotein thiol:quinone oxidoreductase DsbB/Na(+):H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base A (2/2);  total (37/31)
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.03e-01

TCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACAC  >  NC_000913/1233015‑1233295
                                                                                                                                              |                                                                                                                                          
tCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCCCCAGGAGGCAGCAAAAATAGGCCGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtat                                                                                                                                            >  1:808208/1‑143 (MQ=255)
  aaaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGCGGCAGCAAAAATAGGACGCAATTCGACAATCATCGGGCGACAAGGCAAACCGGGCCTTCatat                                                                                                                                          >  2:350056/1‑143 (MQ=255)
       tCGCCACATAATGATGTCCCTGCAATAGGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCg                                                                                                                                     <  2:194450/143‑1 (MQ=255)
       tCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCg                                                                                                                                     <  1:720038/143‑1 (MQ=255)
       tCGCAACAGAAGCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCg                                                                                                                                     <  1:106635/143‑1 (MQ=255)
          cAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                                                                  >  1:208457/1‑143 (MQ=255)
          cAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtt                                                                                                                                  <  2:564946/143‑2 (MQ=255)
           aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGttt                                                                                                                                 <  1:318054/143‑3 (MQ=255)
                   cATTTCCCTTCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTAt                                                                                                                         <  1:219391/143‑1 (MQ=255)
                      ttCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGAGCCACCAGGAGGTAGGAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaa                                                                                                                      >  1:1047669/1‑143 (MQ=255)
                       tCCCTGCAATATGCATATAGTTCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                                                                     <  1:463394/143‑1 (MQ=255)
                       gcctggCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                                                                     <  2:56841/138‑1 (MQ=255)
                       gcccTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                                                                     <  1:1039580/142‑1 (MQ=255)
                         ccTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGc                                                                                                                   >  2:1062902/1‑143 (MQ=255)
                          cTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAACAATAGGACGCAATTGGACAATTATCGGGCGAAAATGCAAACCGGGCATTTATATGCCCGGTACGTTGTTATCACAGCg                                                                                                                  >  2:75799/1‑143 (MQ=255)
                             cAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCTTCATGTGGGGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                                                               >  1:212005/1‑143 (MQ=255)
                                tatGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTa                                                                                                            >  1:246248/1‑143 (MQ=255)
                                     aTAGAGTGCGCAGTTTAAACCAATGCATGCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAg                                                                                                       <  2:739370/143‑1 (MQ=255)
                                        gagTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATGCAGccc                                                                                                    >  2:585074/1‑143 (MQ=255)
                                            gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAttt                                                                                                <  1:492887/143‑1 (MQ=255)
                                            gcgcAGGTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAttt                                                                                                <  1:133462/143‑1 (MQ=255)
                                               cAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCACAGCGTTGCTATCCAGCCCATTTGCa                                                                                             >  1:1062470/1‑143 (MQ=255)
                                                  tttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaa                                                                                          >  1:1046541/1‑143 (MQ=255)
                                                  tttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaa                                                                                          >  1:890662/1‑143 (MQ=255)
                                                      aaCCAATTCATTCGCTGAGCCACCAGGAGGCAGCAAAAATAGGACGCAATTAGACAATGATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCa                                                                                      <  2:287823/143‑1 (MQ=255)
                                                         cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                                                                   >  1:313471/1‑143 (MQ=255)
                                                         cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                                                                   >  1:1076541/1‑143 (MQ=255)
                                                           attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCgg                                                                                 >  1:1023039/1‑143 (MQ=255)
                                                              cattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAACCGGTGCTATCAGGCGCATTTGCATAAACCGTTCGGCTa                                                                              >  2:702437/1‑143 (MQ=255)
                                                               attcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTAc                                                                             >  1:97963/1‑143 (MQ=255)
                                                                ttcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACa                                                                            >  1:772032/1‑143 (MQ=255)
                                                                ttcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACa                                                                            >  1:772028/1‑143 (MQ=255)
                                                                 tcGCGGCGCCAGCAGGAGGCAGCAAAAAGAGGACGCAATTGGACAATGATCGGGCGAAAAGGCAAACCGGGCATTTAGATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATGTGCATAAACCATTCGGTTACAg                                                                           <  2:124807/143‑1 (MQ=255)
                                                                   gCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                                                                         <  2:167707/143‑1 (MQ=255)
                                                                    cTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTACGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCAGAAACCATAAGGTTGCAgcgg                                                                        >  2:947621/1‑143 (MQ=255)
                                                                       cgcgACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCaa                                                                     <  2:725943/143‑1 (MQ=255)
                                                                       cgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTAGCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCaa                                                                     <  1:830719/143‑1 (MQ=255)
                                                                        gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCGAAGGGTTGCTATCCAGCCCATTTGCATAAACCATGCGGTTGCAGGGGGAAg                                                                    >  1:225549/1‑143 (MQ=255)
                                                                         cgagcaGGAGGCAGCAAGAATAGGACGCAAGTGGACAAGTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGc                                                                   <  2:1049491/139‑1 (MQ=255)
                                                                            ccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTa                                                                >  2:919803/1‑143 (MQ=255)
                                                                               ggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                             <  2:212005/143‑1 (MQ=255)
                                                                                gaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACt                                                            >  1:925207/1‑143 (MQ=255)
                                                                                   gcagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa                                                         >  2:790409/1‑143 (MQ=255)
                                                                                      gcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGtt                                                                               <  2:749661/118‑1 (MQ=255)
                                                                                      gcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGtt                                                                               >  1:749661/1‑118 (MQ=255)
                                                                                       caAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa                                                     >  1:1021076/1‑143 (MQ=255)
                                                                                        aagaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag                                                    <  1:439697/140‑1 (MQ=255)
                                                                                        aaaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGCCGCag                                                    >  1:775917/1‑143 (MQ=255)
                                                                                        aaaaaTAGGACGCAATTGAACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag                                                    <  1:3837/143‑1 (MQ=255)
                                                                                         aaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCaga                                                   >  1:79142/1‑143 (MQ=255)
                                                                                             tAGGATGCATTTGGATAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa                                               <  1:924006/143‑1 (MQ=255)
                                                                                             tAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa                                               <  1:751913/143‑1 (MQ=255)
                                                                                                     aaTTGGACACTTATCGGGCGCCAAGGCCAACCGGGCATTTATCTGCCCGGTAAGTTGCTATCAAAGCGTTGCTATCCAGCCCATTTGCATCAACCACTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCcgacga                                       <  2:630343/143‑1 (MQ=255)
                                                                                                          gACAAATATCGGCCGAAAAGGCAAACCTGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTc                                  <  1:885238/143‑1 (MQ=255)
                                                                                                             aaTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATACACCATTCGGTTACAGGGGCAAGCGTAAACTCGACACAGAGCCAGCCGACGAGTGTCAgg                               >  1:574369/1‑143 (MQ=255)
                                                                                                             aaTTATCGGGCGAAAAGGCAAACAGGGCCTTTATATGCCCGCGAAGTTGTTATCAAAGCGTTGCTATCCAGCCCACTTGCATCAACCATTCGGCTACAGGGGCCAGCGCAGACTCGCCGCCGACCCAGCCGACGAGTGTCAgc                               >  2:370859/1‑142 (MQ=255)
                                                                                                                   cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAgg                         <  1:328896/143‑1 (MQ=255)
                                                                                                                        gAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTa                    <  1:919803/143‑1 (MQ=255)
                                                                                                                        gAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTa                    >  1:1100070/1‑143 (MQ=255)
                                                                                                                          agaGGCAAACCGTGCATTTATATGCCCGGTAAGTTGTTAACAAAGCGGTGCAAGCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGCGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAg                  <  1:396406/141‑1 (MQ=255)
                                                                                                                               cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACTATTCGGTTACAGGGGCAAGCGTAAACTCGAAGCAGAGCAAGCCGACGAGTGGCAGGAAGAAGGTGTAAGGCAgg             >  2:797826/1‑143 (MQ=255)
                                                                                                                               cAAACCGGGCATTTATATGCCCGGTAAGTTGTGATCAAAGCGTTGCTATCCAGCCCAGGTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAgg             <  1:896950/143‑1 (MQ=255)
                                                                                                                                 aaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGATACAGGGGCAAGCGTAAACTCGAAGCAGATCAAACAGACGAGTGTCAGGACGAGGGAGTAAGGCAGGGc           >  2:952663/1‑143 (MQ=255)
                                                                                                                                    cGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCAt        >  1:531266/1‑143 (MQ=255)
                                                                                                                                    cGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGCGTCCGGACGAGGGTGTCAGGCAGGGCCAt        >  2:747951/1‑143 (MQ=255)
                                                                                                                                      ggCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGCCGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATcc      >  2:226686/1‑143 (MQ=255)
                                                                                                                                       gCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGAGCAAGCCTAAACTCGACGCAGAGCAAGCCGACGAGTGACAGGACGAGGGTGTAAGGCAGAGCCATCca     >  2:388205/1‑143 (MQ=255)
                                                                                                                                       gCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCAAAAACAGAGGCAAGCGTAAAATCGACGCCAAGCAAGCCGACGAGAGTCAGGACGAGGGTATACGCAAGGGcccacca     >  2:1017186/1‑138 (MQ=255)
                                                                                                                                        cATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATGCGGTTACAGGGGCAAGCGGAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCcac    <  2:1046541/143‑1 (MQ=255)
                                                                                                                                          tttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGAGGCCAGCGTAAACTCGACGCCGAGCAAGCCGCCGAGTGTCAGGACGAGGGCGTAAGGCAGAGCCCTCCacac  >  2:601164/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                          
TCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACAC  >  NC_000913/1233015‑1233295

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: