Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,390,016:1 +C 6.9% intergenic (+52/+32) 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,0161.C 6.9% 491.3 / 32.1 116intergenic (+52/+32)yfbP/nuoNTPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base . (51/57);  new base C (3/5);  total (54/62)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 116 total aligned reads displayed.

GAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTGTAC  >  NC_000913/2389872‑2390159
                                                                                                                                                 |                                                                                                                                                
gAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGTATGTcgg                                                                                                                                              >  4:210944/1‑149 (MQ=255)
gAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcgg                                                                                                                                              >  2:62421/1‑149 (MQ=255)
gAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTc                                                                                                                                                >  1:626429/1‑147 (MQ=255)
gAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTc                                                                                                                                                <  2:626429/147‑1 (MQ=255)
 aaCGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcggc                                                                                                                                             <  4:823951/149‑1 (MQ=255)
   cGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTcggcgg                                                                                                                                           <  1:405721/149‑1 (MQ=255)
   cGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTcggcgg                                                                                                                                           <  2:356186/149‑1 (MQ=255)
   cGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTcggcgg                                                                                                                                           <  1:550438/149‑1 (MQ=255)
    gCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCggcggc                                                                                                                                         >  4:798253/1‑149 (MQ=255)
    gCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTcggcggt                                                                                                                                          <  3:177015/149‑2 (MQ=255)
     cTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTcggcggtt                                                                                                                                         <  3:262638/149‑3 (MQ=255)
      tATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCtt                                                                                                                                       <  4:724686/149‑1 (MQ=255)
      tATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCtt                                                                                                                                       <  3:430883/149‑1 (MQ=255)
      tATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCtt                                                                                                                                       >  1:630487/1‑149 (MQ=255)
           ttttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAGAAACCGCCGAGCATGTCGGCGGCTTTCTGa                                                                                                                                  >  1:580784/1‑149 (MQ=255)
           ttttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGa                                                                                                                                  >  2:253043/1‑149 (MQ=255)
           ttttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGa                                                                                                                                  <  4:13023/149‑1 (MQ=255)
             ttAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACt                                                                                                                                <  3:365155/149‑1 (MQ=255)
                atatCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTAc                                                                                                                             <  3:346928/149‑1 (MQ=255)
                 tatCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACa                                                                                                                            >  2:219032/1‑149 (MQ=255)
                 gagaactGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACa                                                                                                                            <  3:147419/143‑1 (MQ=255)
                   tCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCTCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACaaa                                                                                                                          <  3:77861/149‑1 (MQ=255)
                   tCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACaaa                                                                                                                          <  3:117961/149‑1 (MQ=255)
                   tCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACaaa                                                                                                                          <  2:848070/149‑1 (MQ=255)
                     gATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGt                                                                                                                        >  1:687208/1‑149 (MQ=255)
                     gATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGt                                                                                                                        <  1:219032/149‑1 (MQ=255)
                       tGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTaa                                                                                                                      <  3:798253/149‑1 (MQ=255)
                         aaTGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACa                                                                                                                    <  1:597177/149‑1 (MQ=255)
                         aaTGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACa                                                                                                                    <  4:395485/149‑1 (MQ=255)
                         aaTGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACa                                                                                                                    <  4:42032/149‑1 (MQ=255)
                         aaTGGCGGGAAAAATTAAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGAACACAGTTCTCTTAACGCTCTGCGTCTATCTACCATCCCTTATTTTACAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGCCTTATAAAGAAACa                                                                                                                    <  2:295882/149‑1 (MQ=255)
                            ggCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAt                                                                                                                 <  3:388731/149‑1 (MQ=255)
                             gCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAtt                                                                                                                <  2:670457/149‑1 (MQ=255)
                             gCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAtt                                                                                                                <  2:630487/149‑1 (MQ=255)
                              cGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTa                                                                                                               <  4:194150/149‑1 (MQ=255)
                                ggCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACa                                                                                                             >  2:19979/1‑149 (MQ=255)
                                           ggCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGa                                                                                                                  >  3:31585/1‑133 (MQ=255)
                                           ggCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGa                                                                                                                  <  4:31585/133‑1 (MQ=255)
                                           ggCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATTAGCGGCAtt                                                                                                  <  4:334364/149‑1 (MQ=255)
                                           ggCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCAtt                                                                                                  <  2:541088/149‑1 (MQ=255)
                                           ggCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCAtt                                                                                                  <  4:707258/149‑1 (MQ=255)
                                            gCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTg                                                                                                 >  2:425523/1‑149 (MQ=255)
                                                   tGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACg                                                                                          >  2:832465/1‑149 (MQ=255)
                                                     tttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGca                                                                                        >  1:115549/1‑149 (MQ=255)
                                                      ttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTTCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGcac                                                                                       >  2:143613/1‑149 (MQ=255)
                                                           gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGc                                                                                  >  1:180099/1‑149 (MQ=255)
                                                                         gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGTATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                                                                    <  3:210944/149‑1 (MQ=255)
                                                                                acaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATa                                                             <  3:111332/149‑1 (MQ=255)
                                                                                acaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATa                                                             >  2:743554/1‑149 (MQ=255)
                                                                                acaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATa                                                             <  1:199334/149‑1 (MQ=255)
                                                                                  aGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATaca                                                           >  3:297110/1‑149 (MQ=255)
                                                                                   gTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac                                                          <  1:345285/149‑1 (MQ=255)
                                                                                    ttCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAcc                                                         >  2:832849/1‑149 (MQ=255)
                                                                                      cACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCa                                                       >  3:877380/1‑149 (MQ=255)
                                                                                      cACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCa                                                       >  1:712375/1‑149 (MQ=255)
                                                                                       aCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCAtt                                                                                                  <  1:173525/105‑1 (MQ=255)
                                                                                       aCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCAtt                                                                                                  >  2:173525/1‑105 (MQ=255)
                                                                                        cTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc                                                     >  3:557402/1‑149 (MQ=255)
                                                                                        cTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc                                                     >  4:817834/1‑149 (MQ=255)
                                                                                        cTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc                                                     >  2:509841/1‑149 (MQ=255)
                                                                                        cTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc                                                     >  3:269626/1‑149 (MQ=255)
                                                                                         ttAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTa                                                                                                               <  1:222944/90‑1 (MQ=255)
                                                                                          tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAc                                                   <  1:19979/149‑1 (MQ=255)
                                                                                              gCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGTATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGt                                               <  1:577410/149‑1 (MQ=255)
                                                                                                tctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTAc                                             >  1:419201/1‑149 (MQ=255)
                                                                                                        tatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTg                                     >  2:411502/1‑149 (MQ=255)
                                                                                                        tatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTg                                     >  3:671471/1‑149 (MQ=255)
                                                                                                          tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTa                                              >  3:231385/1‑138 (MQ=255)
                                                                                                          tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTa                                              <  4:231385/138‑1 (MQ=255)
                                                                                                          tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCa                                   <  2:712375/149‑1 (MQ=255)
                                                                                                            tactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCaga                                 <  3:214614/149‑1 (MQ=255)
                                                                                                            tactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCaga                                 >  1:53072/1‑149 (MQ=255)
                                                                                                             actaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCagag                                <  3:817834/149‑1 (MQ=255)
                                                                                                              ctaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAgaga                               >  3:392136/1‑149 (MQ=255)
                                                                                                               taTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAccc                                                        >  1:275573/1‑123 (MQ=255)
                                                                                                               taTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAccc                                                        <  2:275573/123‑1 (MQ=255)
                                                                                                               taTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGAt                              >  4:300777/1‑149 (MQ=255)
                                                                                                                 tCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCagag                                >  3:591805/1‑145 (MQ=255)
                                                                                                                 tCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCCGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCatag                                <  4:591805/145‑1 (MQ=255)
                                                                                                                    cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTaa                                                                               >  3:491615/1‑95 (MQ=255)
                                                                                                                    cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTaa                                                                               <  4:491615/95‑1 (MQ=255)
                                                                                                                    cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTAc                                             <  3:492918/129‑1 (MQ=255)
                                                                                                                    cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTAc                                             >  4:492918/1‑129 (MQ=255)
                                                                                                                    cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGca                         <  2:580784/149‑1 (MQ=255)
                                                                                                                    cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGca                         >  3:45443/1‑149 (MQ=255)
                                                                                                                     ttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTCACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGc                          <  2:255110/147‑1 (MQ=255)
                                                                                                                     ttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGcac                        >  3:695269/1‑149 (MQ=255)
                                                                                                                        tttttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAAc                     <  1:793368/149‑1 (MQ=255)
                                                                                                                        tttttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAAc                     <  1:832849/149‑1 (MQ=255)
                                                                                                                         ttttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACg                    >  1:46393/1‑149 (MQ=255)
                                                                                                                            tCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGtt                                                                      >  4:411303/1‑96 (MQ=255)
                                                                                                                            gtcagaaagCCGCCGACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGAt                  >  4:79932/6‑149 (MQ=38)
                                                                                                                              aGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAcc               >  1:583846/1‑149 (MQ=255)
                                                                                                                               gAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccg              >  1:6833/1‑149 (MQ=255)
                                                                                                                               gAAAAACCGCCAAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccg              >  1:691070/1‑149 (MQ=255)
                                                                                                                                gaaagCCGCCGACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccg              >  1:879354/2‑149 (MQ=255)
                                                                                                                                   agccgccgACATGCTCGGCGGTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccgcccg          >  2:507/3‑149 (MQ=255)
                                                                                                                                     ccgccgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCcgcg        <  1:598227/149‑1 (MQ=255)
                                                                                                                                      cgccgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCgcgc       <  2:687208/149‑1 (MQ=255)
                                                                                                                                           aGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTGTAc  <  2:53072/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTGTAC  >  NC_000913/2389872‑2390159

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