Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,390,012 Δ1 bp 6.0% intergenic (+48/+36) yfbP → / ← nuoN TPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,390,0120G. 6.0% 232.1 / 8.4 67intergenic (+48/+36)yfbP/nuoNTPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base G (28/35);  new base . (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

ATTTCAGGAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCT  >  NC_000913/2389865‑2390155
                                                                                                                                                   |                                                                                                                                                
aTTTCAGGAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGc                                                                                                                                                 <  4:224277/149‑1 (MQ=255)
 tttCAGGAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCa                                                                                                                                                >  3:256997/1‑149 (MQ=255)
  ttCAGGAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCAt                                                                                                                                               >  3:381351/1‑149 (MQ=255)
            cTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCt                                                                                                                                    <  1:654656/149‑1 (MQ=255)
              aTGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGGGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCttt                                                                                                                                  >  3:528312/1‑149 (MQ=255)
               tGATTTTAGTATCGATGAATGGCGGGCAAAATTCAGGAAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTc                                                                                                                                 <  4:274162/149‑1 (MQ=255)
               tGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTc                                                                                                                                 <  1:624767/149‑1 (MQ=255)
                gATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTCGGCGGtttttc                                                                                                                                <  2:527765/149‑3 (MQ=255)
                   tttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGAc                                                                                                                             <  4:664985/149‑1 (MQ=255)
                   tttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGAc                                                                                                                             <  2:91153/149‑1 (MQ=255)
                   tttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGAc                                                                                                                             <  2:205877/149‑1 (MQ=255)
                   tttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGAc                                                                                                                             <  2:1283/149‑1 (MQ=255)
                    ttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACt                                                                                                                            <  3:188071/149‑1 (MQ=255)
                          tCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACaaa                                                                                                                      <  4:416901/149‑1 (MQ=255)
                            gATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGt                                                                                                                    <  2:622945/149‑1 (MQ=255)
                              tGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTaa                                                                                                                  <  1:364722/149‑1 (MQ=255)
                                aaTGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACa                                                                                                                <  2:63986/149‑1 (MQ=255)
                                aaTGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACa                                                                                                                <  1:38673/149‑1 (MQ=255)
                                aaTGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACa                                                                                                                <  1:343010/149‑1 (MQ=255)
                                aaTGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACa                                                                                                                <  4:639821/149‑1 (MQ=255)
                                  tGGCGGGCAAAATTCAGGCAGTTATGTTTAGGGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGa                                                                                                              <  4:528312/149‑1 (MQ=255)
                                   ggCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGATCACTTAACGCGCGGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTTACTTACAAAGTAACAGAt                                                                                                             <  4:282264/149‑1 (MQ=255)
                                            aaTTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGc                                                                                                    >  2:444442/1‑149 (MQ=255)
                                             aTTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCg                                                                                                   >  3:188798/1‑149 (MQ=255)
                                              ttCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCgg                                                                                                  >  4:180636/1‑149 (MQ=255)
                                                     aGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGcc                                                                                           >  1:17743/1‑149 (MQ=255)
                                                                   tATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCt                                                                             <  2:107820/149‑1 (MQ=255)
                                                                   tATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCt                                                                             <  1:399371/149‑1 (MQ=255)
                                                                     ttttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTaa                                                                           >  4:250965/1‑149 (MQ=255)
                                                                           cGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCg                                                                     <  3:25071/149‑1 (MQ=255)
                                                                             aTAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTcggcggtt                                                                                                                                    >  1:649896/1‑82 (MQ=255)
                                                                             aTAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTcggcggtt                                                                                                                                    <  2:649896/84‑3 (MQ=38)
                                                                                gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCATCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGAGTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGTTAATCCGCGGTtgt                                                                >  4:420258/1‑149 (MQ=255)
                                                                                   gTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                                                             >  2:596729/1‑149 (MQ=255)
                                                                                       acaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATa                                                         >  3:311886/1‑149 (MQ=255)
                                                                                         aGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATaca                                                       <  1:145146/149‑1 (MQ=255)
                                                                                          gTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TTGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac                                                      >  4:305564/1‑149 (MQ=255)
                                                                                          gTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac                                                      >  4:587077/1‑149 (MQ=255)
                                                                                           ttCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAcc                                                     >  1:484651/1‑149 (MQ=255)
                                                                                               cTTAACGCTCTGCGTCTATTTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc                                                 >  4:350130/1‑149 (MQ=255)
                                                                                                ttAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCa                                                >  3:295414/1‑149 (MQ=255)
                                                                                                ttAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCa                                                >  4:547998/1‑149 (MQ=255)
                                                                                                 tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAc                                               <  1:588928/149‑1 (MQ=255)
                                                                                                   aCGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCa                                             <  2:241197/149‑1 (MQ=255)
                                                                                                          gCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCaa                                      >  4:364243/1‑149 (MQ=255)
                                                                                                            gTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACa                                    >  4:642482/1‑149 (MQ=255)
                                                                                                              ctatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGt                                  <  1:32482/149‑1 (MQ=255)
                                                                                                                 tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCa                               <  4:453611/149‑1 (MQ=255)
                                                                                                                  ctactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCag                              >  2:417424/1‑149 (MQ=255)
                                                                                                                   tactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATCCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCaga                             <  2:207916/149‑1 (MQ=255)
                                                                                                                      taTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGAt                          >  2:124298/1‑149 (MQ=255)
                                                                                                                       aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                                                             <  2:136402/113‑1 (MQ=255)
                                                                                                                       aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                                                             >  1:136402/1‑113 (MQ=255)
                                                                                                                       aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATc                         <  1:417424/149‑1 (MQ=255)
                                                                                                                       aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATc                         <  1:326049/149‑1 (MQ=255)
                                                                                                                        tCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCa                        >  2:570283/1‑149 (MQ=255)
                                                                                                                        tCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCa                        <  2:17743/149‑1 (MQ=255)
                                                                                                                        ccccTTCTTCTTCCGAAACGCCGCCCAGAACG‑TCGGAGGCCTTCTGACTTACAAAGTAACCGATTACATCAGCGAAATTGTCAAACGCACAATGCTAATCAGCGATTGTGGTCATACACCCAGCAACACTACCATTACTGCACAGATCa                        <  2:118726/148‑1 (MQ=255)
                                                                                                                          ccTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGGCTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGc                      >  2:569123/1‑149 (MQ=255)
                                                                                                                               tttttCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCaaa                                                                                        <  2:120387/78‑1 (MQ=255)
                                                                                                                               tttttCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCaaa                                                                                        >  1:120387/1‑78 (MQ=255)
                                                                                                                               tttttCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAAc                 >  3:80973/1‑149 (MQ=255)
                                                                                                                               tttttCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAAc                 <  3:547998/149‑1 (MQ=255)
                                                                                                                                    tcagaaagCCGCCGA‑CATGCTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAc            >  2:393615/5‑149 (MQ=255)
                                                                                                                                      gAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccg          >  1:627549/1‑149 (MQ=255)
                                                                                                                                         aaaCCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccgccc       >  4:2307/1‑149 (MQ=255)
                                                                                                                                              gccgAGCATG‑TCGGCTGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCt  <  4:276971/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATTTCAGGAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCT  >  NC_000913/2389865‑2390155

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑