Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 273,668 G→C 5.9% R34G (CGC→GGC)  metNB ← methionine ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947273,6680GC 5.9% 273.0 / 6.3 85R34G (CGC→GGC) metNBmethionine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (40/40);  new base C (3/2);  total (43/42)
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.48e-01

CAGGCGGTTGATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAGCTTTGTTCAA  >  NC_002947/273519‑273816
                                                                                                                                                      |                                                                                                                                                    
cAGGCGGTTGATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGAGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAgcg                                                                                                                                                      <  1:533984/150‑1 (MQ=255)
cAGGCGGGTGATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATGTCACCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAgcg                                                                                                                                                      >  2:239479/1‑150 (MQ=255)
     ggTTGATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGt                                                                                                                                                 <  2:1864746/150‑1 (MQ=255)
      gTTGATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGCGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTg                                                                                                                                                <  1:1280399/150‑1 (MQ=255)
       ttGATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTgg                                                                                                                                               <  2:1181650/150‑1 (MQ=255)
         gATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGc                                                                                                                                             >  1:1362165/1‑150 (MQ=255)
         gATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCCCACGTgcgc                                                                                                                                             >  2:1406193/1‑147 (MQ=255)
               gcgcAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCAtt                                                                                                                                       <  2:1043257/150‑1 (MQ=255)
               gcgcAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCAtt                                                                                                                                       <  2:1733215/150‑1 (MQ=255)
               gcgcAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCCCACGTGCGCTTCAtc                                                                                                                                       <  1:1406193/150‑2 (MQ=255)
                cgcAGCAGTTATGACTTACCGGCGCCGCTGCGGCCGATCATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGACGCGCACGTGGGCTTCATTg                                                                                                                                      <  1:140743/150‑1 (MQ=255)
                 gcagcaGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTgg                                                                                                                                     <  1:379344/150‑2 (MQ=255)
                             aCTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGGCGCGCACTTGGGCTTCATTGACGTCGGGGCGCa                                                                                                                         >  2:2109488/1‑150 (MQ=255)
                              cTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAg                                                                                                                        <  1:1291252/150‑1 (MQ=255)
                                  ccGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGcc                                                                                                                    >  2:1646252/1‑150 (MQ=255)
                                  ccGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGcc                                                                                                                    >  1:1069605/1‑150 (MQ=255)
                                     gcgcCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTggg                                                                                                                                              >  1:1042974/1‑121 (MQ=255)
                                     gcgcCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTggg                                                                                                                                              <  2:1042974/121‑1 (MQ=255)
                                       gccgcTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGAGCGCGGGGCCTGTTc                                                                                                               >  2:1676482/1‑150 (MQ=255)
                                          gcTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCt                                                                                                            <  2:1571112/150‑1 (MQ=255)
                                          gcTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCt                                                                                                            >  2:616642/1‑150 (MQ=255)
                                           cTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTg                                                                                                           <  2:1304602/150‑1 (MQ=255)
                                           cTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTg                                                                                                           <  1:1786123/150‑1 (MQ=255)
                                            tGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTgg                                                                                                          >  2:1992852/1‑150 (MQ=255)
                                            tGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTgg                                                                                                          <  2:2170770/150‑1 (MQ=255)
                                             gCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGc                                                                                                         <  1:958344/150‑1 (MQ=255)
                                                   gatgatGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGTCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTgg                                                                                                   <  2:1573113/150‑1 (MQ=255)
                                                   gatgatGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTgg                                                                                                   >  1:2363698/1‑150 (MQ=255)
                                                        tgCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTgggggcg                                                                                              >  2:136016/1‑149 (MQ=255)
                                                          ccGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCATGACCGGTAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCTCTGGCCATTGGCGGCAGt                                                                                            >  1:1956913/1‑150 (MQ=255)
                                                          ccGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAgg                                                                                            >  2:2355775/1‑149 (MQ=255)
                                                               gATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGcc                                                                                                                    <  2:148070/121‑1 (MQ=255)
                                                               gATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGcc                                                                                                                    >  1:148070/1‑121 (MQ=255)
                                                                aTCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTgg                                                                                      <  2:661849/150‑1 (MQ=255)
                                                                 tctcGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTggg                                                                                     <  1:628574/150‑1 (MQ=255)
                                                                   tcGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGgt                                                                                   >  2:821026/1‑150 (MQ=255)
                                                                   tcGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGgt                                                                                   >  1:1882449/1‑150 (MQ=255)
                                                                     gCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGgtgt                                                                                 >  2:1561704/1‑150 (MQ=255)
                                                                     gCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGgtgt                                                                                 >  1:1801858/1‑150 (MQ=255)
                                                                      ccATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTa                                                                                >  1:603954/1‑150 (MQ=255)
                                                                       cATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAg                                                                               <  1:51301/150‑1 (MQ=255)
                                                                        aTGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAgg                                                                              <  2:603954/150‑1 (MQ=255)
                                                                             ggATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGccc                                                                         >  2:1062419/1‑150 (MQ=255)
                                                                               aTGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTg                                                                       >  2:2107230/1‑150 (MQ=255)
                                                                                  ttCAGGTCGATGCATTGCAACGCCGGTTGCGCCTGGCCGGGGTAGTTCTTTCCCAGGCGGATGAAGCGCACGTGGGCTTCATTGACGTAGGGGCGCACGGTCGAACAGTTTAGCTTTTGCGGGAGTGCTTGGGGTGTAGGCGCCCTGAgc                                                                    >  2:529513/1‑150 (MQ=255)
                                                                                     aGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTggg                                                                                     >  2:1067634/1‑130 (MQ=255)
                                                                                     aGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTggg                                                                                     <  1:1067634/130‑1 (MQ=255)
                                                                                       gTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGa                                                                      <  2:1261224/143‑1 (MQ=255)
                                                                                       gTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGa                                                                      >  1:1261224/1‑143 (MQ=255)
                                                                                       gTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGct                                                               <  1:1561704/150‑1 (MQ=255)
                                                                                        tCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGt                                                                                                                 >  1:738637/1‑99 (MQ=255)
                                                                                        tCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGt                                                                                                                 <  2:738637/99‑1 (MQ=255)
                                                                                        tCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTAATGCCATGGCCGAAGACGCGCATGTGGGCTTCATTGACGTCGGGGCTCTGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGctg                                                              <  2:2078420/150‑1 (MQ=255)
                                                                                         cGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGctgg                                                             <  2:408932/150‑1 (MQ=255)
                                                                                          gATGCCTTGCAACGCCGGATGCGCCTGGCCGGGGAGGGTCTTTCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTCGCTGCAGTGCTGGGGGTGTAGGCGCCCTGAGCGCGctggc                                                            >  2:1955163/1‑150 (MQ=255)
                                                                                           aTGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCt                                                                                                            <  2:1507635/101‑1 (MQ=255)
                                                                                           aTGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCt                                                                                                            >  1:1507635/1‑101 (MQ=255)
                                                                                            tGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCt                                                                        >  2:1490141/1‑136 (MQ=255)
                                                                                            tGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCt                                                                        <  1:1490141/136‑1 (MQ=255)
                                                                                             gCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGctggcctg                                                         <  2:1956913/150‑1 (MQ=255)
                                                                                               cTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTggg                                                                                     >  1:933337/1‑120 (MQ=255)
                                                                                               cTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTggg                                                                                     <  2:933337/120‑1 (MQ=255)
                                                                                               cTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGctggcctggc                                                       >  1:175143/1‑150 (MQ=255)
                                                                                                   cAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCAt                                                   >  2:1506612/1‑150 (MQ=255)
                                                                                                                gcCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCa                                      >  1:915442/1‑150 (MQ=255)
                                                                                                                   tGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGcc                                   >  2:125989/1‑150 (MQ=255)
                                                                                                                    ggCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTACCGGCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCg                                  <  1:1955163/150‑1 (MQ=255)
                                                                                                                       cGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCt                               <  2:1069605/150‑1 (MQ=255)
                                                                                                                       cGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGCGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCt                               >  2:2086219/1‑150 (MQ=255)
                                                                                                                        ggggTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCtt                              <  1:2157822/150‑1 (MQ=255)
                                                                                                                          ggTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTgg                            <  2:1067650/150‑1 (MQ=255)
                                                                                                                           gTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGCCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGt                           >  2:868149/1‑150 (MQ=255)
                                                                                                                            tAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTa                          <  2:2363698/150‑1 (MQ=255)
                                                                                                                              ggTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTaga                        <  1:929238/150‑1 (MQ=255)
                                                                                                                                tCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGCCGGGGTGAAGGGCCTGTTCCCTGGCCATTTGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCGTCCCAGCCGGCTTGGTAGAGc                      >  1:1023024/1‑150 (MQ=255)
                                                                                                                                  ttGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCtt                    >  2:1788922/1‑150 (MQ=255)
                                                                                                                                   tGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTg                   <  2:1882449/150‑1 (MQ=255)
                                                                                                                                      ccAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCg                >  2:893099/1‑150 (MQ=255)
                                                                                                                                      ccAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCg                >  2:980725/1‑150 (MQ=255)
                                                                                                                                          ggCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGa             >  2:1592385/1‑150 (MQ=255)
                                                                                                                                          ggCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGa             <  1:896533/150‑1 (MQ=255)
                                                                                                                                           gCCGATGAAGCGCACGCGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAg            <  1:2086219/150‑1 (MQ=255)
                                                                                                                                            tcaaTGAAGCCCACGTGCGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAGc           >  1:152546/4‑150 (MQ=255)
                                                                                                                                              gATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAGCtt         >  1:2230381/1‑150 (MQ=255)
                                                                                                                                              aaTGAAGCCCACGTGCGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCa                                                                                              <  1:645406/64‑1 (MQ=38)
                                                                                                                                              aaTGAAGCCCACGGGCGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCa                                                                                              >  2:645406/2‑65 (MQ=255)
                                                                                                                                                tGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAGCTTTg       <  2:175143/150‑1 (MQ=255)
                                                                                                                                                    gcgcACGTGGGCTTCCTAGAAGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAGCTTTGTTCa   >  2:1420179/1‑150 (MQ=255)
                                                                                                                                                     cgcACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAGCTTTGTTCaa  <  1:821026/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAGGCGGTTGATGGTGCGCAGCAGTGATGACTTACCGGCGCCGCTGCGGCCGATGATGCCGAAGATCTCGCCATGGCGGATGTTCAGGTCGATGCCTTGCAACGCCGGTTGCGCCTGGCCGGGGTAGGTCTTGCCCAGGCCGATGAAGCGCACGTGGGCTTCATTGACGTCGGGGCGCAGGGCCTGTTCCCTGGCCATTGGCGGCAGTGCTTGGGGTGTAGGCGCCCTGAGCGCGCTGGCCTGGCTCATTTCAGCCTTCCCAGCCGGCTTGGTAGAGCTTGCCGTGAGCTTTGTTCAA  >  NC_002947/273519‑273816

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