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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,2490C.97.3% 495.3 / 6.5 110pseudogene (3739/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base . (60/47);  total (62/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 122 total aligned reads displayed.

CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCG  >  NC_000913/1594150‑1594343
                                                                                                   |                                                                                              
cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                                >  2:87318/1‑100 (MQ=255)
cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                                >  1:2785539/1‑100 (MQ=255)
cACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                                <  1:2585125/100‑1 (MQ=255)
 aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc                                                                                               <  2:1741765/100‑1 (MQ=255)
 aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcc                                                                                               <  1:1530334/100‑1 (MQ=255)
  ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                                                              >  1:159773/1‑100 (MQ=255)
   cTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccc                                                                                             <  2:1303321/100‑1 (MQ=255)
     gCAGTGGTGACATTTCCCTACATCACACCTCTGTTATAGATATTAAGGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTGcccccc                                                                                           <  1:739275/100‑1 (MQ=255)
          ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccacga                                                                                      >  1:2663846/1‑96 (MQ=255)
           gtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGat                                                                                    >  2:1125161/1‑100 (MQ=255)
            tgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGata                                                                                   <  2:1294755/100‑1 (MQ=255)
                  ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtc                                                                             <  1:194996/100‑1 (MQ=255)
                  ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtc                                                                             >  1:513603/1‑100 (MQ=255)
                  ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtc                                                                             >  2:255965/1‑100 (MQ=255)
                   tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGCACCGAAAGCAATTTCCCCCCCCACGATATCcgtc                                                                             >  1:541159/1‑100 (MQ=255)
                   tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcg                                                                            <  2:2123474/100‑1 (MQ=255)
                   tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcg                                                                            >  2:1433036/1‑100 (MQ=255)
                   tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcg                                                                            >  1:2053401/1‑100 (MQ=255)
                   tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcg                                                                            <  2:1419158/100‑1 (MQ=255)
                   tCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGACATATTAATGGCAAAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcg                                                                            <  1:639144/100‑1 (MQ=255)
                    cccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCcgtcgt                                                                           >  1:2097539/1‑100 (MQ=255)
                     ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTa                                                                          <  2:2391410/100‑1 (MQ=255)
                     ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTa                                                                          >  2:1839706/1‑100 (MQ=255)
                     ccGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTa                                                                          >  2:819292/1‑100 (MQ=255)
                       gACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTAcc                                                                        >  2:816998/1‑100 (MQ=255)
                       gACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTAcc                                                                        >  1:142189/1‑100 (MQ=255)
                         cATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCACTTTCCCCCCCACGATATCCGTCGTACCCt                                                                      <  1:1684378/100‑1 (MQ=255)
                         cATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCt                                                                      >  1:151602/1‑100 (MQ=255)
                          aTCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCt                                                                      >  2:2779624/1‑100 (MQ=255)
                          aTCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTg                                                                     <  1:1125153/100‑1 (MQ=255)
                           tCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGa                                                                    >  1:2156565/1‑100 (MQ=255)
                           tCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGa                                                                    <  2:2102708/100‑1 (MQ=255)
                            cacaCCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGAt                                                                   <  2:136448/100‑1 (MQ=255)
                              cacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                 >  2:2458316/1‑100 (MQ=255)
                              cacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                 <  1:2349608/100‑1 (MQ=255)
                              cacCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAc                                                                 >  1:45164/1‑100 (MQ=255)
                                ccGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAcgc                                                               >  2:2024414/1‑100 (MQ=255)
                                ccGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAcgc                                                               >  1:972997/1‑100 (MQ=255)
                                   cTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCTGAACAGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAcaccgc                                                            >  1:328466/1‑100 (MQ=255)
                                   cTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATAcgccgc                                                            <  1:274564/100‑1 (MQ=255)
                                       tATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                        >  2:45792/1‑100 (MQ=255)
                                        aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGCAAGCACTTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAAt                                                       >  1:267233/1‑100 (MQ=255)
                                        aTGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAAt                                                       >  1:700955/1‑100 (MQ=255)
                                         tGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATCTCCGTCGTACCCTGATACGCCGCGGAATc                                                      <  1:1152989/100‑1 (MQ=255)
                                          ggATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                     >  2:446586/1‑100 (MQ=255)
                                          ggATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                     <  2:1369866/100‑1 (MQ=255)
                                           gATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgt                                                    >  2:2054115/1‑100 (MQ=255)
                                           gATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCg                                                     <  1:2428124/99‑1 (MQ=255)
                                            atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGACTCctt                                                   >  1:638763/1‑100 (MQ=255)
                                            atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                   >  2:2808377/1‑100 (MQ=255)
                                            atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                   <  1:1333255/100‑1 (MQ=255)
                                            atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                   <  1:1417692/100‑1 (MQ=255)
                                            atatTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgtt                                                   <  1:1333093/100‑1 (MQ=255)
                                              atTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgt                                                 >  2:720670/1‑100 (MQ=255)
                                               ttAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttgtt                                                <  2:27886/100‑1 (MQ=255)
                                                tAATGTGTTGACTTGCCATATTAATGTCAGAGGCGGAACCGAAAGCAATTTCCCCCACACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                               <  1:1323057/100‑1 (MQ=255)
                                                tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                               >  2:2540039/1‑100 (MQ=255)
                                                 aaTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTaa                                              >  2:1913920/1‑100 (MQ=255)
                                                  aTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                             <  2:1685593/100‑1 (MQ=255)
                                                       ttGACTTGCCATATTAATGACAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCa                                        <  2:2197794/100‑1 (MQ=255)
                                                        tGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCa                                        >  2:1270023/1‑99 (MQ=255)
                                                         gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                      >  2:1790906/1‑100 (MQ=255)
                                                         gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                      >  2:1322960/1‑100 (MQ=255)
                                                          aCTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTCAGAGTCAGGa                                     >  2:1301735/1‑100 (MQ=255)
                                                          aCTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAgg                                      >  2:2340408/1‑99 (MQ=255)
                                                            ttGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATc                                   <  1:87317/100‑1 (MQ=255)
                                                               ccATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCg                                <  2:2053416/100‑1 (MQ=255)
                                                                cATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                               <  2:2595292/100‑1 (MQ=255)
                                                                cATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCgg                               <  1:1472321/100‑1 (MQ=255)
                                                                  tatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTa                             >  1:1614198/1‑100 (MQ=255)
                                                                   atTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAc                            <  2:2617664/100‑1 (MQ=255)
                                                                      aaTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCa                         <  2:813474/100‑1 (MQ=255)
                                                                        tGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGt                       >  1:723853/1‑100 (MQ=255)
                                                                         ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                      >  1:590823/1‑100 (MQ=255)
                                                                         ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                      <  1:2249045/100‑1 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                     >  1:1238236/1‑100 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                     >  2:449990/1‑100 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                     >  2:2741125/1‑100 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                     <  2:2672165/100‑1 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                     >  2:1153446/1‑100 (MQ=255)
                                                                          gCAGAGTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt                     <  1:1322952/100‑1 (MQ=255)
                                                                            cgatctggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                   <  1:181451/94‑1 (MQ=255)
                                                                            agagTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                   <  1:39162/100‑1 (MQ=255)
                                                                             gagTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGt                  <  1:1755384/100‑1 (MQ=255)
                                                                              agTCGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTc                 <  2:1869693/100‑1 (MQ=255)
                                                                               gTCGGACCCGAAAGCAATGTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTca                <  1:262824/100‑1 (MQ=255)
                                                                                tCGGAACCCAAAGCTATTTCCCCCCCACGATATCCGTCGGACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATACCGGTACCCAGCTTTGTcaa               <  1:33659/100‑2 (MQ=255)
                                                                                 cGGAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca              <  2:1238244/100‑1 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCGCCCTACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcaca              <  1:1400954/99‑1 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTgagatcg             >  2:181453/1‑94 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcacct             >  2:1039354/1‑98 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt             >  2:2455684/1‑100 (MQ=255)
                                                                                  ggAACCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACAt             >  1:1638037/1‑100 (MQ=255)
                                                                                    aaCCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCa           >  1:945793/1‑100 (MQ=255)
                                                                                     aCCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac          >  2:454851/1‑100 (MQ=255)
                                                                                     aCCGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCac          >  2:1161078/1‑100 (MQ=255)
                                                                                      ccGAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc         >  2:1570218/1‑100 (MQ=255)
                                                                                        gAAAGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgc       >  2:1401405/1‑100 (MQ=255)
                                                                                         aaaGCAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCACGATCACGGTACCCAGTTTTGTCACATCaccgct      >  1:200234/1‑100 (MQ=255)
                                                                                             cAATTTCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctaccg  >  1:2821478/1‑100 (MQ=255)
                                                                                                   |                                                                                              
CACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCG  >  NC_000913/1594150‑1594343

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