Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,090,046 G→A *17* (TAG→TAA hisL → his operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,0460GA97.6% 249.9 / ‑2.2 82*17* (TAG→TAAhisLhis operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (40/40);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGA  >  NC_000913/2089915‑2090180
                                                                                                                                   |                                                                                                                                      
tataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCa                                                                                                                                 >  2:46228/1‑139 (MQ=255)
     aaaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGa                                                                                                                            >  1:395623/1‑139 (MQ=255)
       aaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtg                                                                                                                          >  2:787542/1‑139 (MQ=255)
       aaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtg                                                                                                                          <  1:46228/139‑1 (MQ=255)
        aGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgt                                                                                                                         >  2:220927/1‑139 (MQ=255)
        aGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgt                                                                                                                         >  1:626896/1‑139 (MQ=255)
        aGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTATATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgt                                                                                                                         >  2:137742/1‑139 (MQ=255)
          cccTTGCTTTCTAACGTGACAGCGGCTTAGGTTAAAAGACATCAGTTGAATAAACATTCACATAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACGAATCTTTCAGGCGAtgtgt                                                                                                                       >  2:846463/1‑139 (MQ=255)
          cccTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtgt                                                                                                                       >  2:1151596/1‑139 (MQ=255)
          cccTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtgt                                                                                                                       >  1:145820/1‑139 (MQ=255)
             ttGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCt                                                                                                                                     <  1:1094541/122‑1 (MQ=255)
             ttGCTTTCTAACGTGAAAGCGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCt                                                                                                                                     >  2:1094541/1‑122 (MQ=255)
              tGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                                >  2:64493/1‑126 (MQ=255)
              tGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                                <  1:64493/126‑1 (MQ=255)
                 tttCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                                <  1:220927/139‑1 (MQ=255)
                 tttCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                                <  2:316458/139‑1 (MQ=255)
                        cGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTc                                                                                                         <  2:1199036/139‑1 (MQ=255)
                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                              <  1:386308/139‑1 (MQ=255)
                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                              <  1:787542/139‑1 (MQ=255)
                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                              <  2:367998/139‑1 (MQ=255)
                                    ttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAgt                                                                                             <  2:395623/139‑1 (MQ=255)
                                    ttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAgt                                                                                             <  2:1264402/139‑1 (MQ=255)
                                       ggTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAgtggt                                                                                          <  2:626896/139‑1 (MQ=255)
                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATg                                                                                     >  2:614776/1‑139 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTc                                                                                                 <  2:1136919/121‑1 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTc                                                                                                 >  1:1136919/1‑121 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCa                                                                               >  2:1107087/1‑139 (MQ=255)
                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCa                                                                               >  2:78563/1‑139 (MQ=255)
                                                              aaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                   <  2:229512/139‑1 (MQ=255)
                                                               aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                           >  2:731493/1‑82 (MQ=255)
                                                               aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                           <  1:731493/82‑1 (MQ=255)
                                                               aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                   >  1:1016150/1‑138 (MQ=255)
                                                               aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                   <  2:1016150/138‑1 (MQ=255)
                                                                 aTTCATAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                              <  1:781311/109‑1 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCgg                                                                >  2:286142/1‑139 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                        >  1:152802/1‑83 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                        <  2:152802/83‑1 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                              >  2:781311/1‑109 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCgg                                                                >  1:375684/1‑139 (MQ=255)
                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGACACCgg                                                                >  2:94222/1‑139 (MQ=255)
                                                                    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                              >  2:40582/1‑106 (MQ=255)
                                                                    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                              <  1:40582/106‑1 (MQ=255)
                                                                      caGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGAt                                                           <  1:614776/139‑1 (MQ=255)
                                                                           aCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                                <  2:626052/65‑1 (MQ=255)
                                                                           aCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                                >  1:626052/1‑65 (MQ=255)
                                                                                 ataatacGCGTCCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGATGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                >  2:751387/6‑139 (MQ=255)
                                                                                     cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAgg                                                                                                                               <  2:388431/56‑1 (MQ=255)
                                                                                     cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAgg                                                                                                                               >  1:388431/1‑56 (MQ=255)
                                                                                     cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCt                                            >  1:1149615/1‑139 (MQ=255)
                                                                                         cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttttt                                        <  1:286142/139‑1 (MQ=255)
                                                                                                tttAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                 <  2:375684/139‑1 (MQ=255)
                                                                                                        ccaccaTCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACAtt                                                                                                          <  2:1030415/58‑1 (MQ=255)
                                                                                                        ccaccaTCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACAtt                                                                                                          >  1:1030415/1‑58 (MQ=255)
                                                                                                         caccatcatcaccatcatCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                   >  1:779384/1‑96 (MQ=255)
                                                                                                         caccatcatcaccatcatCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                   <  2:779384/96‑1 (MQ=255)
                                                                                                         caccatcatcaccatcatCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                        <  2:1149615/139‑1 (MQ=255)
                                                                                                         caccatcatcaccatcatCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                        <  1:1110364/139‑1 (MQ=255)
                                                                                                          accatcatcaccatcatcCTGACTAATCTTTCAgg                                                                                                                               >  1:158898/1‑35 (MQ=255)
                                                                                                          accatcatcaccatcatcCTGACTAATCTTTCAgg                                                                                                                               <  2:158898/35‑1 (MQ=255)
                                                                                                              tcatcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGAt                                                           <  1:1250747/99‑1 (MQ=255)
                                                                                                              tcatcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGAt                                                           >  2:1250747/1‑99 (MQ=255)
                                                                                                                 tcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                        >  2:284083/1‑131 (MQ=255)
                                                                                                                 tcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                        <  1:284083/131‑1 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa              <  1:78563/139‑1 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa              <  2:145820/139‑1 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa              <  1:1151596/139‑1 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATAAACGCATGAGAAAGCCCCCGGAAGAt                                                           >  2:806669/1‑94 (MQ=255)
                                                                                                                   accatcatcCTGACTAATCTTTCAGGAGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATAAACGCATGAGAAAGCCCCCGGAAGAt                                                           <  1:806669/94‑1 (MQ=255)
                                                                                                                    ccatgatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGACCAACTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGAc             <  2:1000126/139‑1 (MQ=255)
                                                                                                                    ccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGAc             <  2:230363/139‑1 (MQ=255)
                                                                                                                     catcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                                   >  1:1218072/1‑68 (MQ=255)
                                                                                                                     catcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                                   <  2:1218072/68‑1 (MQ=255)
                                                                                                                     catcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACa            >  1:48300/1‑139 (MQ=255)
                                                                                                                        catcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                              <  2:667412/102‑1 (MQ=255)
                                                                                                                        catcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                              >  1:667412/1‑102 (MQ=255)
                                                                                                                          tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGTTGCATGAACGCATGAGAAAGCCCCCGGAAGAt                                                           >  1:1243775/1‑87 (MQ=255)
                                                                                                                          tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGAt                                                           <  2:1243775/87‑1 (MQ=255)
                                                                                                                           ccTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                          >  1:1275639/1‑71 (MQ=255)
                                                                                                                           ccTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                          <  2:1275639/71‑1 (MQ=255)
                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa                       <  2:337309/118‑1 (MQ=255)
                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa                       >  1:337309/1‑118 (MQ=255)
                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga  >  1:239552/1‑139 (MQ=255)
                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga  >  2:1246985/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGA  >  NC_000913/2089915‑2090180

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: