Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,130,005 G→A 7.0% L57L (TTG→TTA metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,0050GA 7.0% 410.4 / 14.2 129L57L (TTG→TTAmetLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base G (59/61);  new base A (5/4);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 130 total aligned reads displayed.

GGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACC  >  NC_000913/4129857‑4130150
                                                                                                                                                    |                                                                                                                                                 
gggCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTg                                                                                                                                                   >  3:994504/1‑149 (MQ=255)
 ggCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGa                                                                                                                                                  <  2:880770/149‑1 (MQ=255)
  gCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTCCTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGAt                                                                                                                                                 <  2:998222/149‑1 (MQ=255)
  gCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGAt                                                                                                                                                 >  1:382151/1‑149 (MQ=255)
        ggTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTg                                                                                                                                           <  4:296761/149‑1 (MQ=255)
          tcgtcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGt                                                                                                                                         <  1:325202/149‑1 (MQ=255)
            gtcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTg                                                                                                                                       >  2:206/1‑149 (MQ=255)
             tcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGttag                                                                                                                                      <  3:759527/149‑3 (MQ=255)
               aGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGaaa                                                                                                                                    <  2:474054/149‑1 (MQ=255)
                     aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAAc                                                                                                                                             <  4:23269/134‑1 (MQ=255)
                     aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAAc                                                                                                                                             >  3:23269/1‑134 (MQ=255)
                     aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGc                                                                                                                              <  2:184889/149‑1 (MQ=255)
                     aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGc                                                                                                                              >  2:921750/1‑149 (MQ=255)
                     aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGc                                                                                                                              <  2:991647/149‑1 (MQ=255)
                     aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGttagt                                                                                                                                     <  1:846363/142‑4 (MQ=255)
                     aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGttagt                                                                                                                                     >  2:846363/1‑139 (MQ=255)
                       aaaTTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCa                                                                                                                            <  1:91490/149‑1 (MQ=255)
                       aaaTTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCAATATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCa                                                                                                                            <  1:1006165/149‑1 (MQ=255)
                         aTTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCTGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGcc                                                                                                                             >  1:857042/1‑146 (MQ=255)
                         aTTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGcc                                                                                                                             <  2:857042/146‑1 (MQ=255)
                         aTTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGa                                                                                                                          <  3:500305/149‑1 (MQ=255)
                               tggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                    <  3:736806/149‑1 (MQ=255)
                               tggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                    <  4:366228/149‑1 (MQ=255)
                                   agtagtCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCCGTTGATTAACTTGTTGAAACTAAGCCAGACCGATCGtct                                                                                                                >  4:657607/1‑149 (MQ=255)
                                    gtagtCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTaa                                                                                                                                >  1:1013229/1‑132 (MQ=255)
                                    gtagtCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTaa                                                                                                                                <  2:1013229/132‑1 (MQ=255)
                                          tGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTgcg                                                                                                         >  3:510116/1‑149 (MQ=255)
                                           ggCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTgcgc                                                                                                        <  3:598110/149‑1 (MQ=255)
                                           ggCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTgcgc                                                                                                        >  1:1077261/1‑149 (MQ=255)
                                            gCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCa                                                                                                       >  2:432226/1‑149 (MQ=255)
                                             cTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACCGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCAt                                                                                                      <  1:621521/149‑1 (MQ=255)
                                              tGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGTCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGTCCGATCGTCTCTCTGCGCATc                                                                                                     >  2:1077222/1‑149 (MQ=255)
                                                 tgtgAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAgg                                                                                                  <  4:627481/149‑1 (MQ=255)
                                                   tgAAGGGTTATTTCCGTGTGGCGGACATTGTGGGGGAGTATTCTCAGCCTGACGATATGGTGGTGGTTTGCGCCGCCGGTAGCAGCACTAACCAGTTGATTAACTGGTGGAAACTAAGCGAGACCGATCGTCTGTCTGCGCATCAGGtt                                                                                                >  2:940577/1‑149 (MQ=255)
                                                     aaGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGttag                                                                                                                                      <  1:363866/109‑3 (MQ=255)
                                                      aGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaa                                                                                             >  1:957235/1‑149 (MQ=255)
                                                        tgtTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGtctctc                                                                                                             <  4:526453/131‑1 (MQ=255)
                                                        tgtTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaaca                                                                                           >  2:567685/1‑149 (MQ=255)
                                                        tgtTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaaca                                                                                           >  2:797531/1‑149 (MQ=255)
                                                        tgtTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaaca                                                                                           >  2:281446/1‑149 (MQ=255)
                                                             tttGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGc                                                                                      >  4:534183/1‑149 (MQ=255)
                                                             tttGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGc                                                                                      <  4:994504/149‑1 (MQ=255)
                                                                   tgtCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTc                                                                                <  1:232005/149‑1 (MQ=255)
                                                                     tCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGc                                                                              >  2:103409/1‑149 (MQ=255)
                                                                           gCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAg                                                                        >  1:981461/1‑149 (MQ=255)
                                                                             aTTATGGCGGAGTACTCTCAGTCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAAtc                                                                                                                                                >  2:494924/1‑74 (MQ=255)
                                                                             aTTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAAtc                                                                                                                                                <  1:494924/75‑2 (MQ=255)
                                                                             aTTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGGAGCACCACTAACCAGTTGAGTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTg                                                                      >  1:305195/1‑149 (MQ=255)
                                                                              ttATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGaaa                                                                                                                                    >  4:939185/1‑86 (MQ=255)
                                                                              ttATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGaaa                                                                                                                                    <  3:939185/86‑1 (MQ=255)
                                                                              ttATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGc                                                                     <  1:469840/149‑1 (MQ=255)
                                                                              ttATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTCCAACAAACGCTGCGTCGCTATCAGTGc                                                                     >  2:764039/1‑149 (MQ=255)
                                                                                 tGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGAt                                                                  <  3:873552/149‑1 (MQ=255)
                                                                                 tGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGAt                                                                  <  3:160356/149‑1 (MQ=255)
                                                                                  ggCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACAGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATc                                                                 <  2:139639/149‑1 (MQ=255)
                                                                                     ggAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGa                                                              <  4:20793/149‑1 (MQ=255)
                                                                                         tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAg                                                          >  4:133705/1‑149 (MQ=255)
                                                                                         tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAg                                                          >  1:839284/1‑149 (MQ=255)
                                                                                         tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAg                                                          >  2:908540/1‑149 (MQ=255)
                                                                                           ctctCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCg                                                        <  4:773819/149‑1 (MQ=255)
                                                                                           ctctCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCg                                                        <  3:792290/149‑1 (MQ=255)
                                                                                            tctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCgg                                                       >  2:608951/1‑149 (MQ=255)
                                                                                            tctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCgg                                                       >  1:965114/1‑149 (MQ=255)
                                                                                             ctcAGCCTGACGATATGATGGTGGTTTCCGCCGGCGGTAGCACCACTAGCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGt                                                      <  4:504717/149‑1 (MQ=255)
                                                                                              tcAGCCTGGCGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTc                                                     >  4:466236/1‑149 (MQ=255)
                                                                                              tcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTc                                                     >  1:235958/1‑149 (MQ=255)
                                                                                              tcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTc                                                     >  1:37617/1‑149 (MQ=255)
                                                                                                  ccTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCt                                                 <  1:758709/149‑1 (MQ=255)
                                                                                                     gACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAgg                                                                                                  <  2:280859/97‑1 (MQ=255)
                                                                                                        gATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGc                                           <  3:235698/149‑1 (MQ=255)
                                                                                                         atatGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCt                                          <  1:453889/149‑1 (MQ=255)
                                                                                                          tatGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTg                                         >  1:995723/1‑149 (MQ=255)
                                                                                                           atgatgGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAg                                                          >  3:68385/1‑131 (MQ=255)
                                                                                                           atgatgGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAg                                                          <  4:68385/131‑1 (MQ=255)
                                                                                                            tgatgGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaa                                       >  2:435491/1‑149 (MQ=255)
                                                                                                             gatgGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAAc                                                                                        <  1:236807/99‑1 (MQ=255)
                                                                                                               tggtggTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaagaa                                    >  4:117175/1‑149 (MQ=255)
                                                                                                                ggtggtTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGaagaag                                   >  3:481247/1‑149 (MQ=255)
                                                                                                                    gtTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGa                               <  1:432226/149‑1 (MQ=255)
                                                                                                                    gtTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGa                               >  4:170178/1‑149 (MQ=255)
                                                                                                                      ttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATa                             <  1:706796/149‑1 (MQ=255)
                                                                                                                              ccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGc                                  <  3:943280/136‑1 (MQ=255)
                                                                                                                              ccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGc                                  >  4:943280/1‑136 (MQ=255)
                                                                                                                                    gCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCt               <  2:887401/149‑1 (MQ=255)
                                                                                                                                      accacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCttt             >  2:1020451/1‑149 (MQ=255)
                                                                                                                                       ctacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGTCCGACCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtttt            <  2:305195/147‑1 (MQ=255)
                                                                                                                                       ccacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCAGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtttt            <  4:625706/149‑1 (MQ=255)
                                                                                                                                        cacTAACCAGTTGATTAACTGGTTGAAACTATGCCAGACCGATCGTCTCTCTGCGCATCATGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTg           <  4:1060643/149‑1 (MQ=255)
                                                                                                                                        cacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGAt                              >  2:697389/1‑130 (MQ=255)
                                                                                                                                        cacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTg           >  3:197165/1‑149 (MQ=255)
                                                                                                                                         acTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtttt            <  1:78213/147‑1 (MQ=255)
                                                                                                                                         acTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtttt            >  2:78213/1‑147 (MQ=255)
                                                                                                                                           tAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCa        >  2:1097633/1‑149 (MQ=255)
                                                                                                                                           tAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCa        <  4:194044/149‑1 (MQ=255)
                                                                                                                                           caaCCAGTTAATCAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCa        >  4:302347/2‑149 (MQ=255)
                                                                                                                                            aaCCAGTTAATCAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAg       >  3:75602/1‑149 (MQ=255)
                                                                                                                                             aCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                    <  1:1057074/39‑1 (MQ=255)
                                                                                                                                             aCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                    >  2:1057074/1‑39 (MQ=255)
                                                                                                                                             aCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGc      <  4:481247/149‑1 (MQ=255)
                                                                                                                                               cAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGa    <  2:985190/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACC  >  NC_000913/4129857‑4130150

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