Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,900 T→A 5.5% V49D (GTC→GAC)  evgA → response regulator in two‑component regulatory system with EvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,9000TA 5.5% 468.9 / 13.3 127V49D (GTC→GAC) evgAresponse regulator in two‑component regulatory system with EvgS
Reads supporting (aligned to +/- strand):  ref base T (58/62);  new base A (4/3);  total (62/65)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 132 total aligned reads displayed.

ATGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATG  >  NC_000913/2483755‑2484046
                                                                                                                                                 |                                                                                                                                                  
    aCGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCAtt                                                                                                                                             >  4:466566/1‑149 (MQ=255)
       cAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGAt                                                                                                                                          >  1:92892/1‑149 (MQ=255)
        aataatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATg                                                                                                                                         >  3:1082698/1‑149 (MQ=255)
           aatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCg                                                                                                                                      <  3:1397901/149‑1 (MQ=255)
            attattGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                                                     <  2:341131/149‑1 (MQ=255)
                 tgatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATca                                                                                                                                <  1:514541/149‑2 (MQ=255)
                   atgaCCATCCTATTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                                                              >  3:1495085/1‑149 (MQ=255)
                        cATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGt                                                                                                                         >  2:1727698/1‑149 (MQ=255)
                         aTCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGtt                                                                                                                        >  1:1853893/1‑149 (MQ=255)
                          tCCTCTTGCTATCACAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTa                                                                                                                       >  4:666854/1‑149 (MQ=255)
                             tctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                    <  2:503634/149‑1 (MQ=255)
                             tctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcag                                                                                                                               <  3:582663/138‑3 (MQ=255)
                               ttGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGt                                                                                                                  >  4:1871721/1‑149 (MQ=255)
                                 gCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcc                                                                                                                                <  2:1268558/133‑1 (MQ=255)
                                 gCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcc                                                                                                                                >  1:1268558/1‑133 (MQ=255)
                                   tATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATcc                                                                                                              >  1:1523751/1‑149 (MQ=255)
                                     tCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGa                                                                                                                           <  1:345015/134‑1 (MQ=255)
                                     tCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGa                                                                                                                           >  2:345015/1‑134 (MQ=255)
                                     tCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAg                                                                                                            <  2:92892/149‑1 (MQ=255)
                                         agcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtgt                                                                                                        <  1:1727698/149‑1 (MQ=255)
                                          gcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                                                                       <  4:1816803/149‑1 (MQ=255)
                                          gcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                                                                       <  3:1871721/149‑1 (MQ=255)
                                           caATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATg                                                                                                                                         <  3:323038/114‑1 (MQ=255)
                                           caATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                                                                      <  2:1853893/149‑1 (MQ=255)
                                           caATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                                                                      >  3:1279169/1‑149 (MQ=255)
                                             aTTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCACCGGAGTTAACGGTATCCAGGTGTTAGa                                                                                                    >  4:1086358/1‑149 (MQ=255)
                                                 gTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGa                                                                                                                                           <  2:734376/106‑1 (MQ=255)
                                                 gTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGa                                                                                                                                           >  1:734376/1‑106 (MQ=255)
                                                     tttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGa                                                                                            <  3:666854/149‑1 (MQ=255)
                                                     tttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGa                                                                                            <  2:1436524/149‑1 (MQ=255)
                                                      ttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                                                           <  2:1523751/149‑1 (MQ=255)
                                                      ttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                                                           <  4:1495085/149‑1 (MQ=255)
                                                          tGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAg                                                                                       <  1:410315/149‑1 (MQ=255)
                                                           gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                                                     <  1:1600587/102‑1 (MQ=255)
                                                           gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                                                     >  2:1600587/1‑102 (MQ=255)
                                                           gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgc                                                                                      <  3:466566/149‑1 (MQ=255)
                                                              cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                    <  4:724433/116‑1 (MQ=255)
                                                              cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                    >  3:724433/1‑116 (MQ=255)
                                                              cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                                                                   <  2:969196/149‑1 (MQ=255)
                                                               aaaaaCGATATTGACATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                                                                       <  3:1290561/128‑1 (MQ=255)
                                                               aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                                                                       >  4:1290561/1‑128 (MQ=255)
                                                                aaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCaa                                                                                 >  4:1683059/1‑149 (MQ=255)
                                                                         ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGa                                                                        >  2:1082118/1‑149 (MQ=255)
                                                                         ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGa                                                                        >  2:492033/1‑149 (MQ=255)
                                                                          tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                                                     <  1:432696/103‑1 (MQ=255)
                                                                          tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                                                     >  2:432696/1‑103 (MQ=255)
                                                                          tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                                                          >  1:669158/1‑114 (MQ=255)
                                                                          tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                                                          <  2:669158/114‑1 (MQ=255)
                                                                           gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATg                                                                                                                                         >  4:654425/1‑82 (MQ=255)
                                                                           gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATg                                                                                                                                         <  3:654425/82‑1 (MQ=255)
                                                                              aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                                                          <  2:1716890/110‑1 (MQ=255)
                                                                              aTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                                                          >  1:1716890/1‑110 (MQ=255)
                                                                                 ttAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                                                                 >  4:1112741/1‑148 (MQ=255)
                                                                                 ttAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATACCGGAattatta                                                                 <  3:1112741/148‑1 (MQ=255)
                                                                                  tAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGt                                                                                                                                        <  1:1256968/76‑1 (MQ=255)
                                                                                   aGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtgt                                                                                                        <  1:1657563/107‑1 (MQ=255)
                                                                                   aGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtgt                                                                                                        >  2:1657563/1‑107 (MQ=255)
                                                                                    gCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                                                          >  1:49750/1‑104 (MQ=255)
                                                                                         gTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtct                                                        <  3:1683059/149‑1 (MQ=255)
                                                                                               tGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTa                                                  >  3:153102/1‑149 (MQ=255)
                                                                                               tGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTa                                                  >  1:978056/1‑149 (MQ=255)
                                                                                                 aaGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATc                                                                                                               >  4:880480/1‑86 (MQ=255)
                                                                                                 aaGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATc                                                                                                               >  4:387746/1‑86 (MQ=255)
                                                                                                 aaGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATc                                                                                                               <  3:880480/86‑1 (MQ=255)
                                                                                                   ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTaaaaa                                              >  1:480283/1‑149 (MQ=255)
                                                                                                      ggAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGa                                           <  1:1082118/149‑1 (MQ=255)
                                                                                                      ggAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGa                                           <  1:492033/149‑1 (MQ=255)
                                                                                                       gAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAc                                          >  1:214651/1‑149 (MQ=255)
                                                                                                       gAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAc                                          <  2:308204/149‑1 (MQ=255)
                                                                                                           tGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGGTCAtt                                      >  4:1218968/1‑149 (MQ=255)
                                                                                                            gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtctc                                                       <  1:519194/131‑1 (MQ=255)
                                                                                                            gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtctc                                                       >  2:519194/1‑131 (MQ=255)
                                                                                                            gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAttt                                     <  3:1086358/149‑1 (MQ=255)
                                                                                                             ccGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGGTCAtttt                                    >  1:1413123/1‑149 (MQ=255)
                                                                                                               gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                                                               <  2:1378572/88‑1 (MQ=255)
                                                                                                                ttCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTAc                                 <  4:704331/149‑1 (MQ=255)
                                                                                                                ttCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTAc                                 <  2:480283/149‑1 (MQ=255)
                                                                                                                 tCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACg                                <  2:214651/149‑1 (MQ=255)
                                                                                                                    gCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGa                             >  2:1172002/1‑149 (MQ=255)
                                                                                                                          ggAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGa                                                                                            >  1:819106/1‑80 (MQ=255)
                                                                                                                          ggAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGa                                                                                            <  2:819106/80‑1 (MQ=255)
                                                                                                                            aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                                                          <  1:1446463/64‑1 (MQ=255)
                                                                                                                            aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                                                          >  2:1446463/1‑64 (MQ=255)
                                                                                                                            aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGa                                           <  4:1412609/127‑1 (MQ=255)
                                                                                                                              acacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGc                   <  2:978056/149‑1 (MQ=255)
                                                                                                                                acTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTa                                  >  2:590275/1‑132 (MQ=255)
                                                                                                                                acTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTa                                  <  1:590275/132‑1 (MQ=255)
                                                                                                                                 cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGa                <  1:299456/149‑1 (MQ=255)
                                                                                                                                  ttAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGAt               >  2:1465607/1‑149 (MQ=255)
                                                                                                                                        cTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGc                                                                           <  1:1488740/83‑1 (MQ=255)
                                                                                                                                        cTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGc                                                                           >  2:1488740/1‑83 (MQ=255)
                                                                                                                                            tatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATtgt                     <  3:716506/133‑1 (MQ=255)
                                                                                                                                             atCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTaa    <  2:448731/149‑1 (MQ=255)
                                                                                                                                             atCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTaa    >  4:917434/1‑149 (MQ=255)
                                                                                                                                             atCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTaa    >  2:1831598/1‑149 (MQ=255)
                                                                                                                                              tCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAAt   >  4:488341/1‑149 (MQ=255)
                                                                                                                                               cGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATg  <  1:1465607/149‑1 (MQ=255)
                                                                                                                                                 tcatcaTTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                                                            <  3:1620043/89‑1 (MQ=255)
                                                                                                                                                 tcatcaTTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTaa    >  2:136741/1‑145 (MQ=255)
                                                                                                                                                 tcatcaTTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTaa    <  1:136741/145‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATG  >  NC_000913/2483755‑2484046

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