Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,572,928 C→G 15.9% A759P (GCA→CCA)  pqqL ← putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,9280CG15.9% 368.3 / 60.2 138A759P (GCA→CCA) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base C (48/68);  new base G (9/13);  total (57/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

Only 100 of 146 total aligned reads displayed.

GATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGA  >  NC_000913/1572829‑1573027
                                                                                                   |                                                                                                   
gATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGc                                                                                                     >  1:2796635/1‑100 (MQ=255)
 aTCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGcc                                                                                                    <  2:387194/100‑1 (MQ=255)
  tCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGc                                                                                                     >  1:3065810/1‑98 (MQ=255)
   cTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcac                                                                                                  <  2:903909/100‑1 (MQ=255)
    ttttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGAGAAACCTGTGCcaca                                                                                                 <  1:2640310/100‑1 (MQ=255)
     tttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAg                                                                                                <  2:2297907/100‑1 (MQ=255)
      ttGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAgg                                                                                               <  2:439307/100‑1 (MQ=255)
       tGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATCACGCTTCCACTGTGAAACCTGTGCCACAgg                                                                                               >  1:1404926/1‑99 (MQ=255)
        gCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGtt                                                                                             >  2:547833/1‑100 (MQ=255)
        gCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGtt                                                                                             <  1:2732692/100‑1 (MQ=255)
         ccAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGttt                                                                                            <  1:631925/100‑2 (MQ=255)
           aGTGCGCCGTTAAAAGCATCGAGCGCCCTACGCGTCGGCAGCTTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttca                                                                                          <  1:1651310/100‑4 (MQ=255)
            gTGCGACGTTAAAAGTATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                         <  2:3755604/100‑5 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                         <  2:1545153/100‑5 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                                                                         >  2:147318/1‑100 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                                                                         <  1:1821448/100‑1 (MQ=255)
            gTGCGACGTGAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTGCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                         <  1:3271669/100‑5 (MQ=255)
             tGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcaca                                                                                        <  1:2636925/100‑6 (MQ=255)
             tGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcaca                                                                                        <  1:1171817/100‑6 (MQ=255)
              gCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                       >  2:25881/1‑100 (MQ=255)
               cGACGTTAAAAGCATCTAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTg                                                                                      <  1:3546155/100‑1 (MQ=255)
               cGACGTTAAAAGCATCGAGCGCCATACTCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacag                                                                                       <  2:297240/99‑7 (MQ=255)
               cGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacagt                                                                                      <  2:2528949/100‑8 (MQ=255)
               cGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacagt                                                                                      <  1:749374/100‑8 (MQ=255)
                gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacagtg                                                                                     <  2:3751725/100‑9 (MQ=255)
                  cGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTc                                                                                   <  1:3110488/100‑1 (MQ=255)
                   gtaaaaaGCACCGAACGCAAAACGATCCCGCAAATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCGGCAGGTTGAtttttgttg                                                                                  <  1:1650565/97‑6 (MQ=255)
                   gTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCt                                                                                  >  2:3244416/1‑100 (MQ=255)
                   gTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCt                                                                                  >  2:1608910/1‑100 (MQ=255)
                    ttAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtt                                                                                 >  1:1435364/1‑100 (MQ=255)
                       aacgCATCGCGCGCCATACGCGTCGGCAGATTAGCAGGTGTCCGGGAATCATTACGCTTCTACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtttt                                                                               <  1:2627847/96‑1 (MQ=255)
                       aaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTa                                                                              >  2:1952058/1‑100 (MQ=255)
                        aaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                             <  2:2881493/100‑1 (MQ=255)
                        aaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                             <  1:1504638/100‑1 (MQ=255)
                         aGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                             >  1:439446/1‑99 (MQ=255)
                          gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAg                                                                           <  1:1808617/100‑1 (MQ=255)
                          gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAg                                                                           <  2:895067/100‑1 (MQ=255)
                              cGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAAc                                                                       >  2:2448193/1‑100 (MQ=255)
                              cGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAAc                                                                       >  2:1392071/1‑100 (MQ=255)
                               gAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAAcc                                                                      <  1:148855/100‑1 (MQ=255)
                                 gcgcCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACTCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                    <  1:1740940/100‑1 (MQ=255)
                                 gcgcCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGa                                                                    >  1:2023444/1‑100 (MQ=255)
                                   gcCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCCTTTTGTTCTTTTACAGTAACCGATg                                                                  >  1:2137361/1‑100 (MQ=255)
                                    ccATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGc                                                                 <  2:2631707/100‑1 (MQ=255)
                                     cATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGATCCATTTTGTTCTTTTACAGTAACCGATgtg                                                                >  1:2896984/1‑98 (MQ=255)
                                         cgcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGt                                                            <  1:801943/100‑1 (MQ=255)
                                          gcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTc                                                           <  2:1031407/100‑1 (MQ=255)
                                          gcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTc                                                           <  2:1223419/100‑1 (MQ=255)
                                          gcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTc                                                           <  1:3444392/100‑1 (MQ=255)
                                            gTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCg                                                         <  1:1021506/100‑1 (MQ=255)
                                            gTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCg                                                         >  1:936643/1‑100 (MQ=255)
                                            gTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCg                                                         <  1:378465/100‑1 (MQ=255)
                                            gTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCg                                                         <  2:1300293/100‑1 (MQ=255)
                                             tCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTTCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGt                                                        >  1:2042472/1‑100 (MQ=255)
                                               ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTCCAGTAACCGATGCGTTGTCCGTcg                                                      >  2:3231190/1‑100 (MQ=255)
                                               ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcg                                                      <  2:1060635/100‑1 (MQ=255)
                                                   gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgcg                                                  <  2:1248026/100‑1 (MQ=255)
                                                   gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgcg                                                  >  1:862507/1‑100 (MQ=255)
                                                   gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgcg                                                  <  1:176146/100‑1 (MQ=255)
                                                   gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCCTTTACAGTAACCGCTGCGTTGTCCGCcgcgc                                                   >  1:3158519/1‑99 (MQ=255)
                                                       aaCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGtt                                              >  2:3506066/1‑100 (MQ=255)
                                                       aaCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGtt                                              <  2:2854388/100‑1 (MQ=255)
                                                       aaCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGtt                                              <  2:2753672/100‑1 (MQ=255)
                                                           ggTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                          >  2:3599965/1‑100 (MQ=255)
                                                            gtgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATgg                                         <  2:2181064/100‑1 (MQ=255)
                                                             tgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGt                                        <  2:539697/100‑1 (MQ=255)
                                                              gtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGtt                                       >  2:1846201/1‑100 (MQ=255)
                                                                ccGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTa                                     >  2:1033686/1‑100 (MQ=255)
                                                                ccGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTa                                     >  2:2943564/1‑100 (MQ=255)
                                                                 cGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAc                                    >  2:2525217/1‑100 (MQ=255)
                                                                  gggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                   <  1:3474868/100‑1 (MQ=255)
                                                                  gggAATCATAACGCTTCCACTGTGAAACCTGTACCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                   >  2:3620823/1‑100 (MQ=255)
                                                                     aaTCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGc                                <  1:1716709/100‑1 (MQ=255)
                                                                        cATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGc                             <  2:2132970/100‑1 (MQ=255)
                                                                        cATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGc                             <  1:646170/100‑1 (MQ=255)
                                                                        cATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTGGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGc                             >  1:2083094/1‑100 (MQ=255)
                                                                          tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGc                             <  1:2315376/98‑1 (MQ=255)
                                                                          tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTa                           >  2:1000660/1‑100 (MQ=255)
                                                                             cGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATg                        <  2:1462545/100‑1 (MQ=255)
                                                                               cTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGc                      <  1:2525205/100‑1 (MQ=255)
                                                                                ttCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCg                     <  1:1028834/100‑1 (MQ=255)
                                                                                 tCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGa                    <  2:2443415/100‑1 (MQ=255)
                                                                                  ccACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGaa                   >  1:3255360/1‑100 (MQ=255)
                                                                                  aacaaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGCTGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGc                      >  1:2887022/10‑97 (MQ=38)
                                                                                   acaaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAAt                  >  2:707919/9‑100 (MQ=255)
                                                                                    caaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATc                 >  2:3379584/8‑100 (MQ=255)
                                                                                    aCTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATc                 <  1:2810832/100‑1 (MQ=255)
                                                                                     cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCa                <  2:1641395/100‑1 (MQ=255)
                                                                                     cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCa                >  2:883921/1‑100 (MQ=255)
                                                                                     cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCa                >  2:1962719/1‑100 (MQ=255)
                                                                                     aaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATc                 >  1:3039329/7‑99 (MQ=255)
                                                                                      aaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag               >  1:2280008/6‑100 (MQ=255)
                                                                                      aaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag               >  2:3353075/6‑100 (MQ=255)
                                                                                        tgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCagag             <  1:1020327/100‑1 (MQ=255)
                                                                                         tgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgt            >  2:1679045/3‑100 (MQ=255)
                                                                                          aaaCCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgtg           <  2:2278677/100‑1 (MQ=255)
                                                                                                  gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgattg   >  1:885210/1‑100 (MQ=255)
                                                                                                  gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgattg   <  1:1173889/100‑1 (MQ=255)
                                                                                                  gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgattg   <  2:2019476/100‑1 (MQ=255)
                                                                                                   ccACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTCATGGCGAATCAGAGTGTttgcttg   >  2:1649661/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
GATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGA  >  NC_000913/1572829‑1573027

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