Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,414,242 T→G 6.7% intergenic (+51/‑33) yjfC → / → aidB ATP‑Grasp family ATPase/DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,414,2420TG 6.7% 455.6 / 19.9 165intergenic (+51/‑33)yjfC/aidBATP‑Grasp family ATPase/DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (83/70);  new base G (6/5);  total (89/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-02

Only 100 of 166 total aligned reads displayed.

GTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAT  >  NC_000913/4414096‑4414391
                                                                                                                                                   |                                                                                                                                                     
gTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCAAACTAAt                                                                                                                                                    >  2:887811/1‑150 (MQ=255)
 tCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATg                                                                                                                                                   >  1:401236/1‑150 (MQ=255)
 tCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATg                                                                                                                                                   >  1:1283654/1‑150 (MQ=255)
  cGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCAGGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATgg                                                                                                                                                  <  1:919813/150‑1 (MQ=255)
  cGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTATCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATgg                                                                                                                                                  <  1:718593/150‑1 (MQ=255)
     tgatgaAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGa                                                                                                                                               <  1:1183303/150‑1 (MQ=255)
      gatgaAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGAc                                                                                                                                              >  2:713483/1‑150 (MQ=255)
         gaAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGc                                                                                                                                           <  2:1284712/150‑1 (MQ=255)
          aaGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGAGACTGca                                                                                                                                          >  1:1588306/1‑149 (MQ=255)
           aGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCCCGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCa                                                                                                                                         >  2:1681881/1‑150 (MQ=255)
           aGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCAATACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACAATACTAATGGTGACTGCCa                                                                                                                                         >  1:429027/1‑150 (MQ=255)
            gCGTGCGGAATGGGGCTCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCTCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAt                                                                                                                                        <  2:1208784/150‑1 (MQ=255)
            gCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAt                                                                                                                                        <  2:1074044/150‑1 (MQ=255)
             cGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCAtt                                                                                                                                       <  1:874194/150‑1 (MQ=255)
               tGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGa                                                                                                                                     >  2:1485384/1‑150 (MQ=255)
               tGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGa                                                                                                                                     >  2:1476851/1‑150 (MQ=255)
                gCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACAATACTAATGGTGACTGCCATTGAt                                                                                                                                    >  1:1148339/1‑150 (MQ=255)
                  ggAATGGGGATCCGTGAAGATAACCCACTGATCACCAACGCCACCTCCGGTTTCGTTCCGCAGTTCATGGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATgg                                                                                                                                  >  1:145403/1‑150 (MQ=255)
                       ggggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggaggga                                                                                                                             >  2:377948/1‑150 (MQ=255)
                        gggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGCCCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                            >  2:32063/1‑150 (MQ=255)
                          gATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGac                                                                                                                          <  1:580719/150‑1 (MQ=255)
                               gTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGTATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGt                                                                                                                     <  1:1065086/150‑1 (MQ=255)
                                 gAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGc                                                                                                                   >  2:618796/1‑150 (MQ=255)
                                    gATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACt                                                                                                                <  1:1818354/150‑1 (MQ=255)
                                     ctAACCCACTGCTCCCCAAAGCCACCTCACCTTTCGGTCCGCAGTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACAACTGGGGACTGCCATTGATGGAGGGAGACACAGTGCACTg                                                                                                               <  2:1694906/149‑1 (MQ=255)
                                      tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                            <  2:460779/136‑1 (MQ=255)
                                      tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                            >  1:460779/1‑136 (MQ=255)
                                      tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTgg                                                                                                              <  2:1588152/150‑1 (MQ=255)
                                                aTCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcaca                                                                                                    <  2:1014743/150‑1 (MQ=255)
                                                  cACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACAcc                                                                                                  >  2:1300158/1‑150 (MQ=255)
                                                   aCCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCCCTGGCAAACTCGCACCg                                                                                                 >  1:1125080/1‑150 (MQ=255)
                                                   aCCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCg                                                                                                 >  2:1537028/1‑150 (MQ=255)
                                                    ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGt                                                                                                >  2:636821/1‑150 (MQ=255)
                                                     cAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGtt                                                                                               >  2:814970/1‑150 (MQ=255)
                                                       aaGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGtttt                                                                                             <  2:1591598/150‑1 (MQ=255)
                                                        aGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGttttt                                                                                            <  2:1048068/150‑1 (MQ=255)
                                                         gACACCTCACGTTTCGGTCCGCAGTACATTGATTGATAAGAAGTGATACGCAATATCCTTCTGTCATGAATCCATTGCAGTGACCATACTAATGGTGAATTGCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTa                                                                                           <  1:84388/150‑1 (MQ=255)
                                                             ccTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                                       >  1:646960/1‑150 (MQ=255)
                                                             ccTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCACTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                                       >  1:1026069/1‑150 (MQ=255)
                                                                cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                    <  2:920746/150‑1 (MQ=255)
                                                                cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                    >  1:197189/1‑150 (MQ=255)
                                                                cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                    <  2:565728/150‑1 (MQ=255)
                                                                    tttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                                <  1:284988/150‑1 (MQ=255)
                                                                    tttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                                <  2:126223/150‑1 (MQ=255)
                                                                     ttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                               >  1:563663/1‑150 (MQ=255)
                                                                     ttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                               >  1:1212454/1‑150 (MQ=255)
                                                                      tCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAcc                                                                              <  1:1389528/150‑1 (MQ=255)
                                                                        gTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCAt                                                                            >  2:95725/1‑150 (MQ=255)
                                                                           ccGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTaa                                                                         >  1:1236823/1‑150 (MQ=255)
                                                                             gCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                                       >  1:35118/1‑150 (MQ=255)
                                                                             gCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                                       >  2:1707084/1‑150 (MQ=255)
                                                                              cATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATa                                                                      <  1:713483/150‑1 (MQ=255)
                                                                              cATTACATTGCTGGATAAGAATGTTTTAGAAATCTCTATCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATa                                                                      <  1:1604567/150‑1 (MQ=255)
                                                                                ttACATTGCTGGATACGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAATCTCACCCCGTTTTTAATCAACCTATACCATTAAATAAc                                                                    <  1:1261595/150‑1 (MQ=255)
                                                                                 tACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACa                                                                   >  2:532757/1‑150 (MQ=255)
                                                                                  aCATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAg                                                                  <  2:1199886/150‑1 (MQ=255)
                                                                                   cATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGc                                                                 >  2:1635040/1‑150 (MQ=255)
                                                                                        cTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                            <  1:558524/150‑1 (MQ=255)
                                                                                         tGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTta                                                           >  2:497507/1‑150 (MQ=255)
                                                                                               aGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                                     <  1:1528705/150‑1 (MQ=255)
                                                                                                 aaTGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                                   <  1:1730514/150‑1 (MQ=255)
                                                                                                     ttttAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGTGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATAAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                               >  1:90628/1‑150 (MQ=255)
                                                                                                     tattAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                               <  1:42573/148‑1 (MQ=255)
                                                                                                       ttAGCATTCTCTTTATGTCATGTATCCATGGCAGTGACCATACTAATGGTGACTGCCTTTGATGGAGGGAGACACAGTGCACTTGCAAACTCACACCGTTTTTAATCAACCTATACCATTATATAACAGCAACTTATACCTGTCTGATgg                                             <  2:1090528/150‑1 (MQ=255)
                                                                                                       ttAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATgg                                             <  1:259698/150‑1 (MQ=255)
                                                                                                           gCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGAGACAGTGCACTGGCAACCTCACCCCGTTTTTAATCAGCCTATTCCATGAAAGAAGAGCAAGTTATACCTGTGTGATGgcgc                                          >  2:307854/1‑150 (MQ=255)
                                                                                                           gCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgc                                          >  1:992659/1‑150 (MQ=255)
                                                                                                             aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcgc                                        <  2:400408/150‑1 (MQ=255)
                                                                                                              aTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGct                                       >  1:1816057/1‑150 (MQ=255)
                                                                                                                 tctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTg                                    <  1:532757/150‑1 (MQ=255)
                                                                                                                 tctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTCCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTg                                    <  2:109885/150‑1 (MQ=255)
                                                                                                                   tttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCg                                  <  1:923481/150‑1 (MQ=255)
                                                                                                                      cTGTCTTGAATCCATGTCAGTGTCCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCTACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAg                               <  2:429027/150‑1 (MQ=255)
                                                                                                                       tGTCATGAATCCATGGCAGTGACCATACTAATGTTGACTGCCATTGATGGATGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTATATAACTGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGc                              <  1:1540818/150‑1 (MQ=255)
                                                                                                                       tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACCCCGTTTTTAATCACCCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGc                              >  2:43904/1‑150 (MQ=255)
                                                                                                                       tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGc                              <  1:848083/150‑1 (MQ=255)
                                                                                                                        gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAACCTCACCCCGTTTTTACTCAACCGATACCATGTAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                             >  1:1101902/1‑150 (MQ=255)
                                                                                                                        gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                             >  1:1765129/1‑150 (MQ=255)
                                                                                                                          cATGAATCCATGGCAGTGACCATACTTATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATTCCTTTAAATAACAGCAACTTATACCTGTCTGTTGGCGCGCTCTGCGAAGCTGt                           <  1:1650307/150‑1 (MQ=255)
                                                                                                                          cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                           >  2:1450866/1‑150 (MQ=255)
                                                                                                                          cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                           >  1:1725939/1‑150 (MQ=255)
                                                                                                                            tccaTCAATGGCAGTCACCATTAGTATGGTCAATGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAAATTATACCGGTCTGATGGCGCGCTCTGCGAAGCGGTaa                         >  2:1119174/4‑150 (MQ=38)
                                                                                                                            tGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTaa                         >  2:22638/1‑150 (MQ=255)
                                                                                                                              aaTCCATGGCAGTGGCCATACTAATGGTGACTGCCATTGATGGATGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcg                       <  2:1214417/150‑1 (MQ=255)
                                                                                                                               aTCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACCGCAACTTATACCTGTATGATGGCGCGCTCTGCGAAGCGGTAAcgc                      >  1:349739/1‑150 (MQ=255)
                                                                                                                                tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcgcg                     <  2:1798888/150‑1 (MQ=255)
                                                                                                                                 ccATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGt                    >  2:1273615/1‑150 (MQ=255)
                                                                                                                                    tGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGaa                 <  1:1632471/150‑1 (MQ=255)
                                                                                                                                    tGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGaa                 <  2:607741/150‑1 (MQ=255)
                                                                                                                                            aCCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACCGGCAAACTCCCACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctgg         >  1:952315/1‑150 (MQ=255)
                                                                                                                                              cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct       >  1:512105/1‑150 (MQ=255)
                                                                                                                                              cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct       <  2:966492/150‑1 (MQ=255)
                                                                                                                                               aTACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg      >  2:280045/1‑150 (MQ=255)
                                                                                                                                                tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGTTgctgtctag     >  1:260095/1‑148 (MQ=255)
                                                                                                                                                tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg     >  2:1418235/1‑150 (MQ=255)
                                                                                                                                                tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg     >  1:1009265/1‑150 (MQ=255)
                                                                                                                                                tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg     <  1:1537028/150‑1 (MQ=255)
                                                                                                                                                 aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTggg    <  2:197189/150‑1 (MQ=255)
                                                                                                                                                 aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTggg    <  1:83128/150‑1 (MQ=255)
                                                                                                                                                   tAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTTTCTGAGGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAt  <  2:354765/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAT  >  NC_000913/4414096‑4414391

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