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.8% 269.2 / NA 84*434* (TAG→TAArlmD23S rRNA m(5)U1939 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (43/40);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCA  >  NC_000913/2913573‑2913828
                                                                                                                              |                                                                                                                                 
tAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTTCCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAAcgcgc                                                                                                                       <  1:97193/139‑1 (MQ=255)
         gACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGt                                                                                                              >  2:374801/1‑139 (MQ=255)
         gACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGt                                                                                                              >  1:1089603/1‑139 (MQ=255)
          aCTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTa                                                                                                             >  2:823141/1‑139 (MQ=255)
             tGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          <  1:760757/139‑1 (MQ=255)
              gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          <  2:415958/138‑1 (MQ=255)
              gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                          >  1:415958/1‑138 (MQ=255)
              gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         >  2:728274/1‑139 (MQ=255)
              gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         >  1:416511/1‑139 (MQ=255)
              gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         <  2:502512/139‑1 (MQ=255)
              gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         <  2:615276/139‑1 (MQ=255)
              gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                         <  2:1173341/139‑1 (MQ=255)
               cGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATc                                                                                                        >  2:996553/1‑139 (MQ=255)
                  tCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAt                                                                                                     <  2:367156/139‑1 (MQ=255)
                  tCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAt                                                                                                     <  2:1089603/139‑1 (MQ=255)
                            cGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTc                                                                                           <  2:1130676/139‑1 (MQ=255)
                              gATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                                                         <  1:352236/139‑1 (MQ=255)
                               aTCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                                                                        <  1:728274/139‑1 (MQ=255)
                                tCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgt                                                                                       <  1:386793/139‑1 (MQ=255)
                                   aaTTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtgtg                                                                                    <  2:166787/139‑1 (MQ=255)
                                    aTTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtgtgt                                                                                   <  2:584718/139‑1 (MQ=255)
                                                 aTTGATATGTGCACTTCTTACCGCCACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAACGTACCAGCGAGTTCAGATGTCCCGTGTGTGGAAACATAGCCa                                                                      >  1:662049/1‑139 (MQ=255)
                                                        tgtgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGc                                                               >  2:273665/1‑139 (MQ=255)
                                                         gtgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                                                         >  1:810052/1‑112 (MQ=255)
                                                         gtgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                                                         <  2:810052/112‑1 (MQ=255)
                                                          tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaa                                                                                                                  <  1:426461/86‑1 (MQ=255)
                                                          tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaa                                                                                                                  >  2:426461/1‑86 (MQ=255)
                                                          tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                             >  1:935023/1‑139 (MQ=255)
                                                          tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                             >  1:506614/1‑139 (MQ=255)
                                                          tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                             >  1:184741/1‑139 (MQ=255)
                                                          tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                             >  1:781215/1‑139 (MQ=255)
                                                           gCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAg                                                            <  1:374801/139‑1 (MQ=255)
                                                            cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  >  2:1129384/1‑116 (MQ=255)
                                                            cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  <  1:1129384/116‑1 (MQ=255)
                                                            cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt                                                           >  1:1020643/1‑139 (MQ=255)
                                                                tcttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       <  1:273665/139‑1 (MQ=255)
                                                                  ttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     <  2:233336/139‑1 (MQ=255)
                                                                    aCCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaaaa                                                                                                                <  2:181557/78‑1 (MQ=255)
                                                                    aCCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaaaa                                                                                                                >  1:181557/1‑78 (MQ=255)
                                                                    aCCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAAt                                                   >  2:763612/1‑139 (MQ=255)
                                                                         aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAg                                                            >  1:1204829/1‑125 (MQ=255)
                                                                         aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAg                                                            <  2:1204829/125‑1 (MQ=255)
                                                                         aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATg                                                  >  2:579148/1‑135 (MQ=255)
                                                                         aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATg                                                  <  1:579148/135‑1 (MQ=255)
                                                                         aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGtat                                              >  2:325215/1‑139 (MQ=255)
                                                                                      tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaa                                                                                                                  <  1:850375/58‑1 (MQ=255)
                                                                                      tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaa                                                                                                                  >  2:850375/1‑58 (MQ=255)
                                                                                      tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcg                                                         >  1:101475/1‑115 (MQ=255)
                                                                                      tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcg                                                         <  2:101475/115‑1 (MQ=255)
                                                                                           tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            >  2:1126591/1‑139 (MQ=255)
                                                                                           tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            <  1:996553/139‑1 (MQ=255)
                                                                                           tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            >  1:211323/1‑139 (MQ=255)
                                                                                           tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            <  1:763612/139‑1 (MQ=255)
                                                                                             gggACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATCTCCCGTCTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGCTATATCCTGCTTTTAATAACGCt                          >  2:891858/1‑139 (MQ=255)
                                                                                              ggACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCtt                         >  2:1173380/1‑139 (MQ=255)
                                                                                                aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGt                                                                                            <  1:514484/70‑1 (MQ=255)
                                                                                                aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGt                                                                                            >  2:514484/1‑70 (MQ=255)
                                                                                                aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            <  2:1002512/134‑1 (MQ=255)
                                                                                                aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                            >  1:1002512/1‑134 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTaa                                                                                                                            <  1:411616/36‑1 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTaa                                                                                                                            >  2:411616/1‑36 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaa                                                                                                                  <  1:371270/46‑1 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaa                                                                                                                  >  2:371270/1‑46 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAt                                                                                                     <  2:290482/59‑1 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAt                                                                                                     >  1:290482/1‑59 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                                                                             >  1:894787/1‑83 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                                                                             <  2:894787/83‑1 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       <  2:634639/105‑1 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                       >  1:634639/1‑105 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     <  1:1117455/107‑1 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                     >  2:1117455/1‑107 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCt                     >  1:4812/1‑139 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCt                     >  1:1161961/1‑139 (MQ=255)
                                                                                                  cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCt                     >  2:963365/1‑139 (MQ=255)
                                                                                                    gACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTAt                   >  2:195078/1‑139 (MQ=255)
                                                                                                     aCCTGCCGAAATCGGCAAATCGCAATTATTTAAcgcgc                                                                                                                       <  1:354822/38‑1 (MQ=255)
                                                                                                     aCCTGCCGAAATCGGCAAATCGCAATTATTTAAcgcgc                                                                                                                       >  2:354822/1‑38 (MQ=255)
                                                                                                          ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  <  2:948044/70‑1 (MQ=255)
                                                                                                          ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                                                                  >  1:948044/1‑70 (MQ=255)
                                                                                                            gAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                                                                      <  1:29569/80‑1 (MQ=255)
                                                                                                            gAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                                                                      >  2:29569/1‑80 (MQ=255)
                                                                                                              aaTCGTCAAATCGCAATGATTTAACCCGCGACAAAAGTACCACCGATTCCATAGGTCCCGTGTGTCGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAg         <  1:150480/139‑1 (MQ=255)
                                                                                                                     aaaTCGCAATTATTTAACGCGCGATAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACAGACCCATCATATGTAGTCGCGCAACGGTATATGCTGCTTGTACTAACGCTGCGCTATCCCGAGCCAGCGTTGCa  <  1:808548/139‑1 (MQ=255)
                                                                                                                     aaaTCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCa  <  1:325215/139‑1 (MQ=255)
                                                                                                                     aaaTCGCAACTATTTAACGCGCGAGAAAAGTACCCTCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCa  <  2:781215/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                 
TAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCA  >  NC_000913/2913573‑2913828

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