Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,143,125 C→T 100% D415N (GAC→AAC)  rne ← endoribonuclease; RNA‑binding protein;RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,143,1250CT98.9% 272.5 / ‑5.5 85D415N (GAC→AAC) rneendoribonuclease; RNA‑binding protein;RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (27/56);  total (27/58)
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.97e-01

ATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCC  >  NC_000913/1143028‑1143221
                                                                                                 |                                                                                                
aTTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtac                                                                                                <  2:919793/100‑1 (MQ=255)
 ttGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacg                                                                                               <  1:3040442/100‑1 (MQ=255)
 ttGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacg                                                                                               <  2:2773734/100‑1 (MQ=255)
   ggCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgca                                                                                             >  2:791000/1‑100 (MQ=255)
   ggCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgca                                                                                             <  2:2450294/100‑1 (MQ=255)
   ggCACAGGAACAATGGCGTGAACTTACTGGGTGTTCTCTTTCAGCGCTTCTTCTTCCATCAGACGAAGAATAGAGAGCGACAGCGATTCGTTGTtacgca                                                                                             >  2:1040272/1‑100 (MQ=255)
   gccacaGGCACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCCGACGCAGAATAGCGAGCGCCAGCGATTCGTTGTtacgca                                                                                             >  2:2264771/3‑100 (MQ=255)
       caGGAACAATGGCGTGAACTTCGTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTtacgcacg                                                                                           >  2:539456/1‑98 (MQ=255)
       caGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGt                                                                                         <  1:2919837/100‑1 (MQ=255)
        aGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTg                                                                                        <  1:2919841/100‑1 (MQ=255)
           aaCAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCa                                                                                     >  2:1366365/1‑100 (MQ=255)
            aCAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTTTTCTTCGTTCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAg                                                                                    <  1:988871/100‑1 (MQ=255)
            aCAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAg                                                                                    <  2:1330665/100‑1 (MQ=255)
            aCAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCATCTTCGATCAGACGCAGAAGAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAg                                                                                    >  2:1647120/1‑100 (MQ=255)
             cAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGt                                                                                   <  1:1594327/100‑1 (MQ=255)
              aaTGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTa                                                                                  <  1:2962750/100‑1 (MQ=255)
               aTGGCGTGAACTTCCTGTGTGTTCTCTTTTAGCGCTTCTTCTTAGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTa                                                                                  <  2:1738032/99‑1 (MQ=255)
                tGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTAcc                                                                                <  2:2205717/100‑1 (MQ=255)
                 ggCGTTAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCGTCGATCAGTCGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCa                                                                               <  1:1763940/100‑1 (MQ=255)
                  gCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCa                                                                               <  2:1237128/99‑1 (MQ=255)
                  gCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAg                                                                              >  2:2162217/1‑100 (MQ=255)
                   cGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGa                                                                             <  2:642775/100‑1 (MQ=255)
                        aCTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaacaac                                                                        <  2:1174823/100‑1 (MQ=255)
                        aCTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGaacaac                                                                        >  1:2512356/1‑100 (MQ=255)
                          ttCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAAcgc                                                                      >  1:727165/1‑100 (MQ=255)
                          ttCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAAcgc                                                                      <  1:2710751/100‑1 (MQ=255)
                           tCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAAcgcg                                                                     <  1:352591/100‑1 (MQ=255)
                           tCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAAcgcg                                                                     >  2:2961310/1‑100 (MQ=255)
                            ccTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCgg                                                                    <  1:1082423/100‑1 (MQ=255)
                             cTGGGTGTTCTCTTTTTGCGCCTCTTCTTCGATCAGACTCAGAATAGTTATCGACAGCGATTCGTTGTTACGCTCGGTGCCTGTACCACAACAACTCGGa                                                                   <  2:689322/100‑1 (MQ=255)
                             cTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGa                                                                   <  2:1958796/100‑1 (MQ=255)
                              tGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGTTTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGAc                                                                  <  2:1914388/100‑1 (MQ=255)
                              tGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGa                                                                   <  2:1448258/99‑1 (MQ=255)
                              tGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGAc                                                                  <  2:1785855/100‑1 (MQ=255)
                              tGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGAc                                                                  <  2:544665/100‑1 (MQ=255)
                              tGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGAc                                                                  <  1:2648671/100‑1 (MQ=255)
                               gggTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACa                                                                 <  2:2948456/100‑1 (MQ=255)
                               gggTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACa                                                                 >  2:216050/1‑100 (MQ=255)
                                ggTGTTCTCTTTCAGCGCTTCTTCTTCTTTCATACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCATTACCAGAACAACGCGGACaa                                                                <  2:1845046/100‑1 (MQ=255)
                                   gtTCTCTTTCAGCGCTTCTTCTTCTATCTTTTGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACg                                                             <  1:636795/100‑1 (MQ=255)
                                   gtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACg                                                             <  1:1697916/100‑1 (MQ=255)
                                   gtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACg                                                             >  1:1115194/1‑100 (MQ=255)
                                   gtTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACg                                                             >  1:1009989/1‑100 (MQ=255)
                                    ttCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACg                                                             <  2:2997094/99‑1 (MQ=255)
                                    ttCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGt                                                            >  1:2481400/1‑100 (MQ=255)
                                     tctcTTTCAGCGCTTCTTCTTCGATTAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtg                                                           <  1:1187826/100‑1 (MQ=255)
                                     tctcTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtg                                                           <  2:782185/100‑1 (MQ=255)
                                     tctcTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtg                                                           <  2:1792644/100‑1 (MQ=255)
                                       tctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtgat                                                         <  1:2959516/100‑1 (MQ=255)
                                       tctTTCAGCGCTTCTTCTTCGATCAGACGCAGAATACAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGtg                                                           <  2:1646485/98‑1 (MQ=255)
                                          ttCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGAc                                                      >  2:79413/1‑100 (MQ=255)
                                          ttCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGAc                                                      <  1:760811/100‑1 (MQ=255)
                                            cAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGCGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTg                                                    >  2:41572/1‑100 (MQ=255)
                                            cAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTg                                                    >  1:1752758/1‑100 (MQ=255)
                                               cgcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGAt                                                 >  1:1704217/1‑100 (MQ=255)
                                                gcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGAtt                                                <  1:1227251/100‑1 (MQ=255)
                                                gcTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAAAGCGGACAAACGTGATGACTGGAtt                                                >  1:1244854/1‑100 (MQ=255)
                                                 cttcttcttCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTc                                               >  1:1320723/1‑100 (MQ=255)
                                                   tcttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAc                                             >  1:591934/1‑100 (MQ=255)
                                                   tcttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAc                                             <  2:1189671/100‑1 (MQ=255)
                                                   tcttcttcGATCAGACGCAGAATAGAAAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAATCGTGATGACTGGATTCAc                                             >  2:136428/1‑100 (MQ=255)
                                                    cttgttcGATCAGGCGCAGAAGAGAGAGCGACAGCGATTCGTTGTTACTCACGGTTCCAGTACCAGAACAACTCGGACAATCGTGATGACTGGATTCAcc                                            <  2:1133328/100‑1 (MQ=255)
                                                    cttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAcc                                            <  2:1364905/100‑1 (MQ=255)
                                                     ttcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAccc                                           <  1:2733134/100‑1 (MQ=255)
                                                      tcttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCAccc                                           >  2:433856/1‑99 (MQ=255)
                                                        ttcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGt                                        <  2:2151894/100‑1 (MQ=255)
                                                         tcGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTg                                       <  1:2603360/100‑1 (MQ=255)
                                                         tcGATCAGACGCAGAATAGAGAGCGACAGCGATACGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGCACAAACGTGATGACTGGATTCACCCAGTg                                       >  2:1825362/1‑100 (MQ=255)
                                                          cGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGa                                      <  2:1320711/100‑1 (MQ=255)
                                                           gATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGAt                                     >  2:273786/1‑100 (MQ=255)
                                                            atgagaggaaGAATAGAGGGCGACAGAGATGCGTGGTGACGCACGGTGCCGGTGCCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATg                                    <  2:2455686/91‑1 (MQ=255)
                                                               aGACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGc                                 >  2:116116/1‑100 (MQ=255)
                                                                gACGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGCACAAACGTGATGACTGGATTCACCCAGTGATGGGc                                 >  2:319281/1‑99 (MQ=255)
                                                                 aCGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTc                               <  1:1466805/100‑1 (MQ=255)
                                                                  cGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCa                              <  1:2286729/100‑1 (MQ=255)
                                                                  cGCAGAATAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCa                              <  1:2350425/100‑1 (MQ=255)
                                                                    cAGAATAGAGAGCGACAGCGATTCGTTGTTATGTATGGTGTCAGTAACAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAgn                            <  1:2705906/100‑2 (MQ=255)
                                                                      gAATAGAGAGCGACAGCGATACGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGgcg                          >  2:1847399/1‑100 (MQ=255)
                                                                         tAGAGAGCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTg                       <  2:2078558/100‑1 (MQ=255)
                                                                           gagagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGAc                     <  2:245088/100‑1 (MQ=255)
                                                                            agagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACg                    <  2:399962/100‑1 (MQ=255)
                                                                             gagCGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAAAGTGATGACTGGATTCACCCAGTGATGAGCTCAGGCGCTGACgg                   >  2:1938411/1‑100 (MQ=255)
                                                                                ctaCAGCGATTCGTTGTTACGCACGGTTCCACTACCAGAACAACGCGAACAAACGTGATTACTGGATTCACCCATTGATGGGCTCAGGCGCTGACGGGAc                <  2:349467/98‑1 (MQ=255)
                                                                                cGACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGAc                <  1:2961343/100‑1 (MQ=255)
                                                                                 gACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACa               >  1:1567683/1‑100 (MQ=255)
                                                                                 gACAGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACa               >  2:2214835/1‑100 (MQ=255)
                                                                                    aGCGATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACAttt            <  1:791004/100‑1 (MQ=255)
                                                                                       gATTCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCca         <  2:1551082/100‑1 (MQ=255)
                                                                                       aaTTCGTGGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCTGGCGCTGACGGGACATTTCca         <  2:1424469/99‑1 (MQ=255)
                                                                                        aTTCGTTGTCACGCACGGTGCCAGTACCAGAATAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcag        <  1:2488180/100‑1 (MQ=255)
                                                                                         ttCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGAGATTTCcagc       <  2:402948/100‑1 (MQ=255)
                                                                                         ttCGTTGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcagc       <  1:1645923/100‑1 (MQ=255)
                                                                                            gttgttACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCcagcag     <  2:152650/99‑1 (MQ=255)
                                                                                              tGTTACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGcc  <  2:1072172/100‑1 (MQ=255)
                                                                                                 |                                                                                                
ATTGGCACAGGAACAATGGCGTGAACTTCCTGGGTGTTCTCTTTCAGCGCTTCTTCTTCGATCAGACGCAGAATAGAGAGCGACAGCGATTCGTTGTCACGCACGGTGCCAGTACCAGAACAACGCGGACAAACGTGATGACTGGATTCACCCAGTGATGGGCTCAGGCGCTGACGGGACATTTCCAGCAGGCC  >  NC_000913/1143028‑1143221

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