Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,597,313 (C)7→8 pseudogene (675/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,597,3061.C91.4% 267.4 / 21.1 82pseudogene (682/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base C (32/42);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACC  >  NC_000913/1597208‑1597399
                                                                                                     |                                                                                                
gagTCTGTGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca‑‑a‑‑c                                                                                                <  2:728362/100‑1 (MQ=255)
 agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca‑‑a‑‑cc                                                                                               >  1:2179348/1‑100 (MQ=255)
  gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑ccc                                                                                             >  1:1188107/1‑100 (MQ=255)
  gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑ccc                                                                                             <  2:1728837/100‑1 (MQ=255)
      gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑ccccccc                                                                                        <  1:1394975/100‑1 (MQ=255)
       cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑cccccccc                                                                                       >  2:1714398/1‑99 (MQ=255)
         tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑cccccccctc                                                                                     <  1:995864/100‑4 (MQ=255)
         tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑cccccccctc                                                                                     >  2:159138/1‑97 (MQ=255)
         tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑cccccccctc                                                                                     <  2:1821846/100‑4 (MQ=255)
          aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑cccccccctcg                                                                                    >  2:2035798/1‑96 (MQ=255)
            cGTGATGACAGAGTCTTGATCCTGCACATTCACGGTTCCTATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTcgcg                                                                                   <  1:1565028/100‑1 (MQ=255)
             gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGa                                                                                  <  1:2145055/100‑1 (MQ=255)
             gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGa                                                                                  >  2:2105959/1‑100 (MQ=255)
                   aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTc                                                                            <  2:556614/100‑1 (MQ=255)
                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCg                                                                           >  2:78428/1‑100 (MQ=255)
                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCg                                                                           >  2:87245/1‑100 (MQ=255)
                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCg                                                                           >  2:1356358/1‑100 (MQ=255)
                    cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCg                                                                           <  2:1286792/100‑1 (MQ=255)
                     agagTCTTGATCCTGCACATTCAGGGTTCCGGTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACCC‑‑CCCCCCCTCGCGAATAg                                                                              <  1:2194588/97‑1 (MQ=255)
                        gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTc                                                                       >  1:783241/1‑100 (MQ=255)
                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTccc                                                                     >  1:243404/1‑100 (MQ=255)
                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTccc                                                                     >  2:160921/1‑100 (MQ=255)
                          cTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTccc                                                                     >  2:1871780/1‑100 (MQ=255)
                           ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTcccc                                                                    <  1:1121428/100‑1 (MQ=255)
                                 cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                              >  1:1019880/1‑100 (MQ=255)
                                 cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                              >  1:618963/1‑100 (MQ=255)
                                  tGCACATTCATGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                             <  2:989846/100‑1 (MQ=255)
                                  tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                             >  2:1231406/1‑100 (MQ=255)
                                  tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                             <  2:616062/100‑1 (MQ=255)
                                   gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                            <  2:1339942/100‑1 (MQ=255)
                                   gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                            <  2:1336909/100‑1 (MQ=255)
                                   gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                            <  1:78428/100‑1 (MQ=255)
                                    cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                           <  1:2105947/100‑1 (MQ=255)
                                    cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CCTCCCCCTCGCGAATATTCTCCCCCCCAGCTCCTTg                                                           >  2:2087420/1‑100 (MQ=255)
                                      caTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAt                                                         >  2:1401042/1‑100 (MQ=255)
                                       aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                        <  1:2221595/100‑1 (MQ=255)
                                       aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                        >  2:1560969/1‑100 (MQ=255)
                                         tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                      <  1:2023363/100‑1 (MQ=255)
                                         tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                     <  2:436867/100‑1 (MQ=255)
                                          cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                     <  1:2035786/100‑1 (MQ=255)
                                          cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                     <  2:1391268/100‑1 (MQ=255)
                                          cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                     <  2:1661083/100‑1 (MQ=255)
                                           aGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                                                    >  1:9278/1‑100 (MQ=255)
                                                 ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CCCCTCCCTCGCGAATAGTCGCCTCCCCACTTCCTTAATTACCAAATTGa                                              >  2:245096/1‑100 (MQ=255)
                                                  cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGan                                             <  2:618964/100‑2 (MQ=255)
                                                  cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                             <  2:438772/100‑1 (MQ=255)
                                                     tACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                          <  2:143653/100‑1 (MQ=255)
                                                      aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAttn                                         <  1:679438/100‑2 (MQ=255)
                                                       ccGGTTCCAGTGCCACCGATAATAGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                         <  2:1993859/99‑1 (MQ=255)
                                                       ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                        >  2:1096120/1‑100 (MQ=255)
                                                         ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                      <  1:1356356/100‑1 (MQ=255)
                                                         ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                     >  1:817184/1‑100 (MQ=255)
                                                             gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                  >  2:10987/1‑100 (MQ=255)
                                                              cATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgn                                 <  1:1155020/100‑2 (MQ=255)
                                                              cATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                 <  1:1882912/100‑1 (MQ=255)
                                                              cATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGa                                >  2:1561251/1‑100 (MQ=255)
                                                                ttGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ATCCCCCCCCTCGCGAATAGTCGCCTCTCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                 <  1:687491/99‑1 (MQ=255)
                                                                ttGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                               <  2:1621129/100‑1 (MQ=255)
                                                                  gCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCTCCTCCCCCGTTCCTTGATTACCAATTTGAACTTCAATTGAGGAATc                             >  1:1928830/1‑100 (MQ=255)
                                                                  gCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                             >  1:714819/1‑100 (MQ=255)
                                                                   ccACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CACCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAGTACCAATTTGAAATTCAATTGAGGAATCat                           <  1:2192808/100‑1 (MQ=255)
                                                                    cACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCGCATTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                           >  1:1561197/1‑100 (MQ=255)
                                                                    cACCGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                           <  2:417229/100‑1 (MQ=255)
                                                                       cGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                        >  2:1394551/1‑100 (MQ=255)
                                                                       cGATAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                        <  1:760171/100‑1 (MQ=255)
                                                                          tAATCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                     >  2:1861176/1‑100 (MQ=255)
                                                                           aaTCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                    <  2:2179360/100‑1 (MQ=255)
                                                                           aaTCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                    <  1:1526685/100‑1 (MQ=255)
                                                                            aTCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                   >  1:1971421/1‑100 (MQ=255)
                                                                             tCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                  <  1:1267092/100‑1 (MQ=255)
                                                                             tCGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGCTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                   >  1:363148/1‑99 (MQ=255)
                                                                              cGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCGCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                 <  1:87245/100‑1 (MQ=255)
                                                                              cGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                 <  1:1817339/100‑1 (MQ=255)
                                                                              cGTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                 <  1:1258859/100‑1 (MQ=255)
                                                                               gTATTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                >  1:351570/1‑100 (MQ=255)
                                                                               gTATTTCCAGCCGTAACCA‑‑ACCCCCCCCTCGCAAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                <  1:917960/100‑1 (MQ=255)
                                                                                 aTTTTCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg              <  1:159138/100‑1 (MQ=255)
                                                                                  ttttCAGCCGTAACCA‑‑ACCCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg              <  2:105851/100‑1 (MQ=255)
                                                                                  ttttCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTTTTAGc             >  1:810682/1‑100 (MQ=255)
                                                                                   tttCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc            <  1:733358/100‑1 (MQ=255)
                                                                                   tttCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc            <  2:1904431/100‑1 (MQ=255)
                                                                                   tttCAGCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc            >  2:271032/1‑100 (MQ=255)
                                                                                      cAGCCGTAACCA‑‑ACCCCCCCCTAGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCCATTTGAAATTCAATTGAAGAATCATTATTTTTTATTAGCCACt         >  2:1729370/1‑100 (MQ=255)
                                                                                        gCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg       >  2:403417/1‑100 (MQ=255)
                                                                                        gCCGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg       <  2:1245181/100‑1 (MQ=255)
                                                                                         ccGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTccc      >  1:1202434/1‑98 (MQ=255)
                                                                                         ccGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc      >  2:1578178/1‑100 (MQ=255)
                                                                                          cGTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc     >  1:193035/1‑100 (MQ=255)
                                                                                           gTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGGTTACCAAGTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGGCACTCGCCa    >  2:1889889/1‑100 (MQ=255)
                                                                                           gTAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCGCCCCAGCTCCTTTATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc     <  1:1560961/99‑1 (MQ=255)
                                                                                            tAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc   <  2:1986437/100‑1 (MQ=255)
                                                                                            tAACCA‑‑ACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc   <  2:193035/100‑1 (MQ=255)
                                                                                             aacca‑‑accCCCCCCTCGCGAATAGTCGCCCCCCCAGTTCCTTGATTACCCATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAcc  >  1:926490/1‑100 (MQ=255)
                                                                                                     |                                                                                                
GAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑‑A‑‑CCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACC  >  NC_000913/1597208‑1597399

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