Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,256 (C)8→9 pseudogene (3831/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,2481.C99.5% 802.8 / ‑3.1 191pseudogene (3839/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (126/64);  total (127/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

Only 100 of 220 total aligned reads displayed.

ATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCG  >  NC_000913/1594148‑1594343
                                                                                                     |                                                                                               
atCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAttt                                                                                                  >  1:2334946/1‑101 (MQ=255)
  cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc                                                                                               <  2:2334587/101‑1 (MQ=255)
  cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc                                                                                               >  2:2334947/1‑101 (MQ=255)
     cTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccc                                                                                            >  1:2334952/1‑101 (MQ=255)
      tGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccc                                                                                           >  2:2334953/1‑101 (MQ=255)
         aGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccccc                                                                                        >  2:2334958/1‑100 (MQ=255)
         aGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccccc                                                                                        <  1:2334681/101‑2 (MQ=255)
         aGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccccc                                                                                        >  1:2334959/1‑100 (MQ=255)
          gtggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccccca                                                                                       <  1:2334647/101‑3 (MQ=255)
          gtggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccccca                                                                                       >  2:2334962/1‑99 (MQ=255)
              tgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGat                                                                                    <  1:2334722/101‑1 (MQ=255)
               gACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGata                                                                                   >  1:2334969/1‑101 (MQ=255)
               gACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGata                                                                                   >  1:2334967/1‑101 (MQ=255)
                aCATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGatat                                                                                  >  2:2334971/1‑101 (MQ=255)
                 cATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATc                                                                                 >  1:2334974/1‑101 (MQ=255)
                   tttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcg                                                                               <  2:2334670/101‑1 (MQ=255)
                   tttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccccc                                                                                        >  1:2334976/1‑90 (MQ=255)
                    ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcgt                                                                              >  2:2334977/1‑101 (MQ=255)
                    ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcgt                                                                              <  1:2334716/101‑1 (MQ=255)
                        cGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTa                                                                          >  2:2334981/1‑101 (MQ=255)
                        cGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTa                                                                          <  1:2334586/101‑1 (MQ=255)
                         gACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATcc                                                                                >  1:2334984/1‑94 (MQ=255)
                          aCATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTAcc                                                                        <  1:2334833/101‑1 (MQ=255)
                             tCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTg                                                                     >  2:2334986/1‑101 (MQ=255)
                               acacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGAt                                                                   >  2:2334988/1‑101 (MQ=255)
                                 acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                 >  2:2334992/1‑101 (MQ=255)
                                 acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                 >  1:2334991/1‑101 (MQ=255)
                                  ccGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGAt                                                                   <  2:2334728/98‑1 (MQ=255)
                                   cGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAcgc                                                               <  2:2334677/101‑1 (MQ=255)
                                    gCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAcgcc                                                              >  2:2334996/1‑101 (MQ=255)
                                    gCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAcgcc                                                              >  1:2334995/1‑101 (MQ=255)
                                      tGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGa                                                                    >  2:2334999/1‑93 (MQ=255)
                                       gTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                                           >  1:2335001/1‑101 (MQ=255)
                                          aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                        <  1:2334521/101‑1 (MQ=255)
                                          aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                        >  1:2335004/1‑101 (MQ=255)
                                          aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                        <  1:2334695/101‑1 (MQ=255)
                                           tGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAAt                                                       >  2:2335007/1‑101 (MQ=255)
                                           tGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccc                                                                                           >  2:2335009/1‑64 (MQ=255)
                                            ggATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                      >  1:2335011/1‑101 (MQ=255)
                                             gATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                      <  1:2334750/100‑1 (MQ=255)
                                             gATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                     >  1:2335013/1‑101 (MQ=255)
                                              atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAcg                                                                >  1:2335017/1‑89 (MQ=255)
                                              atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgt                                                    >  1:2335015/1‑101 (MQ=255)
                                               tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                      >  1:2335020/1‑98 (MQ=255)
                                               tatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                   >  2:2335019/1‑101 (MQ=255)
                                                atTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                  >  2:2335023/1‑101 (MQ=255)
                                                  tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                >  1:2335026/1‑101 (MQ=255)
                                                   aaTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                               <  1:2334617/101‑1 (MQ=255)
                                                    aTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTaa                                              >  2:2335028/1‑101 (MQ=255)
                                                       tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAagag                                           >  1:2335029/1‑101 (MQ=255)
                                                         ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTc                                         <  1:2334691/101‑1 (MQ=255)
                                                         ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTc                                         >  1:2335031/1‑101 (MQ=255)
                                                         ttGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCATACCCTGATACGCCGCGGAATCgttgtt                                                <  2:2334638/94‑1 (MQ=255)
                                                          tGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCa                                        >  1:2335033/1‑101 (MQ=255)
                                                          tGACTGGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCa                                        <  1:1511618/101‑1 (MQ=255)
                                                           gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                       >  1:2335037/1‑101 (MQ=255)
                                                            aCTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                      <  1:2334736/101‑1 (MQ=255)
                                                            aCTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                      >  2:2335040/1‑101 (MQ=255)
                                                            aCTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                      <  2:2334712/101‑1 (MQ=255)
                                                              ttGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGAt                                    >  1:2335042/1‑101 (MQ=255)
                                                                 ccATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATccc                                 >  2:2335044/1‑101 (MQ=255)
                                                                 ccATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATccc                                 >  1:2335045/1‑101 (MQ=255)
                                                                  cATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                               >  1:2335049/1‑101 (MQ=255)
                                                                   atatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                               >  2:2335051/1‑101 (MQ=255)
                                                                   atatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                               <  2:2334803/101‑1 (MQ=255)
                                                                   atatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                               >  1:2335052/1‑101 (MQ=255)
                                                                    tatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGt                              >  2:2335054/1‑101 (MQ=255)
                                                                       tAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAcc                           >  1:2335057/1‑101 (MQ=255)
                                                                        aaTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAccc                          >  2:2335060/1‑101 (MQ=255)
                                                                         aTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCa                         >  1:2335062/1‑101 (MQ=255)
                                                                          tGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGAACCCAg                        <  1:2334971/101‑1 (MQ=255)
                                                                           tgCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATc                                                                                 >  2:2335068/2‑43 (MQ=255)
                                                                           ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGt                       >  2:2335066/1‑101 (MQ=255)
                                                                            gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                      <  2:2334690/101‑1 (MQ=255)
                                                                            gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                      <  1:2334734/101‑1 (MQ=255)
                                                                            gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGGAcc                           >  2:2335071/1‑96 (MQ=255)
                                                                              agagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                                           <  2:2334530/62‑1 (MQ=255)
                                                                              agagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtttt                    <  2:2334762/101‑1 (MQ=255)
                                                                              agagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtttt                    >  2:2335073/1‑101 (MQ=255)
                                                                               gagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                   <  2:2334861/101‑1 (MQ=255)
                                                                               gagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                   <  2:2334752/101‑1 (MQ=255)
                                                                               gagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                   >  1:2335074/1‑101 (MQ=255)
                                                                               gagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                   <  1:2334602/101‑1 (MQ=255)
                                                                                 gTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtttt                    >  2:2335076/1‑98 (MQ=255)
                                                                                   cGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac               >  1:2335079/1‑101 (MQ=255)
                                                                                    ggAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca              >  2:2335082/1‑101 (MQ=255)
                                                                                    ggAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca              >  1:2335083/1‑101 (MQ=255)
                                                                                    ggAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca              <  1:2334806/101‑1 (MQ=255)
                                                                                    ggAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca              <  1:2334669/101‑1 (MQ=255)
                                                                                     gAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt             >  2:2335086/1‑101 (MQ=255)
                                                                                      aaCCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATc            >  1:2335087/1‑101 (MQ=255)
                                                                                       aCCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCa           >  2:2335090/1‑101 (MQ=255)
                                                                                         cGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTca                <  2:2334672/94‑1 (MQ=255)
                                                                                         cGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc         >  2:2335094/1‑101 (MQ=255)
                                                                                         cGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc         >  1:2335093/1‑101 (MQ=255)
                                                                                             aGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATAcgccgc                                                            >  2:2335098/1‑46 (MQ=39)
                                                                                             aGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgcta     <  2:2334821/101‑1 (MQ=255)
                                                                                             aGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgcta     >  1:2335097/1‑101 (MQ=255)
                                                                                                aaTTTCCCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                 >  1:2335102/1‑38 (MQ=37)
                                                                                                aaTTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctaccg  <  2:2334625/101‑1 (MQ=255)
                                                                                                     |                                                                                               
ATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCG  >  NC_000913/1594148‑1594343

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