Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,800,844 C→T 6.7% intergenic (+54/+25) PP_3360 → / ← PP_3361 membrane protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,800,8440CT 6.7% 267.1 / 5.8 90intergenic (+54/+25)PP_3360/PP_3361membrane protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (46/37);  new base T (4/2);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

CCGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGATC  >  NC_002947/3800701‑3800992
                                                                                                                                                |                                                                                                                                                    
ccGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACg                                                                                                                                                 >  4:201138/1‑149 (MQ=255)
 cGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGc                                                                                                                                                >  6:335305/1‑149 (MQ=255)
       gcgcgATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGcc                                                                                                                                          <  1:207158/149‑1 (MQ=255)
            aTCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCgcgc                                                                                                                                     <  5:97770/149‑1 (MQ=255)
               aaGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTg                                                                                                                                   >  1:385986/1‑148 (MQ=255)
               aaGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTg                                                                                                                                   >  1:380091/1‑148 (MQ=255)
                 ggTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATCCGCACGCCCTGCCCGCGCTGaaa                                                                                                                                >  2:321713/1‑149 (MQ=255)
                   tCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACg                                                                                                                              >  5:271685/1‑149 (MQ=255)
                     aCGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGGCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGc                                                                                                                             <  3:335138/148‑1 (MQ=255)
                     aCGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGc                                                                                                                             >  5:5074/1‑148 (MQ=255)
                     aCGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGc                                                                                                                             >  6:171190/1‑148 (MQ=255)
                     aCGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCt                                                                                                                            >  5:57068/1‑149 (MQ=255)
                          aTTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCGCATGTGCACGCCCt                                                                                                                                             <  2:322726/127‑1 (MQ=255)
                          aTTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGAGCACCGCAGGGCGTGCGCATGTGCACGCCCt                                                                                                                                             >  1:322726/1‑127 (MQ=255)
                          aTTGCCAGAGGCGTCCGCCCTGCATCGCCAGGCTCGCGAAGGCCTGAGCATCACGGGCGCATGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGAGCACAAGCGCACGCCCTGCCCGCCCTGAAACGCTaata                                                                                                                        >  2:121578/1‑145 (MQ=255)
                            tgcatgagGCGTCCGGCCTGCATGCCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCa                                                                                                                      <  2:380091/143‑1 (MQ=255)
                            tGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCa                                                                                                                      <  2:385986/148‑1 (MQ=255)
                               cAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAACGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGCTAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAgcgc                                                                                                                  >  5:234609/1‑149 (MQ=255)
                                   ggCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCcgcg                                                                                                              <  3:98636/149‑1 (MQ=255)
                                      gTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTg                                                                                                           <  5:171190/149‑1 (MQ=255)
                                        ccGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTAttt                                                                                                                        >  5:2993/1‑134 (MQ=255)
                                        ccGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTAttt                                                                                                                        <  6:2993/134‑1 (MQ=255)
                                         cGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAg                                                                                                        <  6:271685/149‑1 (MQ=255)
                                         cGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAg                                                                                                        >  1:43929/1‑149 (MQ=255)
                                            ccTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTg                                                                                                     <  3:98112/149‑1 (MQ=255)
                                             cTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGt                                                                                                    >  1:343848/1‑149 (MQ=255)
                                              tGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTc                                                                                                   >  4:201259/1‑149 (MQ=255)
                                                cATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCtt                                                                                                 <  8:52616/149‑1 (MQ=255)
                                                cATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCt                                                                                                  <  2:343848/148‑1 (MQ=255)
                                                 aTGGCCAGGCTCGCGAAGGCCTGAGCATCTGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATATGTTCACCGCAGGGCGTGCTCATGGGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCttt                                                                                                <  1:121578/149‑1 (MQ=255)
                                                  tGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAATCTGTCTTTc                                                                                               <  4:320132/149‑1 (MQ=255)
                                                        ggCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCAcc                                                                                         <  6:28117/149‑1 (MQ=255)
                                                        ggCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCAcc                                                                                         >  5:28117/1‑149 (MQ=255)
                                                           tCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCACGCCAGATAGGTGCACCGCAGCGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCCCCCCGGTGAAGCTGTCTTTCAGCACCTc                                                                                       >  5:202877/1‑148 (MQ=255)
                                                             gcgAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCt                                                                                                      <  7:242835/131‑1 (MQ=255)
                                                             gcgAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGACCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCt                                                                                                      >  8:242835/1‑131 (MQ=255)
                                                                 aGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaagg                                                                                <  6:5074/149‑1 (MQ=255)
                                                                 aGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaagg                                                                                <  6:57068/149‑1 (MQ=255)
                                                                 aGGCCGGAGCATCAGGGGCCCTTGATGCGGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaagg                                                                                <  4:203030/149‑1 (MQ=255)
                                                                   gCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaaggaa                                                                              <  4:82302/149‑1 (MQ=255)
                                                                   gCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaagg                                                                                >  1:393291/1‑147 (MQ=255)
                                                                     cTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGt                                                                            <  3:372521/149‑1 (MQ=255)
                                                                            tCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCAtt                                                                     >  8:30341/1‑149 (MQ=255)
                                                                                  gcgcTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCg                                                               >  8:301744/1‑149 (MQ=255)
                                                                                    gcTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGc                                                             <  6:234609/149‑1 (MQ=255)
                                                                                         aTGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaa                                                         <  4:206749/148‑1 (MQ=255)
                                                                                         aTGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaa                                                         >  3:206749/1‑148 (MQ=255)
                                                                                         aTGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATCCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaa                                                        <  1:321713/149‑1 (MQ=255)
                                                                                          tGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaa                                                         >  1:80928/1‑148 (MQ=255)
                                                                                          tGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaaa                                                       >  7:288312/1‑149 (MQ=255)
                                                                                          tGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaaa                                                       <  5:335305/149‑1 (MQ=255)
                                                                                          tGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaaa                                                       >  7:37273/1‑149 (MQ=255)
                                                                                           gCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGaaaa                                                       >  1:98857/1‑148 (MQ=255)
                                                                                            cTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCgg                                                                                                             >  8:324104/1‑93 (MQ=255)
                                                                                            cTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCgg                                                                                                             <  7:324104/93‑1 (MQ=255)
                                                                                             tGGCGACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAatc                                                     >  7:107472/1‑148 (MQ=255)
                                                                                             tGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAatca                                                    <  2:393291/149‑1 (MQ=255)
                                                                                               gCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAatcatc                                                  >  2:10568/1‑149 (MQ=255)
                                                                                                  aCCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGaa                                               <  8:37273/149‑1 (MQ=255)
                                                                                                     tGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGc                                            >  5:208168/1‑149 (MQ=255)
                                                                                                     tGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCCCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGc                                            >  1:352088/1‑149 (MQ=255)
                                                                                                         cGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCCCCTCAAGAAGGCTGTCCGCATTGGCCCGGCTGCAAATCATCGAAGGCCGc                                         >  1:329350/1‑148 (MQ=255)
                                                                                                         cGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCa                                        >  4:165739/1‑149 (MQ=255)
                                                                                                             tCGCCAGATAGGTGCACCGCAGGGCGTGCGCATGCGCGCGCCCTGCCCGCCCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTCTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCGGCAtttt                                    >  4:270552/1‑149 (MQ=255)
                                                                                                               gCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCAtttt                                    >  6:261226/1‑147 (MQ=255)
                                                                                                               gCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCa                                  <  3:201259/149‑1 (MQ=255)
                                                                                                                 cAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAg                                 >  2:66859/1‑148 (MQ=255)
                                                                                                                  aGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCa                               <  2:98857/149‑1 (MQ=255)
                                                                                                                   gATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCAc                              >  7:115170/1‑149 (MQ=255)
                                                                                                                      aGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGtt                           >  7:179827/1‑149 (MQ=255)
                                                                                                                        gTGCAGCGCAGGGCGGGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGATGGAAGTCCGCATTGGCGCGGCTGAAAATCATCGACGGCCGGAtttt                                    <  1:280556/138‑1 (MQ=255)
                                                                                                                        gTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCa                                                                                              >  6:195958/1‑80 (MQ=255)
                                                                                                                        gTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCa                                                                                              <  5:195958/80‑1 (MQ=255)
                                                                                                                        gTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTCTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGCAGGCCCCAattt                                    >  2:280556/1‑138 (MQ=255)
                                                                                                                         tGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTc                        <  3:165739/149‑1 (MQ=255)
                                                                                                                          gCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCa                       >  7:273089/1‑149 (MQ=255)
                                                                                                                          gCACCGCAGGGCGTGCACATGCGCACGCCCTGCACGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGATGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGACCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCa                       >  6:307722/1‑149 (MQ=255)
                                                                                                                             ccGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGc                    <  3:201138/149‑1 (MQ=255)
                                                                                                                              cGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaaggaa                                                                              <  2:360098/90‑1 (MQ=255)
                                                                                                                              cGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAgaaggaa                                                                              >  1:360098/1‑90 (MQ=255)
                                                                                                                                 aGGGCGTGCGCATGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGccc                                                                <  2:110119/101‑1 (MQ=255)
                                                                                                                                 aGGGCGTGCGCATGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGccc                                                                >  1:110119/1‑101 (MQ=255)
                                                                                                                                 aGGGCGTGCGCATGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGAATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCggg                 >  4:299321/1‑148 (MQ=255)
                                                                                                                                 aGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCAGTTTCAGCACGTTGTCATGCGGGc                >  6:225117/1‑149 (MQ=255)
                                                                                                                                    gCGTGCGCATGTGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCtt             >  1:381870/1‑149 (MQ=255)
                                                                                                                                    gCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCtt             <  2:80928/149‑1 (MQ=255)
                                                                                                                                     cGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTc            >  5:19104/1‑149 (MQ=255)
                                                                                                                                           acaTGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGc       <  7:301744/149‑1 (MQ=255)
                                                                                                                                             aTGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGcc      >  2:27729/1‑148 (MQ=255)
                                                                                                                                              tGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGa    <  6:202877/149‑1 (MQ=255)
                                                                                                                                               gcgcACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGAt   >  5:200632/1‑149 (MQ=255)
                                                                                                                                               gcgcACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGAt   >  1:114916/1‑149 (MQ=255)
                                                                                                                                                cgcACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGATc  <  1:66859/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCGGTGCGCGCGATCAAGGTCACGGCATTGCCAGAGGCGTCCGGCCTGCATGGCCAGGCTCGCGAAGGCCTGAGCATCAGGGGCGCTTGATGCTGGCAACCTGCCCGGCTCGCCAGATAGGTGCACCGCAGGGCGTGCACATGCGCACGCCCTGCCCGCGCTGAAACGCTATTTCAGCGCCGCGGTGAAGCTGTCTTTCAGCACCTCAAGAAGGAAGTCCGCATTGGCCCGGCTGAAAATCATCGAAGGCCGCATTTTCAGCACGTTGTCATGCGGGCCTTCGGTGCCGATC  >  NC_002947/3800701‑3800992

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