Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,935 T→A 6.0% K302* (AAG→TAG)  nadB ← L‑aspartate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,625,9350TA 6.0% 228.4 / 2.5 77K302* (AAG→TAG) nadBL‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base T (40/31);  new base A (3/2);  total (44/33)
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.86e-01

CGCAGGCACCACCGGGATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAGCTCTTGGC  >  NC_002947/1625792‑1626082
                                                                                                                                                |                                                                                                                                                   
cgcAGGCACCACCGGGATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCAtct                                                                                                                                                <  2:340677/149‑1 (MQ=255)
 gcAGGCACCACCGGGATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCAtctc                                                                                                                                               <  2:39100/149‑1 (MQ=255)
     gCACCACCGGGATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTgg                                                                                                                                           <  2:344237/149‑1 (MQ=255)
         caccGGGATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGa                                                                                                                                       >  1:202701/1‑149 (MQ=255)
         caccGGGATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGa                                                                                                                                       >  1:208437/1‑149 (MQ=255)
              cgATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGCGGAAGTGGCTCTTGATGAAATCTCCAGGCTTGAGAGTGATGTCAACGAATACGAAGTCCACGCCAAGGCCCTACATCTCGTGGTCGATGgcg                                                                                                                                  >  8:279685/2‑149 (MQ=255)
                aTCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGcgcg                                                                                                                                <  2:202701/149‑1 (MQ=255)
                aTCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGcgcg                                                                                                                                <  2:208437/149‑1 (MQ=255)
                            gTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTc                                                                                                                    <  3:20903/149‑1 (MQ=255)
                                 atCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTg                                                                                                               >  3:370877/1‑149 (MQ=255)
                                 aacGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGAGGAAGTTGATCTTGATGAAAACTGCATGCTTGTGAGTGATGTCCAGGAATACGCAGTCCACGCCCAGGCCCTTCATCACGTGGTCGATGGCGCCGCCCACCATGTCCCGTg                                                                                                               >  8:159702/3‑149 (MQ=255)
                                  tCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAAATCTGAAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTg                                                                                                               >  1:96297/1‑149 (MQ=255)
                                    gATGCCCAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCCGGCTTGTGAGTGATGGCCAGGTATCCGCAGTCCACGCCCAGGCGCTGCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCgggggg                                                                                                            >  5:114179/1‑149 (MQ=255)
                                    gATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCgg                                                                                                                               <  4:370162/130‑1 (MQ=255)
                                    gATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCgg                                                                                                                               >  3:370162/1‑130 (MQ=255)
                                        ccAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCg                                                                                                                                              <  8:200631/111‑1 (MQ=255)
                                        ccAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCg                                                                                                                                              >  7:200631/1‑111 (MQ=255)
                                          aaaGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGct                                                                                                      >  2:61565/1‑149 (MQ=255)
                                              gCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTc                                                                                                  >  5:230382/1‑149 (MQ=255)
                                              gCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTc                                                                                                  <  6:230382/149‑1 (MQ=255)
                                                cAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGt                                                                                                                                             >  4:81090/1‑104 (MQ=255)
                                                cAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGt                                                                                                                                             <  3:81090/104‑1 (MQ=255)
                                                  ggCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGt                                                                                              <  2:56273/149‑1 (MQ=255)
                                                  ggCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTACATCTCGTGGt                                                                                                                                          <  2:347069/105‑1 (MQ=255)
                                                  ggCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTACATCTCGTGGt                                                                                                                                          <  2:352631/105‑1 (MQ=255)
                                                  ggCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTACATCTCGTGGt                                                                                                                                          >  1:352631/1‑105 (MQ=255)
                                                  ggCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTACATCTCGTGGt                                                                                                                                          >  1:347069/1‑105 (MQ=255)
                                                            gTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGc                                                                                    >  7:313768/1‑149 (MQ=255)
                                                                 cGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGc                                                                               <  7:279685/149‑1 (MQ=255)
                                                                  ggTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCa                                                                              >  3:76006/1‑149 (MQ=255)
                                                                   gTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCAt                                                                             <  4:76006/149‑1 (MQ=255)
                                                                        gAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAAc                                                                        >  6:71355/1‑149 (MQ=255)
                                                                        gAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAAc                                                                        >  3:180992/1‑149 (MQ=255)
                                                                        gAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAAc                                                                        >  5:219018/1‑149 (MQ=255)
                                                                           gTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGtt                                                                     >  5:151362/1‑149 (MQ=255)
                                                                            tGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTc                                                                    <  1:61565/149‑1 (MQ=255)
                                                                                tctTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGTGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTc                                                                    >  1:292496/1‑145 (MQ=255)
                                                                                tctTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTc                                                                    <  2:292496/145‑1 (MQ=255)
                                                                                tctTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCg                                                                >  3:321215/1‑149 (MQ=255)
                                                                                  ttGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTc                                                                                                  <  3:177475/113‑1 (MQ=255)
                                                                                  ttGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTc                                                                                                  >  4:177475/1‑113 (MQ=255)
                                                                                  ttGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCg                                                                <  4:370877/147‑1 (MQ=255)
                                                                                         aaaTCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAg                                                        <  6:219018/149‑1 (MQ=255)
                                                                                          aaTCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAg                                                        <  4:321215/148‑1 (MQ=255)
                                                                                              tGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGgcg                                                     >  5:94317/1‑147 (MQ=255)
                                                                                                 aGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCgcag                                                  >  8:227898/1‑147 (MQ=255)
                                                                                                   gCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGt                                                                                                                                          >  4:233691/1‑57 (MQ=255)
                                                                                                   gCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGt                                                                                                                                          <  3:233691/57‑1 (MQ=255)
                                                                                                        tgagtgaTGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAgg                                                       >  6:98593/1‑135 (MQ=255)
                                                                                                        tgagtgaTGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAgg                                                       <  5:98593/135‑1 (MQ=255)
                                                                                                          agtgaTGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGccc                                        <  8:313768/148‑1 (MQ=255)
                                                                                                               tGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTcgccg                                  >  2:85589/1‑149 (MQ=255)
                                                                                                                gTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTcgccgc                                 >  1:271738/1‑149 (MQ=255)
                                                                                                                  ccAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCcgcg                                >  6:100222/1‑148 (MQ=255)
                                                                                                                   cAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCcgcg                                >  1:163105/1‑147 (MQ=255)
                                                                                                                   cAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCcgcg                                >  5:288121/1‑147 (MQ=255)
                                                                                                                     ggTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAgg                            <  6:94317/149‑1 (MQ=255)
                                                                                                                           cGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTg                                                                                 <  4:85143/90‑1 (MQ=255)
                                                                                                                           cGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTg                                                                                 >  3:85143/1‑90 (MQ=255)
                                                                                                                           cGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCg                      >  5:153288/1‑149 (MQ=255)
                                                                                                                           cGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCg                      <  2:271738/149‑1 (MQ=255)
                                                                                                                            gCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCgg                     >  1:365942/1‑149 (MQ=255)
                                                                                                                            gCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCACGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCgg                     >  4:248698/1‑149 (MQ=255)
                                                                                                                            gCAGTCCACGCCCAGGCGCTTCATCTCCAGGTCGATCGCCCGGGCCACGATGTCCACTGCGCCCAGCTCTTCGCGGGGCACGAAGCGTGGCATGAAACGATCGCCGTTGGGCACGCGCAGCAGAGCGCCCTAGCCGCGAAGGGCTTCgg                     >  6:117452/1‑149 (MQ=255)
                                                                                                                                  cACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGtt                                                              <  3:331366/102‑1 (MQ=255)
                                                                                                                                  cACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGtt                                                              >  4:331366/1‑102 (MQ=255)
                                                                                                                                  cACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCa               <  6:151362/149‑1 (MQ=255)
                                                                                                                                   aCGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAg              >  7:257649/1‑149 (MQ=255)
                                                                                                                                    cGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGTGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAgg             <  5:71355/149‑1 (MQ=255)
                                                                                                                                     gCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAg              <  5:220839/147‑1 (MQ=255)
                                                                                                                                     gCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAg              >  6:220839/1‑147 (MQ=255)
                                                                                                                                       ccAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAg          >  2:70323/1‑149 (MQ=255)
                                                                                                                                       ccAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAg          >  8:226490/1‑149 (MQ=255)
                                                                                                                                        cAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAGc         >  6:290227/1‑149 (MQ=255)
                                                                                                                                           gcgcTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAGctct      >  3:220626/1‑149 (MQ=255)
                                                                                                                                            cgcTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCAGGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAGCTCtt     <  7:159702/149‑1 (MQ=255)
                                                                                                                                             gcTTCAGCTCGTGGTCGCTGGCGCGGGCCACGAGGTCCCGTGGGGCCAGCTCTTCGCGTGGGGCGAAGCGGGGCATGACCCGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGAGCAGGAAGCTCTTg    <  6:114179/149‑1 (MQ=255)
                                                                                                                                               ttCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAGCTCTTGGc  <  1:85589/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGCAGGCACCACCGGGATCGGCTGACGGGTGATATCGATGCCAAAGGCCAGGCAGCGCTCGTACACGGTGGGGAAGTGGCTCTTGATGAAATCTGCAGGCTTGTGAGTGATGTCCAGGTATACGCAGTCCACGCCCAGGCGCTTCATCTCGTGGTCGATGGCGCGGGCCACGATGTCCCGTGGGGCCAGCTCTTCGCGTGGGTCGAAGCGTGGCATGAAACGTTCGCCGTTGGGCAGGCGCAGCAGGGCGCCCTCGCCGCGCAGGGCTTCGGTGATCAGGAAGCTCTTGGC  >  NC_002947/1625792‑1626082

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