Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,414,249 G→C 6.6% intergenic (+58/‑26) 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,2490GC 6.6% 516.7 / 21.2 163intergenic (+58/‑26)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 G (80/72);  new base C (6/5);  total (86/77)
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.36e-01

Only 100 of 164 total aligned reads displayed.

TGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAT  >  NC_000913/4414101‑4414397
                                                                                                                                                     |                                                                                                                                                    
tgatgaAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGa                                                                                                                                                     <  1:1183303/150‑1 (MQ=255)
 gatgaAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGAc                                                                                                                                                    >  2:713483/1‑150 (MQ=255)
     aaGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGAGACTGca                                                                                                                                                >  1:1588306/1‑149 (MQ=255)
      aGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCCCGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCa                                                                                                                                               >  2:1681881/1‑150 (MQ=255)
       gCGTGCGGAATGGGGCTCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCTCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAt                                                                                                                                              <  2:1208784/150‑1 (MQ=255)
        cGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCAtt                                                                                                                                             <  1:874194/150‑1 (MQ=255)
        cGTGCGGAATGGGGAGCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTTGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAtg                                                                                                                                             <  2:1664973/150‑2 (MQ=255)
          tGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGa                                                                                                                                           >  2:1476851/1‑150 (MQ=255)
          tGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGa                                                                                                                                           >  2:1485384/1‑150 (MQ=255)
           gCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACAATACTAATGGTGACTGCCATTGAt                                                                                                                                          >  1:1148339/1‑150 (MQ=255)
            cGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATg                                                                                                                                         >  1:1535690/1‑150 (MQ=255)
             ggAATGGGGATCCGTGAAGATAACCCACTGATCACCAACGCCACCTCCGGTTTCGTTCCGCAGTTCATGGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATgg                                                                                                                                        >  1:145403/1‑150 (MQ=255)
                  ggggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTT‑GCAATCTCTTTCTGTCATGAATCAATGGCAGTGACCATACTAATGGTGACTGCCATTGATggaggga                                                                                                                                   >  2:319693/1‑150 (MQ=255)
                   gggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGCCCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                                  >  2:32063/1‑150 (MQ=255)
                     gATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGac                                                                                                                                <  1:580719/150‑1 (MQ=255)
                          gTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGTATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGt                                                                                                                           <  1:1065086/150‑1 (MQ=255)
                            gAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGc                                                                                                                         >  2:618796/1‑150 (MQ=255)
                             aaGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCa                                                                                                                        >  2:106134/1‑150 (MQ=255)
                               gATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACt                                                                                                                      <  1:1818354/150‑1 (MQ=255)
                                 tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                                  >  1:460779/1‑136 (MQ=255)
                                 tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                                  <  2:460779/136‑1 (MQ=255)
                                 tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTgg                                                                                                                    <  2:1588152/150‑1 (MQ=255)
                                 tAAAACACTGATCACCAAAGACACCTCACGTTTAGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTCCCATTGATGGAGGTAGACACAGTGCACTgg                                                                                                                    >  2:380609/1‑150 (MQ=255)
                                             cACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACAcc                                                                                                        >  2:1300158/1‑150 (MQ=255)
                                              aCCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCCCTGGCAAACTCGCACCg                                                                                                       >  1:1125080/1‑150 (MQ=255)
                                              aCCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCg                                                                                                       >  2:1537028/1‑150 (MQ=255)
                                               ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGt                                                                                                      >  2:1765884/1‑150 (MQ=255)
                                               ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGt                                                                                                      >  1:1798888/1‑150 (MQ=255)
                                               ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGt                                                                                                      >  2:636821/1‑150 (MQ=255)
                                                cAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGtt                                                                                                     >  2:814970/1‑150 (MQ=255)
                                                   aGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGttttt                                                                                                  <  2:1048068/150‑1 (MQ=255)
                                                    gACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTa                                                                                                 <  2:451639/150‑1 (MQ=255)
                                                    gACACCTCACGTTTCGGTCCGCAGTACATTGATTGATAAGAA‑GTGATACGCAATATCCTTCTGTCATGAATCCATTGCAGTGACCATACTAATGGTGAATTGCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTa                                                                                                 <  1:84388/150‑1 (MQ=255)
                                                        ccTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCACTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                                             >  1:1026069/1‑150 (MQ=255)
                                                           cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                          >  1:197189/1‑150 (MQ=255)
                                                           cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                          <  2:920746/150‑1 (MQ=255)
                                                           cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                                          <  2:565728/150‑1 (MQ=255)
                                                                ttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                                     >  1:1212454/1‑150 (MQ=255)
                                                                ttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                                     >  1:563663/1‑150 (MQ=255)
                                                                  cGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCa                                                                                   >  2:855851/1‑150 (MQ=255)
                                                                   gTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCAt                                                                                  >  2:95725/1‑150 (MQ=255)
                                                                   gTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCAt                                                                                  >  1:554868/1‑150 (MQ=255)
                                                                      ccGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTaa                                                                               >  1:1236823/1‑150 (MQ=255)
                                                                        gCATTACATTGCTGGATTAGAATTTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATAATAATTGTGACTGACATTGATGGAGGTATACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                                             <  1:260972/150‑1 (MQ=255)
                                                                        gCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                                             >  1:35118/1‑150 (MQ=255)
                                                                         cATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATa                                                                            <  1:713483/150‑1 (MQ=255)
                                                                           ttACATTGCTGGATACGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAATCTCACCCCGTTTTTAATCAACCTATACCATTAAATAAc                                                                          <  1:1261595/150‑1 (MQ=255)
                                                                            tACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACa                                                                         >  2:532757/1‑150 (MQ=255)
                                                                                   cTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                                  <  1:558524/150‑1 (MQ=255)
                                                                                   cTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                                  >  1:453263/1‑150 (MQ=255)
                                                                                    tGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTta                                                                 >  2:497507/1‑150 (MQ=255)
                                                                                     ggATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGAACTGGCAAACTAACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTtat                                                                >  2:1084567/1‑150 (MQ=255)
                                                                                          aGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                                           <  1:1528705/150‑1 (MQ=255)
                                                                                            aaTGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                                         <  1:1730514/150‑1 (MQ=255)
                                                                                            aaTGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                                         <  2:1001928/150‑1 (MQ=255)
                                                                                                ttttA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGTGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATAAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                                     >  1:90628/1‑150 (MQ=255)
                                                                                                ttttA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                                     >  1:65595/1‑150 (MQ=255)
                                                                                                tattA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                                     <  1:42573/148‑1 (MQ=255)
                                                                                                      gCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgc                                                >  1:992659/1‑150 (MQ=255)
                                                                                                       cAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcg                                               >  1:215481/1‑150 (MQ=255)
                                                                                                        aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcgc                                              <  2:400408/150‑1 (MQ=255)
                                                                                                        aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcgc                                              >  2:1625246/1‑150 (MQ=255)
                                                                                                          tctcTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGctc                                            <  1:377948/150‑1 (MQ=255)
                                                                                                            tctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTCCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTg                                          <  2:109885/150‑1 (MQ=255)
                                                                                                             ctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGc                                         <  1:1591640/150‑1 (MQ=255)
                                                                                                               ttCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAAATCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGa                                       >  1:165329/1‑150 (MQ=255)
                                                                                                                 cTGTCTTGAATCCATGTCAGTGTCCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCTACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAg                                     <  2:429027/150‑1 (MQ=255)
                                                                                                                  tGTCATGAATCCATGGCAGTGACCATACTAATGTTGACTGCCATTGATGGATGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTATATAACTGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGc                                    <  1:1540818/150‑1 (MQ=255)
                                                                                                                  tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGc                                    <  1:848083/150‑1 (MQ=255)
                                                                                                                   gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAACCTCACCCCGTTTTTACTCAACCGATACCATGTAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                                   >  1:1101902/1‑150 (MQ=255)
                                                                                                                   gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                                   >  2:1700076/1‑150 (MQ=255)
                                                                                                                   gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                                   >  1:1765129/1‑150 (MQ=255)
                                                                                                                     cctccaTCAATGGCAGTCACCATTAGTATGGTCACTGCCATGGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                                 >  1:29656/6‑150 (MQ=37)
                                                                                                                     cATGAATCCATGGCAGTGACCATACTTATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATTCCTTTAAATAACAGCAACTTATACCTGTCTGTTGGCGCGCTCTGCGAAGCTGt                                 <  1:1650307/150‑1 (MQ=255)
                                                                                                                     cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                                 >  2:1450866/1‑150 (MQ=255)
                                                                                                                     cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGAAGGAGGGAGACACAGTGCCCTGGCAAACTCACACCGTTTTTCATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                                 <  2:1103258/150‑1 (MQ=255)
                                                                                                                     cATGAATCCATGGCAGTGACCATACTAATGGTGACGGCCATAGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTGTTTAATCAACCAAGACCATTAAATAACAGCAACTTATCCCTGTCTGCTGGCCCGCTCTGAGAAGCGGt                                 >  1:920351/1‑150 (MQ=255)
                                                                                                                       tGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTaa                               >  2:22638/1‑150 (MQ=255)
                                                                                                                         aaTCCATGGCAGTGGCCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcg                             <  2:1220395/150‑1 (MQ=255)
                                                                                                                           tcaaTGGCAGTCACCATTAGTATGGTCACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcgcg                           >  2:1211218/4‑150 (MQ=39)
                                                                                                                           tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGTCGCGCTCTGCGAAGCGGTAAcgcg                           >  1:1296712/1‑150 (MQ=255)
                                                                                                                           tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcgcg                           <  2:1798888/150‑1 (MQ=255)
                                                                                                                            ccATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGt                          >  2:1273615/1‑150 (MQ=255)
                                                                                                                             aaTGGCAGTCACCATTAGTATGGTCACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTg                         >  2:118435/2‑150 (MQ=255)
                                                                                                                               tGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGaa                       <  1:1632471/150‑1 (MQ=255)
                                                                                                                               tGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGaa                       <  2:607741/150‑1 (MQ=255)
                                                                                                                                 gCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAgg                     >  2:205761/1‑150 (MQ=255)
                                                                                                                                 gCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAgg                     >  1:935669/1‑150 (MQ=255)
                                                                                                                                         cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct             <  2:966492/150‑1 (MQ=255)
                                                                                                                                         cATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct             >  1:512105/1‑150 (MQ=255)
                                                                                                                                          aTACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggctg            >  2:280045/1‑150 (MQ=255)
                                                                                                                                           tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGTTgctgtctag           >  1:260095/1‑148 (MQ=255)
                                                                                                                                           tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg           >  1:1009265/1‑150 (MQ=255)
                                                                                                                                           tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg           >  2:1418235/1‑150 (MQ=255)
                                                                                                                                           tACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGctggctgg           <  2:1283654/150‑1 (MQ=255)
                                                                                                                                            aCTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTggg          <  1:83128/150‑1 (MQ=255)
                                                                                                                                              tAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTTTCTGAGGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGAt        <  2:354765/150‑1 (MQ=255)
                                                                                                                                               aaTGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATa       >  1:1466084/1‑150 (MQ=255)
                                                                                                                                                 tGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGc     <  2:1588306/150‑1 (MQ=255)
                                                                                                                                                    tGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAt  <  1:636821/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTA‑GCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGAT  >  NC_000913/4414101‑4414397

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: ‑