Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,573,188 G→A 17.4% pseudogene (775/1503 nt) yjiT → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,573,1880GA17.4% 277.5 / 45.6 121pseudogene (775/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base G (46/53);  new base A (12/9);  total (58/63)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 122 total aligned reads displayed.

CTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTT  >  NC_000913/4573089‑4573286
                                                                                                   |                                                                                                  
cTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTg                                                                                                    <  1:356947/100‑1 (MQ=255)
 ttGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                   <  1:2095727/100‑1 (MQ=255)
  tGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTa                                                                                                  <  1:2384821/100‑1 (MQ=255)
     gaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                                                               >  2:2192416/1‑100 (MQ=255)
       aGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGc                                                                                             <  1:987804/100‑1 (MQ=255)
        ggTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCt                                                                                            >  2:424743/1‑100 (MQ=255)
        ggTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCt                                                                                            >  1:2403130/1‑100 (MQ=255)
          ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                          >  1:960691/1‑100 (MQ=255)
          ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                          <  1:743156/100‑1 (MQ=255)
          ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgc                                                                                          <  1:1955463/100‑3 (MQ=255)
          tgcACCCTACCTGACGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                          >  1:1727459/3‑100 (MQ=255)
           tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                         <  2:769631/100‑1 (MQ=255)
           tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2231539/100‑4 (MQ=255)
           tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2316714/100‑4 (MQ=255)
           tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2317231/100‑4 (MQ=255)
           tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:2214680/100‑4 (MQ=255)
           tCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcg                                                                                         <  1:602995/100‑4 (MQ=255)
            cACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                        >  1:1841131/1‑100 (MQ=255)
            cACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                        <  2:2540445/100‑1 (MQ=255)
             aCCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGcc                                                                                       >  1:2489445/1‑100 (MQ=255)
              cccTACCTTGCGACGCGAGATTCAAAAACCGGTAGAATATCAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGACTAGGcc                                                                                       >  2:456915/1‑99 (MQ=255)
               ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATGTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgcata                                                                                     <  1:1474808/100‑8 (MQ=255)
                cTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                    <  1:1508372/100‑1 (MQ=255)
                cTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                    >  2:1835159/1‑100 (MQ=255)
                 tACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgcatagc                                                                                   <  1:673469/100‑10 (MQ=255)
                   ccTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACt                                                                                 >  2:1191932/1‑100 (MQ=255)
                     tGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTAc                                                                               <  2:1177803/100‑1 (MQ=255)
                      ggCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACa                                                                              >  2:1075229/1‑100 (MQ=255)
                          aaGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaa                                                                          <  2:2275467/100‑1 (MQ=255)
                             cGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAAt                                                                       >  2:2348631/1‑100 (MQ=255)
                                gTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAag                                                                    <  1:904664/100‑1 (MQ=255)
                                  ttCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATCagag                                                                  >  2:238160/1‑100 (MQ=255)
                                    ccAATACCGTTAGATGATGAAACAGGAACGCATTTTTTTAATGGACTGTTATGCGCAGCCACTGTAGATGCTAAGGCTAGACTACAAAAAAATAAGAGt                                                                 >  2:555069/1‑99 (MQ=255)
                                     cAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTAGACTACAAAAAAATAAGAGTAc                                                               <  2:1489424/100‑1 (MQ=255)
                                      aaTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcg                                                              >  1:109175/1‑100 (MQ=255)
                                       aTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgc                                                             >  2:1834482/1‑100 (MQ=255)
                                           cGTTAGATGATGAAACAGGAACTCATGTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                         <  1:853260/100‑1 (MQ=255)
                                           cGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                         >  2:1561641/1‑100 (MQ=255)
                                            gTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                        >  1:1381395/1‑100 (MQ=255)
                                              tAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCa                                                      >  1:1133274/1‑100 (MQ=255)
                                               aGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                     >  2:2185352/1‑100 (MQ=255)
                                                gatgatGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAt                                                    >  2:2043126/1‑100 (MQ=255)
                                                 atgatgAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                   >  2:1426653/1‑100 (MQ=255)
                                                  tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  >  1:1933007/1‑100 (MQ=255)
                                                  tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  <  1:1165942/100‑1 (MQ=255)
                                                  tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  >  1:1254679/1‑100 (MQ=255)
                                                  tgatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                  <  1:324349/100‑1 (MQ=255)
                                                   gatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCt                                                 <  1:552399/100‑1 (MQ=255)
                                                      gaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                              <  1:592886/100‑1 (MQ=255)
                                                      gaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                              <  2:769639/100‑1 (MQ=255)
                                                       aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg                                             >  2:386063/1‑100 (MQ=255)
                                                       aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg                                             <  1:2329601/100‑1 (MQ=255)
                                                       aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg                                             <  1:2211135/100‑1 (MQ=255)
                                                        aaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTgg                                            >  1:626659/1‑100 (MQ=255)
                                                         aCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGt                                           >  2:1790264/1‑100 (MQ=255)
                                                          cAGTAACGCATTTCATGAATTGACTGCTATACGCTGCCTCTGTAGAGGCTAAGCATCAACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGtt                                          <  2:129968/100‑2 (MQ=255)
                                                          cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                          <  1:946334/100‑1 (MQ=255)
                                                          cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                          >  1:2034501/1‑100 (MQ=255)
                                                          cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                          <  2:110699/100‑1 (MQ=255)
                                                           aGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCa                                         <  1:283427/100‑1 (MQ=255)
                                                            ggAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAg                                        <  1:1881014/100‑1 (MQ=255)
                                                              aaCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaa                                      >  1:95732/1‑100 (MQ=255)
                                                               aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaa                                     >  2:840121/1‑100 (MQ=255)
                                                               aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaa                                     <  1:134161/100‑1 (MQ=255)
                                                                cGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                    <  2:727451/100‑1 (MQ=255)
                                                                  cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGc                                  >  1:1585397/1‑100 (MQ=255)
                                                                    ttttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAt                                >  1:800168/1‑100 (MQ=255)
                                                                      ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAACAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATcc                              >  2:1226465/1‑100 (MQ=255)
                                                                      ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCCGAAAAGCATcc                              <  2:1967115/100‑1 (MQ=255)
                                                                       tttGAATGGACTGCTATGCGCAGCCTCTTTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCa                             <  2:324874/100‑1 (MQ=255)
                                                                       tttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCa                             >  1:1444316/1‑100 (MQ=255)
                                                                        ttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAg                            >  2:2025751/1‑100 (MQ=255)
                                                                          gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                          >  1:571591/1‑100 (MQ=255)
                                                                          gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                          >  1:1233303/1‑100 (MQ=255)
                                                                          gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGCCAGAAAAGCATCCAGAt                          >  1:2005827/1‑100 (MQ=255)
                                                                           aaTGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGTGTACGCGTTGCCAATTCTATTGGTCAGACAAGCATCCAGacg                         <  1:1737074/100‑3 (MQ=255)
                                                                             tGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAg                       <  1:1043333/100‑1 (MQ=255)
                                                                                 cTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTa                   <  2:1257616/100‑1 (MQ=255)
                                                                                  gtcgaggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAg                  >  2:306769/11‑100 (MQ=38)
                                                                                    cTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAggg                <  2:437852/100‑1 (MQ=255)
                                                                                     tATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt               >  1:2303735/1‑100 (MQ=255)
                                                                                      aggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt               >  1:1451532/7‑99 (MQ=255)
                                                                                      aggcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTg              >  1:2358617/7‑100 (MQ=255)
                                                                                      aTGCGCAGCCTCTGAAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt               <  1:1015553/99‑1 (MQ=255)
                                                                                       tGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTg              >  1:1122795/1‑99 (MQ=255)
                                                                                       tGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa             <  1:2134458/100‑1 (MQ=255)
                                                                                       tGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa             >  1:373783/1‑100 (MQ=255)
                                                                                        gcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGAt            >  2:1954447/5‑100 (MQ=255)
                                                                                         cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATCAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATCAGCTTAGGGTGATa           >  2:16231/4‑100 (MQ=255)
                                                                                         cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  1:1616619/4‑100 (MQ=255)
                                                                                         cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  1:1757760/4‑100 (MQ=255)
                                                                                         cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  1:1409087/4‑100 (MQ=255)
                                                                                         cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           >  2:5043/4‑100 (MQ=255)
                                                                                         cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAATACGCGTTGCCAATTATATTGATCAGAAAAGCATCCAAATAAGCTTAggaaaata           >  2:12718/4‑94 (MQ=255)
                                                                                          gcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATa           <  2:2034487/99‑1 (MQ=255)
                                                                                          gcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAg          <  1:1526316/100‑1 (MQ=255)
                                                                                            tctctctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTa        <  1:1562020/97‑1 (MQ=255)
                                                                                               ctctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTaaga     >  2:1562009/1‑97 (MQ=255)
                                                                                                 ctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATctct   <  2:1333093/100‑1 (MQ=255)
                                                                                                  tGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCtt  >  1:2382758/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTT  >  NC_000913/4573089‑4573286

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: ‑