Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,041,974 G→T 10.8% G221V (GGT→GTT)  ygfZ → folate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,041,9740GT10.8% 280.3 / 16.7 90G221V (GGT→GTT) ygfZfolate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (36/44);  new base T (5/5);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

TATGCTGACCGATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGCTGAAAATGGGCG  >  NC_000913/3041825‑3042123
                                                                                                                                                     |                                                                                                                                                     
tatGCTGACCGATAAACTGCGCTGTCAAGCGGAACTGAAGAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGTGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCgg                                                                                                                                                       <  2:1241377/150‑1 (MQ=255)
tatGCTGACCGATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCgg                                                                                                                                                       <  1:1607829/150‑1 (MQ=255)
tatGCTGACCGATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCgg                                                                                                                                                       <  2:224078/150‑1 (MQ=255)
  tGCTGACCGATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTa                                                                                                                                                     >  1:706756/1‑150 (MQ=255)
       aCCGATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTTGGCGGTATCAGc                                                                                                                                                <  2:208083/150‑1 (MQ=255)
        ccGTTAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCt                                                                                                                                               >  1:809014/1‑150 (MQ=255)
         cGATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCtt                                                                                                                                              <  1:727249/150‑1 (MQ=255)
          gATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCttt                                                                                                                                             >  1:410504/1‑150 (MQ=255)
          gATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCttt                                                                                                                                             <  2:1503262/150‑1 (MQ=255)
          gATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCttt                                                                                                                                             <  1:884875/150‑1 (MQ=255)
             aaaCTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAg                                                                                                                                          >  2:1562240/1‑150 (MQ=255)
                  gcgcGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAgg                                                                                                                                     <  1:535471/150‑1 (MQ=255)
                        ggAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTta                                                                                                                               <  2:791397/149‑1 (MQ=255)
                         gAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTtat                                                                                                                              >  1:1209814/1‑150 (MQ=255)
                           aGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATAc                                                                                                                            >  1:1335350/1‑150 (MQ=255)
                           aGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATAc                                                                                                                            >  1:619630/1‑150 (MQ=255)
                           aGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATAc                                                                                                                            >  2:907931/1‑150 (MQ=255)
                              ggAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCgg                                                                                                                         <  1:174579/150‑1 (MQ=255)
                                aaCTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGAc                                                                                                                       <  2:801347/150‑1 (MQ=255)
                                   tGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAag                                                                                                                    <  2:1177754/150‑1 (MQ=255)
                                   tGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAag                                                                                                                    >  1:463572/1‑150 (MQ=255)
                                      aCAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAgaga                                                                                                                 >  2:76417/1‑150 (MQ=255)
                                      aCAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAgaga                                                                                                                 <  1:1473404/150‑1 (MQ=255)
                                      aCAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAgaga                                                                                                                 >  1:1108889/1‑150 (MQ=255)
                                         aTAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGAtgg                                                                                                              <  1:431024/150‑1 (MQ=255)
                                         aTAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGAtgg                                                                                                              >  2:636607/1‑150 (MQ=255)
                                         aTAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGAtgg                                                                                                              >  1:1197179/1‑150 (MQ=255)
                                              cAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGgcg                                                                                                         <  2:1567973/150‑1 (MQ=255)
                                              cAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCAAACAGCGGGCAGATTATCCCACAGGCGACCAAACTCCAGACGCTCGCCGGTATAAGCTTTAAGAAACGCTGTAATGCCGGAAAAGTGAAGCTGgcg                                                                                                         >  2:688015/1‑150 (MQ=255)
                                                aCAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGcgcg                                                                                                       >  2:281138/1‑150 (MQ=255)
                                                   gTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCGGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGc                                                                                                    <  1:464745/150‑1 (MQ=255)
                                                      gctggcATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCATAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCaaa                                                                                                 <  1:1141160/150‑1 (MQ=255)
                                                             ttAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCgt                                                                                          <  2:323015/150‑1 (MQ=255)
                                                              tAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCgtg                                                                                         >  1:419073/1‑150 (MQ=255)
                                                              tAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCgtg                                                                                         >  1:1486904/1‑150 (MQ=255)
                                                              tAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCgtg                                                                                         <  2:1209814/150‑1 (MQ=255)
                                                               aaaCATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCgtgg                                                                                        >  1:1276494/1‑150 (MQ=255)
                                                                aaCATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCgtggt                                                                                       <  2:893306/150‑1 (MQ=255)
                                                                aaCATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCgtggt                                                                                       <  1:1562240/150‑1 (MQ=255)
                                                                  cATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGc                                                                                     <  2:419073/150‑1 (MQ=255)
                                                                  cATTGAAGCGGGGTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTGTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGc                                                                                     <  2:1057064/150‑1 (MQ=255)
                                                                   aTTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGcc                                                                                    >  1:1185622/1‑150 (MQ=255)
                                                                   aTTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGcc                                                                                    <  1:53329/150‑1 (MQ=255)
                                                                    tttaaGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTAGCCCTCAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGAGAGCAAAATTCCGTGGTGCCa                                                                                   <  2:652338/147‑1 (MQ=255)
                                                                       aaGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATa                                                                                <  1:1368655/150‑1 (MQ=255)
                                                                         gCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATaaa                                                                              >  1:1366562/1‑150 (MQ=255)
                                                                         gCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATaaa                                                                              >  2:1165405/1‑150 (MQ=255)
                                                                          cGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAAc                                                                             <  2:410504/150‑1 (MQ=255)
                                                                             gTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTg                                                                          <  2:313439/150‑1 (MQ=255)
                                                                               ttCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTgcg                                                                        <  2:1185622/150‑1 (MQ=255)
                                                                                tCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTgcgc                                                                       <  1:636607/150‑1 (MQ=255)
                                                                                 cccGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGct                                                                      >  2:1115528/1‑150 (MQ=255)
                                                                                    ggTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTg                                                                   <  2:1335350/150‑1 (MQ=255)
                                                                                             tGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGc                                                          >  2:497249/1‑150 (MQ=255)
                                                                                              gccgccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCa                                                         >  1:1060558/1‑150 (MQ=255)
                                                                                               ccgccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAg                                                        <  2:815871/150‑1 (MQ=255)
                                                                                               ccgccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAg                                                        >  1:154682/1‑150 (MQ=255)
                                                                                                cgccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAgg                                                       >  2:1087317/1‑150 (MQ=255)
                                                                                                cgccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAgg                                                       >  2:1240686/1‑150 (MQ=255)
                                                                                                cgccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAgg                                                       <  2:1348749/150‑1 (MQ=255)
                                                                                                 gccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGt                                                      <  2:1486904/150‑1 (MQ=255)
                                                                                                 gccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGt                                                      >  1:1313534/1‑150 (MQ=255)
                                                                                                  ccAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTa                                                     >  2:1593237/1‑150 (MQ=255)
                                                                                                     aCAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAgcg                                                  >  2:289948/1‑150 (MQ=255)
                                                                                                       aGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGcc                                                <  1:187381/150‑1 (MQ=255)
                                                                                                               gTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACt                                        <  2:463572/150‑1 (MQ=255)
                                                                                                                tttATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGATGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTg                                       >  2:1446073/1‑150 (MQ=255)
                                                                                                                tttATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTg                                       >  1:1371897/1‑150 (MQ=255)
                                                                                                                  tATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGcc                                     >  1:1171874/1‑150 (MQ=255)
                                                                                                                   aTCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCg                                    <  2:1481683/150‑1 (MQ=255)
                                                                                                                    tCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCgg                                   <  1:1033182/150‑1 (MQ=255)
                                                                                                                    tCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCgg                                   <  2:809014/150‑1 (MQ=255)
                                                                                                                     cccACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGa                                  >  2:68466/1‑150 (MQ=255)
                                                                                                                       cacaGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAg                                >  2:379307/1‑150 (MQ=255)
                                                                                                                           ggCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTgg                            >  1:1619017/1‑150 (MQ=255)
                                                                                                                           ggCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTgg                            >  1:1015123/1‑150 (MQ=255)
                                                                                                                             cGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTg                          <  1:341633/150‑1 (MQ=255)
                                                                                                                             cGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTg                          <  2:677732/150‑1 (MQ=255)
                                                                                                                                ccAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAg                       <  2:1313534/150‑1 (MQ=255)
                                                                                                                                ccAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAg                       <  2:1542275/150‑1 (MQ=255)
                                                                                                                                  aaTCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAACGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGAc                     <  2:1355301/150‑1 (MQ=255)
                                                                                                                                  aaTCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGAc                     <  2:44762/150‑1 (MQ=255)
                                                                                                                                    tctcCAGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACtt                   <  1:76417/150‑1 (MQ=255)
                                                                                                                                     ctcCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTa                  <  1:907931/150‑1 (MQ=255)
                                                                                                                                      tcCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTag                 <  1:935744/150‑1 (MQ=255)
                                                                                                                                       ccAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTaga                >  2:1434768/1‑150 (MQ=255)
                                                                                                                                       ccAGGCGCTGGGCGGTATAAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTaga                >  2:1349965/1‑150 (MQ=255)
                                                                                                                                       caaGGCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTagc                >  1:977691/3‑149 (MQ=255)
                                                                                                                                         atgcgcTGGGCGGTATCAGCTTTAATAAAGTCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGc              <  1:688015/148‑1 (MQ=255)
                                                                                                                                         aGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGc              <  1:289948/150‑1 (MQ=255)
                                                                                                                                          ggCGCTGGGCGTTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGCt             <  1:174544/150‑1 (MQ=255)
                                                                                                                                          ggCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGCt             >  1:949029/1‑150 (MQ=255)
                                                                                                                                                  gCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGCTGAAAATgg     <  2:619630/150‑1 (MQ=255)
                                                                                                                                                    ggTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGCTGAAAATGGGc   >  1:1257897/1‑150 (MQ=255)
                                                                                                                                                     gTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGCTGAAAATGGGCg  <  2:1171874/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATGCTGACCGATAAACTGCGCGGTGAAGCGGAACTGAACAATAGCCAACAGTGGCTGGCATTAAACATTGAAGCGGGTTTCCCGGTGATTGATGCCGCCAACAGCGGGCAGTTTATCCCACAGGCGACCAATCTCCAGGCGCTGGGCGGTATCAGCTTTAAGAAAGGCTGTTATACCGGACAAGAGATGGTGGCGCGAGCAAAATTCCGTGGTGCCAATAAACGTGCGCTCTGGTTGCTGGCAGGTAGCGCCAGCCGACTGCCGGAAGCTGGTGAAGACTTAGAGCTGAAAATGGGCG  >  NC_000913/3041825‑3042123

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