Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,440,866 A→G 9.7% L247P (CTG→CCG) ‡ fabB ← 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,440,8660AG 9.7% 246.4 / 14.9 82L247P (CTG→CCG) ‡fabB3‑oxoacyl‑[acyl‑carrier‑protein] synthase I
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base G (5/3);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

TCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAA  >  NC_000913/2440720‑2441012
                                                                                                                                                  |                                                                                                                                                  
tCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAgc                                                                                                                                                  >  4:28456/1‑149 (MQ=255)
    gTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCgt                                                                                                                                              <  3:266704/149‑1 (MQ=255)
     tAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCgtg                                                                                                                                             >  3:603814/1‑149 (MQ=255)
       aTCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGtt                                                                                                                                           >  4:326279/1‑149 (MQ=255)
           aTTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAg                                                                                                                                       <  4:482426/149‑1 (MQ=255)
            ttGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGc                                                                                                                                      <  4:307860/149‑1 (MQ=255)
               gggTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGctct                                                                                                                                   <  2:33688/149‑1 (MQ=255)
               gggTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGctct                                                                                                                                   <  1:205090/149‑1 (MQ=255)
                    tCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAAc                                                                                                                              >  4:476474/1‑149 (MQ=255)
                          ccATGCATCGCCATCTTCATACAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTAc                                                                                                                        <  3:510080/149‑1 (MQ=255)
                              gCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATa                                                                                                                    >  4:475870/1‑149 (MQ=255)
                                   gCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgcc                                                                                                               >  1:524272/1‑149 (MQ=255)
                                    ccATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATAccgccg                                                                                                              >  1:424170/1‑149 (MQ=255)
                                        cTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACcgccgccgc                                                                                                          >  3:147294/1‑149 (MQ=255)
                                        cTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACcgccgccgc                                                                                                          >  3:513520/1‑149 (MQ=255)
                                            aTGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGc                                                                                                      >  3:216185/1‑149 (MQ=255)
                                            aTGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGc                                                                                                      >  3:248950/1‑149 (MQ=255)
                                              gCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCTGGCGCgtggt                                                                                                                                           <  2:38610/110‑3 (MQ=255)
                                              gCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCTGGCGCgtggt                                                                                                                                           >  1:38610/1‑108 (MQ=255)
                                                 gCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATaa                                                                                                 <  3:28456/149‑1 (MQ=255)
                                                     aCTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGaa                                                                                             <  1:82277/149‑1 (MQ=255)
                                                      cTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGaaa                                                                                            <  3:350795/149‑1 (MQ=255)
                                                      cTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGaaa                                                                                            <  1:483418/149‑1 (MQ=255)
                                                       tGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAAc                                                                                           <  4:336099/149‑1 (MQ=255)
                                                            cTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCa                                                                                      <  3:476474/149‑1 (MQ=255)
                                                                  cAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGa                                                                                <  3:383477/149‑1 (MQ=255)
                                                                    gACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGc                                                                              >  2:109201/1‑149 (MQ=255)
                                                                    gACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGc                                                                              >  1:3178/1‑149 (MQ=255)
                                                                    gACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGc                                                                              <  1:560614/149‑1 (MQ=255)
                                                                     aCGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGcg                                                                             <  3:413170/149‑1 (MQ=255)
                                                                          gCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTa                                                                        >  4:141242/1‑149 (MQ=255)
                                                                          gCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTa                                                                        >  4:479585/1‑149 (MQ=255)
                                                                           cAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTaa                                                                       <  3:141242/149‑1 (MQ=255)
                                                                           cAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTaa                                                                       >  2:433306/1‑149 (MQ=255)
                                                                             aCCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGt                                                                     >  3:417100/1‑149 (MQ=255)
                                                                              ccATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCa                                                                                                                                <  1:565195/89‑1 (MQ=255)
                                                                              ccATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCa                                                                                                                                >  2:565195/1‑89 (MQ=255)
                                                                              ccATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACCACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTa                                                                    >  4:279257/1‑149 (MQ=255)
                                                                               cATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTAc                                                                   >  1:287053/1‑149 (MQ=255)
                                                                                   tCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGa                                                               >  1:608882/1‑149 (MQ=255)
                                                                                   tCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGa                                                               >  3:236782/1‑149 (MQ=255)
                                                                                    cTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAg                                                              <  2:424170/149‑1 (MQ=255)
                                                                                    cTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAg                                                              >  3:8216/1‑149 (MQ=255)
                                                                                    cTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCAGCCGCCAGCGATAACGAAACCGTCCCGGTGAGCGTCGTAAGTACGGGAg                                                              >  1:200293/1‑149 (MQ=255)
                                                                                     tGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAgg                                                             <  1:433306/149‑1 (MQ=255)
                                                                                         ccATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCttt                                                         >  1:494670/1‑149 (MQ=255)
                                                                                         ccATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCttt                                                         <  4:603814/149‑1 (MQ=255)
                                                                                          cATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCtttt                                                        <  2:3178/149‑1 (MQ=255)
                                                                                               gagGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCgg                                                   >  1:372956/1‑149 (MQ=255)
                                                                                               gagGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCgg                                                   <  2:287053/149‑1 (MQ=255)
                                                                                                     gcgcCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTc                                             >  1:425171/1‑149 (MQ=255)
                                                                                                      cgcCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCg                                            <  4:227542/149‑1 (MQ=255)
                                                                                                          gTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGt                                        <  1:93265/149‑1 (MQ=255)
                                                                                                            aGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTAt                                      <  2:372956/149‑1 (MQ=255)
                                                                                                               cAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCTGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTta                                   <  3:524955/149‑1 (MQ=255)
                                                                                                               cAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTta                                   <  2:65487/149‑1 (MQ=255)
                                                                                                               cAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTta                                   >  3:34059/1‑149 (MQ=255)
                                                                                                                aaCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGc                                                            >  1:419829/1‑123 (MQ=255)
                                                                                                                aaCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGc                                                            <  2:419829/123‑1 (MQ=255)
                                                                                                                aaCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtag                                  <  3:326279/149‑1 (MQ=255)
                                                                                                                 aCGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtagt                                 >  2:311253/1‑149 (MQ=255)
                                                                                                                  cGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGCGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtagta                                <  4:598163/149‑1 (MQ=255)
                                                                                                                   gATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTtagtag                               >  3:576635/1‑149 (MQ=255)
                                                                                                                           cATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCa                       >  1:560578/1‑149 (MQ=255)
                                                                                                                                 tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCAtt                 <  1:311253/149‑1 (MQ=255)
                                                                                                                                 tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCAtt                 <  4:513520/149‑1 (MQ=255)
                                                                                                                                 tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCAtt                 <  4:147294/149‑1 (MQ=255)
                                                                                                                                 tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCAtt                 <  2:425171/149‑1 (MQ=255)
                                                                                                                                    gAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCg              >  1:315082/1‑149 (MQ=255)
                                                                                                                                     tggaacacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGt             >  3:569674/6‑149 (MQ=255)
                                                                                                                                      ggaacacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc            <  4:569674/145‑1 (MQ=255)
                                                                                                                                      gCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc            <  4:248950/149‑1 (MQ=255)
                                                                                                                                      gCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc            <  4:216185/149‑1 (MQ=255)
                                                                                                                                         ccacGCGCCAGCGCCTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGaa         >  3:46883/1‑149 (MQ=255)
                                                                                                                                         acacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc            <  1:165833/145‑1 (MQ=255)
                                                                                                                                         acacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTc            >  2:165833/2‑146 (MQ=255)
                                                                                                                                         acacGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGaa         >  4:61067/2‑149 (MQ=255)
                                                                                                                                          cacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAAt        >  3:225442/1‑149 (MQ=255)
                                                                                                                                          cacGCGCCAGCGCGTGGTCCATCTCCTCACCCACTATCACACAGCCGTCGCGAGCGATAACGAAACCGTCACGGCGAGCGTCGTAAGTACGGGCGGCCTTTTCCGGGGGGTCGTTATATTTAGTAGACATCGCACCCATTGCGTCGaaa        >  4:587344/1‑148 (MQ=255)
                                                                                                                                           acGCGCTGGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGCTAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAAtt       >  1:85254/1‑149 (MQ=255)
                                                                                                                                                ccAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACTGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCaa  <  2:608882/149‑1 (MQ=255)
                                                                                                                                                ccAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCaa  <  3:475870/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCAGGTAATCGATTGGGGTATCAACGCCATGCATCGCCATCTTCATGCAGCGTACTGCGCCTTCGCCAGACGGAGCAACCATGTCTGCACCATCAGAGGTTGCGCCGTAGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAA  >  NC_000913/2440720‑2441012

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