Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,313,434 T→C intergenic (+258/‑54) micF → / → rcsD sRNA antisense regulator affecting ompF expression, Hfq‑dependent/phosphotransfer intermediate protein in two‑component regulatory system with RcsBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,313,4340TC94.7% 258.5 / 2.5 78intergenic (+258/‑54)micF/rcsDsRNA antisense regulator affecting ompF expression, Hfq‑dependent/phosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (39/33);  total (42/36)
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.93e-01

TGGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGTACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATC  >  NC_000913/2313306‑2313568
                                                                                                                                |                                                                                                                                      
tgGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTtata                                                                                                                              <  2:196055/139‑1 (MQ=255)
tgGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTtata                                                                                                                              <  1:719523/139‑1 (MQ=255)
   gCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGTACCCTTTATACTg                                                                                                                           >  1:1351546/1‑139 (MQ=255)
                 aTTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTcgc                                                                                                             >  1:905022/1‑139 (MQ=255)
                           cAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGTACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTaa                                                                                                   <  2:1351546/139‑1 (MQ=255)
                                tGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAAc                                                                                                  <  1:1218740/135‑1 (MQ=255)
                                tGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAAc                                                                                                  >  2:1218740/1‑135 (MQ=255)
                                tGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGt                                                                                              >  1:1568251/1‑139 (MQ=255)
                                     cTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTta                                                                                                                                >  1:1226315/1‑100 (MQ=255)
                                     cTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTta                                                                                                                                <  2:1226315/100‑1 (MQ=255)
                                     cTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATaa                                                                                         >  1:743343/1‑139 (MQ=255)
                                       tGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAac                                                                                       >  1:1048966/1‑139 (MQ=255)
                                         tACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAacac                                                                                     >  1:161334/1‑139 (MQ=255)
                                           cATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGa                                                                                   <  2:743343/139‑1 (MQ=255)
                                            aTCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGTACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGaa                                                                                  <  1:1249774/139‑1 (MQ=255)
                                                 tttCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGt                                                                             <  1:830980/139‑1 (MQ=255)
                                                  ttCCCCCAATGCAGGATGATAAATATCATGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTc                                                                            >  2:428289/1‑139 (MQ=255)
                                                        cAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAag                                                                      >  1:1476090/1‑139 (MQ=255)
                                                           tGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAgaga                                                                   >  1:1366013/1‑139 (MQ=255)
                                                           tGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAgaga                                                                   >  2:1544832/1‑139 (MQ=255)
                                                            gCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGAc                                                                  <  1:1544832/139‑1 (MQ=255)
                                                            gCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGAc                                                                  <  2:85367/139‑1 (MQ=255)
                                                              aGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACaa                                                                <  2:1476090/139‑1 (MQ=255)
                                                                     tAAATATCACGGGAGAATAGAGAATCATCAATCATGTAAGAGTCTGGAAGTTCACACTG‑‑CCTTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGa                                                       >  1:1065884/1‑139 (MQ=255)
                                                                     tAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCAc                                                         >  2:1214124/1‑139 (MQ=255)
                                                                     tAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCAc                                                         >  1:664438/1‑139 (MQ=255)
                                                                     tAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAA‑‑ACGAATGCGTCAGAAAGAGACAACGGCCACGa                                                       >  2:1506303/1‑139 (MQ=255)
                                                                          atCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGAccc                                                    <  2:1366013/139‑1 (MQ=255)
                                                                            cACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTcgc                                                                                                             <  2:676744/80‑1 (MQ=255)
                                                                            cACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTcgc                                                                                                             >  1:676744/1‑80 (MQ=255)
                                                                                        gagaATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTa                                      <  2:823945/139‑1 (MQ=255)
                                                                                         agaATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTAc                                     >  1:866361/1‑139 (MQ=255)
                                                                                         agaATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTAc                                     >  1:629956/1‑139 (MQ=255)
                                                                                         agaATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTAc                                     >  2:728543/1‑139 (MQ=255)
                                                                                           aaTCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTc                                                                                             <  2:1477086/81‑1 (MQ=255)
                                                                                           aaTCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTc                                                                                             >  1:1477086/1‑81 (MQ=255)
                                                                                           aaTCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATaaa                                                                                        >  2:1263367/1‑86 (MQ=255)
                                                                                           aaTCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATaaa                                                                                        <  1:1263367/86‑1 (MQ=255)
                                                                                           aaTCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCtt                                                <  1:1220478/126‑1 (MQ=255)
                                                                                           aaTCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCtt                                                >  2:1220478/1‑126 (MQ=255)
                                                                                           aaTCATCAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCg                                   >  1:673331/1‑139 (MQ=255)
                                                                                            atcatcAATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCgg                                  <  2:866361/139‑1 (MQ=255)
                                                                                             tcatcaATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCggg                                 <  2:1567587/139‑1 (MQ=255)
                                                                                              catcaATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCtt                                                <  1:1063856/123‑1 (MQ=255)
                                                                                              catcaATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCgggg                                >  2:213062/1‑139 (MQ=255)
                                                                                              catcaATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCgggg                                >  2:133884/1‑139 (MQ=255)
                                                                                              catcaATCAGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGACACGACCCGCtt                                                >  2:1063856/1‑123 (MQ=255)
                                                                                                tcaatcaGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAg                              >  1:596455/1‑139 (MQ=255)
                                                                                                tcaatcaGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAg                              <  2:1568251/139‑1 (MQ=255)
                                                                                                   atcaGGTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCAt                           >  2:1296863/1‑139 (MQ=255)
                                                                                                       ggTAAGAGTCTGGAATTTCACACTGCACCCTTTATAc                                                                                                                             <  1:235931/37‑1 (MQ=255)
                                                                                                       ggTAAGAGTCTGGAATTTCACACTGCACCCTTTATAc                                                                                                                             >  2:235931/1‑37 (MQ=255)
                                                                                                        gTAAGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTAccc                      >  2:105482/1‑139 (MQ=255)
                                                                                                          aaGAGTCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAA‑‑ACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGctt                  <  1:1506303/139‑1 (MQ=255)
                                                                                                               tCTGGAATTTCACACTGTACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGcttctt               >  2:1470907/1‑139 (MQ=255)
                                                                                                               tCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCg                     <  2:12768/133‑1 (MQ=255)
                                                                                                               tCTGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCg                     >  1:12768/1‑133 (MQ=255)
                                                                                                                 tGGAATTTCACTCTG‑‑CCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTGCCGGGGAGCATTACCCGCTTCTTTTTa           <  2:1065884/139‑1 (MQ=255)
                                                                                                                 tGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACaa                                                                <  2:741692/88‑1 (MQ=255)
                                                                                                                 tGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACaa                                                                >  1:741692/1‑88 (MQ=255)
                                                                                                                 tGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCtttt             >  1:769001/1‑139 (MQ=255)
                                                                                                                 tGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCtttt             <  1:428289/139‑1 (MQ=255)
                                                                                                                 tGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCtttt             <  2:629956/139‑1 (MQ=255)
                                                                                                                 tGGAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCtttt             >  1:813502/1‑139 (MQ=255)
                                                                                                                  ggAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCggg                                 >  1:778809/1‑118 (MQ=255)
                                                                                                                  ggAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCggg                                 <  2:778809/118‑1 (MQ=255)
                                                                                                                  ggAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCttttt            >  1:1082538/1‑139 (MQ=255)
                                                                                                                  ggAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCttttt            >  2:257192/1‑139 (MQ=255)
                                                                                                                  ggAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCttttt            >  2:1477111/1‑139 (MQ=255)
                                                                                                                   gAATTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTa           <  1:1214124/139‑1 (MQ=255)
                                                                                                                     aTTTCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACt         <  2:664438/139‑1 (MQ=255)
                                                                                                                        tCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGt                                                                             >  2:559420/1‑68 (MQ=255)
                                                                                                                        tCACACTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGt                                                                             <  1:559420/68‑1 (MQ=255)
                                                                                                                         cacaCTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTg     >  1:22300/1‑139 (MQ=255)
                                                                                                                           cacTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGAt   <  1:133884/139‑1 (MQ=255)
                                                                                                                           cacTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGAt   <  2:673331/139‑1 (MQ=255)
                                                                                                                            acTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATc  <  1:213062/139‑1 (MQ=255)
                                                                                                                            acTGCACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATc  <  2:905022/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                      
TGGGCTTTTTGTAGGTTATTTGTACAGCAAAATGGCGCTTGTACATCTATTTCCCCCAATGCAGGATGATAAATATCACGGGAGAATAGAGAATCATCAATCAGGTAAGAGTCTGGAATTTCACACTGTACCCTTTATACTGCCCTATCACTTCGCGAAGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATC  >  NC_000913/2313306‑2313568

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