Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,913,699 C→T *434* (TAG→TAA rlmD ← 23S rRNA m(5)U1939 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,913,6990CT98.9% 274.1 / NA 87*434* (TAG→TAArlmD23S rRNA m(5)U1939 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (40/46);  minor base G (0/1);  total (40/47)
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.99e-01

GCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCA  >  NC_000913/2913568‑2913828
                                                                                                                                   |                                                                                                                                 
gctggTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTaa                                                                                                                            <  2:1356889/139‑1 (MQ=255)
     tAATCCCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAAcgcgc                                                                                                                       <  2:1484987/139‑1 (MQ=255)
          cccAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaa                                                                                                                  <  2:846734/139‑1 (MQ=255)
              gACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGt                                                                                                              >  1:845528/1‑139 (MQ=255)
               aCTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTaa                                                                                                                            >  2:211173/1‑124 (MQ=255)
               aCTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTaa                                                                                                                            <  1:211173/124‑1 (MQ=255)
               aCTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTa                                                                                                             >  1:1478016/1‑139 (MQ=255)
                  tGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          <  1:1218857/139‑1 (MQ=255)
                   gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         <  1:758218/139‑1 (MQ=255)
                    cGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATc                                                                                                        <  2:1478016/139‑1 (MQ=255)
                    cGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATc                                                                                                        >  2:1499972/1‑139 (MQ=255)
                                 cGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          <  2:861146/124‑1 (MQ=255)
                                 cGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          >  1:861146/1‑124 (MQ=255)
                                 cGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTc                                                                                           <  1:1499972/139‑1 (MQ=255)
                                                 gCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACat                                                                           >  1:949225/1‑139 (MQ=255)
                                                  ccTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACata                                                                          >  1:826039/1‑139 (MQ=255)
                                                      aTTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                                                                      <  2:1525822/139‑1 (MQ=255)
                                                            atGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAAcgcgc                                                                                                                       <  1:1362857/84‑1 (MQ=255)
                                                            atGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAAcgcgc                                                                                                                       >  2:1362857/1‑84 (MQ=255)
                                                            atGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTcc                                                                                          <  2:1022879/113‑1 (MQ=255)
                                                            atGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTcc                                                                                          >  1:1022879/1‑113 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTGCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                             >  2:1040225/1‑139 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTc                                                                                           <  2:243735/109‑1 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTc                                                                                           >  1:243735/1‑109 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  <  2:27477/118‑1 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  <  2:1384974/118‑1 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  >  1:1384974/1‑118 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  >  1:27477/1‑118 (MQ=255)
                                                               tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                             >  2:86844/1‑139 (MQ=255)
                                                                gCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAg                                                            <  1:493468/139‑1 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCga                                                                                                                     >  2:258751/1‑81 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCga                                                                                                                     <  1:258751/81‑1 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGa                                                                                               >  1:1121713/1‑103 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGa                                                                                               <  2:1121713/103‑1 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  >  2:764011/1‑116 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  <  1:764011/116‑1 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTTTGGGACCATATCCAGCATCGCCAGt                                                           >  2:79604/1‑139 (MQ=255)
                                                                 cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt                                                           >  2:44566/1‑139 (MQ=255)
                                                                  aCTTCTTACCGCAACCATCGGCCTCGTCTTTAGGGACCAGACCTGCCGAAATCGTCAAAACGCAAGTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTc                                                          <  1:1476299/139‑1 (MQ=255)
                                                                     tcttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       <  2:617242/139‑1 (MQ=255)
                                                                      cttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTAcc                                                                                                           >  2:156552/1‑86 (MQ=255)
                                                                      cttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTAcc                                                                                                           <  1:156552/86‑1 (MQ=255)
                                                                      cttACCGCAACCATCGTCCTCTCCTTTAGCGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCgcgc                                                      <  1:1040225/139‑1 (MQ=255)
                                                                       ttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     <  2:949225/139‑1 (MQ=255)
                                                                       ttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     <  2:826039/139‑1 (MQ=255)
                                                                              aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTggg                                                                                <  2:626530/105‑1 (MQ=255)
                                                                              aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTggg                                                                                >  1:626530/1‑105 (MQ=255)
                                                                              aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                                                                      >  2:1093311/1‑115 (MQ=255)
                                                                              aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                                                                      <  1:1093311/115‑1 (MQ=255)
                                                                                ccATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTatat                                            >  2:365860/1‑139 (MQ=255)
                                                                                      tCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCt                                      <  1:86844/139‑1 (MQ=255)
                                                                                           tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          <  2:118145/66‑1 (MQ=255)
                                                                                           tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          >  1:118145/1‑66 (MQ=255)
                                                                                              tttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                                                                             >  1:454097/1‑92 (MQ=255)
                                                                                              tttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                                                                             <  2:454097/92‑1 (MQ=255)
                                                                                              tttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTtaataa                              <  2:845528/139‑1 (MQ=255)
                                                                                               ttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAAc                             >  1:773377/1‑139 (MQ=255)
                                                                                                tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                                                                 <  1:809206/70‑1 (MQ=255)
                                                                                                tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                                                                 >  2:809206/1‑70 (MQ=255)
                                                                                                tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                                                                        >  2:779248/1‑79 (MQ=255)
                                                                                                tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                                                                        <  1:779248/79‑1 (MQ=255)
                                                                                                tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            <  2:773377/139‑1 (MQ=255)
                                                                                                   ggACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     >  2:970038/1‑111 (MQ=255)
                                                                                                   ggACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     <  1:970038/111‑1 (MQ=255)
                                                                                                     aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         >  1:561887/1‑57 (MQ=255)
                                                                                                     aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         <  2:561887/57‑1 (MQ=255)
                                                                                                     aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCg                       >  1:320630/1‑139 (MQ=255)
                                                                                                       cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAt                                                                                                     <  1:1344825/59‑1 (MQ=255)
                                                                                                       cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAt                                                                                                     >  2:1344825/1‑59 (MQ=255)
                                                                                                       cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGa                                                                                               >  2:1371611/1‑65 (MQ=255)
                                                                                                       cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGa                                                                                               <  1:1371611/65‑1 (MQ=255)
                                                                                                              gCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       >  1:70633/1‑98 (MQ=255)
                                                                                                              gCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       <  2:70633/98‑1 (MQ=255)
                                                                                                               ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATc                                                                                                        <  2:1356413/48‑1 (MQ=255)
                                                                                                               ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATc                                                                                                        >  1:1356413/1‑48 (MQ=255)
                                                                                                               ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       <  2:348153/97‑1 (MQ=255)
                                                                                                               ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       >  1:348153/1‑97 (MQ=255)
                                                                                                               ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            >  2:405108/1‑124 (MQ=255)
                                                                                                               ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            <  1:405108/124‑1 (MQ=255)
                                                                                                                    aTCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       >  2:1020817/1‑92 (MQ=255)
                                                                                                                    aTCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       <  1:1020817/92‑1 (MQ=255)
                                                                                                                    aTCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     <  2:821718/94‑1 (MQ=255)
                                                                                                                    aTCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     >  1:821718/1‑94 (MQ=255)
                                                                                                                         cAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCCATTCCAGCTGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACCCTTCGCTATCCCGAGCCAGCGTTGc   <  1:79604/139‑1 (MQ=255)
                                                                                                                          aaaTCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCa  <  1:628558/139‑1 (MQ=255)
                                                                                                                             tCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAg                                                            <  2:977935/78‑1 (MQ=255)
                                                                                                                             tCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAGTCCAGATGTCGCGTGTGTGGGAACATGTCCAGCATCGCCAg                                                            >  1:977935/1‑78 (MQ=255)
                                                                                                                                   |                                                                                                                                 
GCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCA  >  NC_000913/2913568‑2913828

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