Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,630,663 G→A 6.1% A48T (GCC→ACC)  putA → proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,630,6630GA 6.1% 219.9 / 2.8 76A48T (GCC→ACC) putAproline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base A (3/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATGGCGACTACCACCCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTGCGCTCGGCGATCACCGCCGCCTA  >  NC_002947/5630519‑5630811
                                                                                                                                                |                                                                                                                                                    
gAAATGGCGACTACCACCCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCAcc                                                                                                                                                 <  1:952640/150‑1 (MQ=255)
gAAATGGCGACTACCACCCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCAcc                                                                                                                                                 >  2:2169542/1‑150 (MQ=255)
  aaTGGCGACTACCACCCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCt                                                                                                                                               <  2:1547017/150‑1 (MQ=255)
    tGGCGACTACCACCCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGa                                                                                                                                             >  2:118305/1‑150 (MQ=255)
           tACCACCCTTGGGGCCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACt                                                                                                                                      >  2:1295916/1‑150 (MQ=255)
            accaccCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTc                                                                                                                                     <  2:746544/150‑1 (MQ=255)
             ccaccCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGCGGTGCCACCCTGACCGAACTCa                                                                                                                                    >  1:1415197/1‑150 (MQ=255)
                 ccTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCTATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACgg                                                                                                                                <  2:1559484/150‑1 (MQ=255)
                   ttGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTc                                                                                                                              >  2:1048599/1‑150 (MQ=255)
                    tGGGGTCAAACTCGACTACCCGACCCGTGAGCGCATCAAAGCAGCTGCGCAGTCAATCGCCCACACGCCTACCACGTTGATCAAGCGAACGATCGGCAGTTACCGTGAGGGGCACGAGGGTGGTTCCACCCTGACCGACCTCAACGATCa                                                                                                                             >  2:879893/1‑150 (MQ=255)
                       ggTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATga                                                                                                                          >  2:1261529/1‑149 (MQ=255)
                         tCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCa                                                                                                                        <  1:1295916/150‑1 (MQ=255)
                           aaaCTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGc                                                                                                                      >  2:623957/1‑150 (MQ=255)
                           aaaCTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGc                                                                                                                      >  2:1575404/1‑150 (MQ=255)
                            aaCTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCa                                                                                                                     <  1:623957/150‑1 (MQ=255)
                             aCTCGACGACCCGCCCCGTGACCGACTCCAAGCCGCTGCGCCGTCCAACCACCGCCCGCATCCCTGGTTGATCAAGCAAGCGCACTTCAATTACCTGGAGAAGCCCGAGGGTGATACCACCCTAACCGAACTCAACGGACATGCCACCaa                                                                                                                    >  2:1543482/1‑150 (MQ=255)
                                cgacgaCCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGTAGAAGCTCGAGGGTGGTGCCCACCTGACCGAAGTCAACGGTCATGCCAGCAAccc                                                                                                                 <  1:1261529/150‑1 (MQ=255)
                                 gacgacCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCg                                                                                                                >  1:461663/1‑150 (MQ=255)
                                       ccgacccgTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCcgac                                                                                                          >  1:1813507/1‑150 (MQ=255)
                                       ccgacccgTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCcgac                                                                                                          <  2:461663/150‑1 (MQ=255)
                                         gacccgTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCcgacga                                                                                                        <  1:118305/150‑1 (MQ=255)
                                         gacccgTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCcgacga                                                                                                        <  1:1048599/150‑1 (MQ=255)
                                                   cGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGaa                                                                                              <  1:1864838/150‑1 (MQ=255)
                                                    gACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAg                                                                                             >  1:1601716/1‑150 (MQ=255)
                                                    gACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGACCGACGATGCCGGTGAAg                                                                                             >  2:2103915/1‑150 (MQ=255)
                                                      cTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTc                                                                                           >  1:518407/1‑150 (MQ=255)
                                                        cAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGCGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCa                                                                                         <  2:1415197/150‑1 (MQ=255)
                                                         aaaGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAg                                                                                        >  2:2235111/1‑150 (MQ=255)
                                                          cagcagcTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAgg                                                                                       >  2:1658131/2‑150 (MQ=255)
                                                          aaGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAgg                                                                                       >  2:50906/1‑150 (MQ=255)
                                                            gcagcTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGcc                                                                                     <  2:1441295/150‑1 (MQ=255)
                                                               gcTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGAc                                                                                  <  1:2169542/150‑1 (MQ=255)
                                                                 tGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGCACCACCCtcga                                                                                                                                            <  2:1040050/90‑4 (MQ=255)
                                                                 tGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGCACCACCCtcga                                                                                                                                            >  1:1040050/1‑87 (MQ=255)
                                                                   cgcAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACcacg                                                                              <  2:462868/150‑2 (MQ=255)
                                                                     cAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGc                                                                            <  1:388218/150‑1 (MQ=255)
                                                                      aGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGcc                                                                           >  2:2071468/1‑150 (MQ=255)
                                                                        tCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGccac                                                                         <  1:273669/150‑1 (MQ=255)
                                                                             cGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGGCCACAGCCACCAGTg                                                                    <  2:1601716/150‑1 (MQ=255)
                                                                                 cGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTTGAGAATCTCGAGGGTGGTTCTACCCTGACCGATCTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCGTGCCACAGCCACCAGTGCTTc                                                                >  2:606128/1‑150 (MQ=255)
                                                                                 cGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGCGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTc                                                                >  1:1499347/1‑150 (MQ=255)
                                                                                    aCGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTc                                                             <  1:377208/150‑1 (MQ=255)
                                                                                    aCGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGCGGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTc                                                             <  1:1234611/150‑1 (MQ=255)
                                                                                       ccTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGaa                                                          >  2:1149043/1‑150 (MQ=255)
                                                                                       ccTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGaa                                                          <  1:131507/150‑1 (MQ=255)
                                                                                        cTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAAt                                                         <  2:2174022/150‑1 (MQ=255)
                                                                                         tCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAAtt                                                        >  2:1500567/1‑150 (MQ=255)
                                                                                          cACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAAccc                                                                                                                 <  2:1376598/92‑1 (MQ=255)
                                                                                          cACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGCTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAAccc                                                                                                                 >  1:1376598/1‑92 (MQ=255)
                                                                                            cTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGc                                                     >  1:2021614/1‑150 (MQ=255)
                                                                                               gTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCGCAGCCACCAGTGCTTCCTCGAATTTGCCga                                                  <  1:1575404/150‑1 (MQ=255)
                                                                                               gTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCga                                                  <  2:2021614/150‑1 (MQ=255)
                                                                                               gTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGACCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCga                                                  <  1:2103915/150‑1 (MQ=255)
                                                                                                 tGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCgaga                                                <  1:1658131/150‑1 (MQ=255)
                                                                                                  gATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGagag                                               >  2:719577/1‑150 (MQ=255)
                                                                                                  gATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGagag                                               >  2:874330/1‑150 (MQ=255)
                                                                                                    tCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCTAATTTGCCGAGATCa                                             >  2:420208/1‑150 (MQ=255)
                                                                                                      aagcaagcGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATc                                           <  1:719577/150‑1 (MQ=255)
                                                                                                          aagcGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGc                                       >  2:1379459/1‑150 (MQ=255)
                                                                                                           agcGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCACCCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGcc                                      <  1:1529040/150‑1 (MQ=255)
                                                                                                           agcGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGcc                                      <  1:1149043/150‑1 (MQ=255)
                                                                                                              gATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGt                                                                     <  1:277808/116‑1 (MQ=255)
                                                                                                              gATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGt                                                                     >  2:277808/1‑116 (MQ=255)
                                                                                                                  ttCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCg                               <  1:874330/150‑1 (MQ=255)
                                                                                                                   tCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCgg                              <  1:2235111/150‑1 (MQ=255)
                                                                                                                     aaTTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGc                                                                                      <  1:1439757/92‑1 (MQ=255)
                                                                                                                     aaTTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGc                                                                                      >  2:1439757/1‑92 (MQ=255)
                                                                                                                           cTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTgcgc                      >  2:1701932/1‑150 (MQ=255)
                                                                                                                                gaAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCgg                              >  1:776436/1‑137 (MQ=255)
                                                                                                                                gaAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCgg                              <  2:776436/137‑1 (MQ=255)
                                                                                                                                gaAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTGCGCTCGGc                 <  2:1813507/150‑1 (MQ=255)
                                                                                                                                   gCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTcc                                                                                          <  1:983825/74‑1 (MQ=255)
                                                                                                                                   gCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTcc                                                                                          >  2:983825/1‑74 (MQ=255)
                                                                                                                                      tcaGGGTGGCACCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGTTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTgcg                       >  1:2072384/3‑138 (MQ=255)
                                                                                                                                      tcaGGGTGGCACCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGTTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTgcg                       <  2:2072384/136‑1 (MQ=255)
                                                                                                                                        aGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTGCGCTCGGCGATCAccg         <  2:1499347/150‑1 (MQ=255)
                                                                                                                                               tgCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTGCGCTCGGCGATCACCGCCGCCTa  <  1:50906/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GAAATGGCGACTACCACCCTTGGGGTCAAACTCGACGACCCGACCCGTGAGCGACTCAAAGCAGCTGCGCAGTCCATCGACCGCACGCCTCACTGGTTGATCAAGCAAGCGATCTTCAATTACCTGGAGAAGCTCGAGGGTGGTGCCACCCTGACCGAACTCAACGGTCATGCCAGCAACCCGGCCGACGATGCCGGTGAAGTCCAGGCCGACCACAGCCACCAGTGCTTCCTCGAATTTGCCGAGAGCATCTTGCCGCAGTCGGTACTGCGCTCGGCGATCACCGCCGCCTA  >  NC_002947/5630519‑5630811

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: