Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,966,460 T→A 6.1% R70R (CGA→CGT PP_2595 ← ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,966,4600TA 6.1% 200.6 / 3.5 66R70R (CGA→CGTPP_2595ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base A (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

CCTGGGTGATCACATGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACTGAATACCCCGTACAGCACGG  >  NC_002947/2966315‑2966602
                                                                                                                                                 |                                                                                                                                              
gcagggTGATCACATGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTGCCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGTCGTCGcc                                                                                                                                            >  1:410299/4‑150 (MQ=255)
 cTGGGTGATCACATGCCCGAGCCGCGCGGTGTTGTCAGCAGACAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCa                                                                                                                                           <  1:374127/150‑1 (MQ=255)
    ggTGATCACATGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGACGCCAggt                                                                                                                                        <  1:543011/150‑2 (MQ=255)
           acaTGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGcac                                                                                                                                 >  1:9103/1‑150 (MQ=255)
           acaTGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTCTTCCACCTGGCGTCGCCAGGCCCAGcac                                                                                                                                 >  1:1047186/1‑150 (MQ=255)
             aTGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAg                                                                                                                               >  1:1306286/1‑150 (MQ=255)
              tGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGc                                                                                                                              <  1:1179885/150‑1 (MQ=255)
              tGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGc                                                                                                                              <  2:1306286/150‑1 (MQ=255)
                   gAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCaa                                                                                                                         <  1:47462/150‑1 (MQ=255)
                   gAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCaa                                                                                                                         >  2:47462/1‑150 (MQ=255)
                        gcgcgGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCcgg                                                                                                                    >  1:69037/1‑150 (MQ=255)
                             gtgtTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGcc                                                                                                               >  2:1258756/1‑150 (MQ=255)
                                  gtCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCAc                                                                                                          <  2:69037/150‑1 (MQ=255)
                                   tCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACg                                                                                                         >  1:192048/1‑150 (MQ=255)
                                    cAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGa                                                                                                        <  1:695869/150‑1 (MQ=255)
                                    cAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGa                                                                                                        <  2:344643/150‑1 (MQ=255)
                                     aGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGac                                                                                                       >  2:231136/1‑150 (MQ=255)
                                         gTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGaccaac                                                                                                   <  2:9103/150‑1 (MQ=255)
                                            aaccaaccCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAg                                                                                                >  1:420648/1‑150 (MQ=255)
                                            aaccaaccCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAg                                                                                                >  1:949003/1‑150 (MQ=255)
                                            aaccaaccCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAg                                                                                                >  2:724981/1‑150 (MQ=255)
                                                  cccAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCa                                                                                                           <  2:477834/133‑1 (MQ=255)
                                                  cccAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCa                                                                                                           >  1:477834/1‑133 (MQ=255)
                                                  cccAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTg                                                                                          <  2:192048/150‑1 (MQ=255)
                                                   ccAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCg                                                                                                                             <  2:1286300/114‑1 (MQ=255)
                                                   ccAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCg                                                                                                                             >  1:1286300/1‑114 (MQ=255)
                                                      gTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAg                                                                                      <  2:949003/150‑1 (MQ=255)
                                                      gTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAg                                                                                      >  2:523553/1‑150 (MQ=255)
                                                                   tACCTGGGCGCCAATGTGCTGGCGCGAGTCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCt                                                                         >  1:749958/1‑150 (MQ=255)
                                                                     ccTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGt                                                                       >  2:1364828/1‑150 (MQ=255)
                                                                      cTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCg                                                                                    >  1:365938/1‑136 (MQ=255)
                                                                      cTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCg                                                                                    <  2:365938/136‑1 (MQ=255)
                                                                         ggCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAg                                                                   <  1:1143925/150‑1 (MQ=255)
                                                                           cgcCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTa                                                                 <  1:523553/150‑1 (MQ=255)
                                                                           cgcCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTa                                                                 <  2:1047186/150‑1 (MQ=255)
                                                                              cAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGt                                                              >  2:1270895/1‑150 (MQ=255)
                                                                                   tgCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTgcagc                                                         >  2:1174104/1‑150 (MQ=255)
                                                                                    gCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTgcagca                                                        <  2:343599/150‑1 (MQ=255)
                                                                                        gcggcCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAg                                                    <  1:1270895/150‑1 (MQ=255)
                                                                                          ggcCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTTGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAgcg                                                  >  2:86813/1‑150 (MQ=255)
                                                                                           gccgccTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCAGTTAAGTAGGTGCAGCATGAgcgg                                                 >  2:578615/1‑150 (MQ=255)
                                                                                             cgccTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTg                                                                        <  2:262824/125‑1 (MQ=255)
                                                                                             cgccTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTg                                                                        >  1:262824/1‑125 (MQ=255)
                                                                                                cTTGCAGCACCGCTACCGCCAGGCTTATGNCGGCCTTTTNCACCNCGCNTCGCCAGGCCCAGCACCGCGTCACGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCa                                            >  1:1293976/1‑150 (MQ=255)
                                                                                                 ttGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAg                                           <  1:724981/150‑1 (MQ=255)
                                                                                                        aCCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTg                                    >  1:888432/1‑150 (MQ=255)
                                                                                                          cGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGAt                                  <  1:1224105/150‑1 (MQ=255)
                                                                                                           gCTACCGCCAGGCTCATGCCCGCCCTTTCCACCTGCCGTCGCCAGGCCCCGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCACCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATg                                 <  2:35982/150‑1 (MQ=255)
                                                                                                                cGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTgc                                                            <  1:199735/118‑1 (MQ=255)
                                                                                                                cGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTgc                                                            >  2:199735/1‑118 (MQ=255)
                                                                                                                cGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAg                            <  2:396284/150‑1 (MQ=255)
                                                                                                                  ccAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGTCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGc                          <  1:372914/150‑1 (MQ=255)
                                                                                                                  ccAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGc                          >  1:1218162/1‑150 (MQ=255)
                                                                                                                    aGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTgca                                                           <  1:922574/115‑1 (MQ=255)
                                                                                                                    aGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTgca                                                           >  2:922574/1‑115 (MQ=255)
                                                                                                                     ggCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCAc                       <  1:1021430/150‑1 (MQ=255)
                                                                                                                      gCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCACGCCGGCGGACACCTGTGCTCTTCAGGTACTGGCAGCGTGTATTTCGCTTTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACt                      <  1:957563/150‑1 (MQ=255)
                                                                                                                          aTGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACTGAAt                  <  1:231136/150‑1 (MQ=255)
                                                                                                                              cGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCt                                                                         >  2:1194007/1‑91 (MQ=255)
                                                                                                                              cGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCt                                                                         <  1:1194007/91‑1 (MQ=255)
                                                                                                                              cGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACTGAATAccc              >  1:1060024/1‑150 (MQ=255)
                                                                                                                                        ggccTGGCGACGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATg                                 <  2:1272034/119‑1 (MQ=255)
                                                                                                                                        ggccTGGCGACGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATg                                 >  1:1272034/3‑121 (MQ=255)
                                                                                                                                        ggccTGGCGACGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACTGAATACCCCGTACAGCAc    >  1:480527/3‑150 (MQ=255)
                                                                                                                                         aCCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACTGAATACCCCGTACAGCACg   <  1:1364828/150‑1 (MQ=255)
                                                                                                                                          ccTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAACATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACTGAATACCCCGTACAGCACgg  <  2:410299/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCTGGGTGATCACATGCCCGAGCCGCGCGGTGTTGTCAGGGGTGAACCAACCCAGTGGCAGCTGTGCTACCTGGGCGCCAATGTGCTGGCGGCCGCCTTGCAGCACCGCTACCGCCAGGCTTATGCCGGCCTTTTCCACCTGGCGTCGCCAGGCCCAGCACAGCGTCAAGCCGGCGGCCAGCACGACCAACCAGGCACTGGCAGCGTGTATTTCGCTGTTCAGTAGGTGCAGCATGAGCGGCGCCAGGGTACTGATGGTCAGGCCACTGAATACCCCGTACAGCACGG  >  NC_002947/2966315‑2966602

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