Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,414,234 G→C 7.0% intergenic (+43/‑41) 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,2340GC 7.0% 504.0 / 22.9 157intergenic (+43/‑41)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/66);  new base C (6/5);  total (86/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

Only 100 of 158 total aligned reads displayed.

GTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGC  >  NC_000913/4414085‑4414381
                                                                                                                                                      |                                                                                                                                                   
  aGCTAGATTGTCGATGATGAAGCGTGTGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGAc                                                                                                                                                   <  2:1196502/150‑1 (MQ=255)
   gCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGAcc                                                                                                                                                  >  1:1693449/1‑150 (MQ=255)
            tCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATg                                                                                                                                         >  1:1283654/1‑150 (MQ=255)
             cGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCAGGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATgg                                                                                                                                        <  1:919813/150‑1 (MQ=255)
             cGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTGTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATgg                                                                                                                                        >  2:1295415/1‑150 (MQ=255)
             cGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTATCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATgg                                                                                                                                        <  1:718593/150‑1 (MQ=255)
               aggatgtATCGTGTGGTATGGGGAGCCGTGAAGATAACCAACTGATGACCAAAGACACCTCACGGTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTg                                                                                                                                      <  1:1575267/148‑1 (MQ=255)
                    gaAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACCCCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTGTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGc                                                                                                                                 >  2:329277/1‑150 (MQ=255)
                      aGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCCCGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCa                                                                                                                               >  2:1681881/1‑150 (MQ=255)
                      aGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCAATACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACAATACTAATGGTGACTGCCa                                                                                                                               >  1:429027/1‑150 (MQ=255)
                       gCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTCACCATTAGTATGGTCACTGCCAt                                                                                                                              <  2:1074044/150‑1 (MQ=255)
                        cGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCAtt                                                                                                                             <  1:874194/150‑1 (MQ=255)
                          tGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGa                                                                                                                           >  2:1476851/1‑150 (MQ=255)
                           gCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACAATACTAATGGTGACTGCCATTGAt                                                                                                                          >  1:1148339/1‑150 (MQ=255)
                            cGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATg                                                                                                                         <  1:419242/150‑1 (MQ=255)
                             ggAATGGGGATCCGTGAAGATAACCCACTGATCACCAACGCCACCTCCGGTTTCGTTCCGCAGTTCATGGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATgg                                                                                                                        >  1:145403/1‑150 (MQ=255)
                                  ggggATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTTGCAATCTCTTTCTGTCATGAATCAATGGCAGTGACCATACTAATGGTGACTGCCATTGATggaggga                                                                                                                   >  2:319693/1‑150 (MQ=255)
                                     gATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGac                                                                                                                <  1:580719/150‑1 (MQ=255)
                                         cGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAg                                                                                                            >  1:1188610/1‑150 (MQ=255)
                                            gAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGc                                                                                                         >  2:618796/1‑150 (MQ=255)
                                             aaGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCa                                                                                                        >  2:106134/1‑150 (MQ=255)
                                                ctAACCCACTGCTCCCCAAAGCCACCTCACCTTTCGGTCCGCAGTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACAACTGGGGACTGCCATTGATGGAGGGAGACACAGTGCACTg                                                                                                     <  2:1694906/149‑1 (MQ=255)
                                                 tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggagggag                                                                                                                  <  2:460779/136‑1 (MQ=255)
                                                 tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTgg                                                                                                    <  2:691285/150‑1 (MQ=255)
                                                 tAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTgg                                                                                                    <  2:1588152/150‑1 (MQ=255)
                                                 tAAAACACTGATCACCAAAGACACCTCACGTTTAGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTCCCATTGATGGAGGTAGACACAGTGCACTgg                                                                                                    >  2:380609/1‑150 (MQ=255)
                                                  aaCACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGc                                                                                                   >  2:1425791/1‑150 (MQ=255)
                                                    cacacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGAACATACTAATGGTGACTGCCATTTATTGATGGAGACACAGTGCAATGTCaa                                                                                                 >  1:1212743/1‑150 (MQ=255)
                                                           aTCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTcaca                                                                                          <  2:1014743/150‑1 (MQ=255)
                                                              aCCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTTGATAAGAATGTTTTAGCAATCACTTTCTGTCATGAATCCATGGCAGTGACCATCCTAATGGTGAATGACATTGATGTATGGATACACAGTGCACTGGCAAACTcacaact                                                                                       >  1:1125768/1‑147 (MQ=255)
                                                               ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGt                                                                                      >  1:1798888/1‑150 (MQ=255)
                                                               ccAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGt                                                                                      >  2:1765884/1‑150 (MQ=255)
                                                                cAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGtt                                                                                     >  2:814970/1‑150 (MQ=255)
                                                                cAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGtt                                                                                     >  2:959375/1‑150 (MQ=255)
                                                                   ttacacCTCACGTTTCGTTCGGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCTTGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGttttt                                                                                  <  2:149841/148‑1 (MQ=255)
                                                                    gACACCTCACGTTTCGGTCCGCAGTACATTGATTGATAAGAAGTGATACGCAATATCCTTCTGTCATGAATCCATTGCAGTGACCATACTAATGGTGAATTGCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTa                                                                                 <  1:84388/150‑1 (MQ=255)
                                                                       acCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATc                                                                              <  1:1635479/150‑1 (MQ=255)
                                                                        ccTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                             >  1:646960/1‑150 (MQ=255)
                                                                        ccTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCACTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                             >  1:1026069/1‑150 (MQ=255)
                                                                          tCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAc                                                                           <  2:256731/150‑1 (MQ=255)
                                                                           cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                          <  2:920746/150‑1 (MQ=255)
                                                                           cACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAAcc                                                                          <  2:565728/150‑1 (MQ=255)
                                                                               tttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                      <  2:126223/150‑1 (MQ=255)
                                                                               tttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCtata                                                                      <  1:284988/150‑1 (MQ=255)
                                                                                ttCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAc                                                                     >  1:563663/1‑150 (MQ=255)
                                                                                 tCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAcc                                                                    <  1:1389528/150‑1 (MQ=255)
                                                                                  cGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCa                                                                   >  2:855851/1‑150 (MQ=255)
                                                                                      ccGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTaa                                                               >  1:1236823/1‑150 (MQ=255)
                                                                                        gCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                             >  2:1707084/1‑150 (MQ=255)
                                                                                        gCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAAt                                                             >  1:35118/1‑150 (MQ=255)
                                                                                         cATTACATTGCTGGATAAGAATGTTTTAGAAATCTCTATCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATa                                                            <  1:1604567/150‑1 (MQ=255)
                                                                                           ttACATTGCTTGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAAc                                                          <  1:861517/150‑1 (MQ=255)
                                                                                           ttACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAAc                                                          <  1:763844/150‑1 (MQ=255)
                                                                                            tACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACa                                                         >  2:532757/1‑150 (MQ=255)
                                                                                             aCATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAg                                                        <  2:1199886/150‑1 (MQ=255)
                                                                                             aCATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAg                                                        <  2:1566377/150‑1 (MQ=255)
                                                                                             aCATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAg                                                        >  1:194464/1‑150 (MQ=255)
                                                                                              cATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGc                                                       >  2:1635040/1‑150 (MQ=255)
                                                                                                   cTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                  <  1:558524/150‑1 (MQ=255)
                                                                                                   cTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACtt                                                  >  1:453263/1‑150 (MQ=255)
                                                                                                    tGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTta                                                 >  2:497507/1‑150 (MQ=255)
                                                                                                          aGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                           <  1:106134/150‑1 (MQ=255)
                                                                                                          aGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                           <  1:1528705/150‑1 (MQ=255)
                                                                                                          aGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTg                                           >  1:472723/1‑150 (MQ=255)
                                                                                                            aaTGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                         <  1:1730514/150‑1 (MQ=255)
                                                                                                            aaTGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                         <  2:1001928/150‑1 (MQ=255)
                                                                                                            aaTGTTTTAGCAAACTCTTTCTGTCATGAATACATGGCAGTGACCATACTAATTGTGACTGACATTTTTGGATGGATACACATTTCACTTGCAAACTAACACCGTTTTTTATCAACCTATACCATTAAATAACAGCAACTTATACCTGTc                                         <  1:1122670/150‑1 (MQ=255)
                                                                                                               gTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGa                                      >  1:1272845/1‑150 (MQ=255)
                                                                                                                ttttAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGAt                                     >  1:65595/1‑150 (MQ=255)
                                                                                                                  ttAGCATTCTCTTTATGTCATGTATCCATGGCAGTGACCATACTAATGGTGACTGCCTTTGATGGAGGGAGACACAGTGCACTTGCAAACTCACACCGTTTTTAATCAACCTATACCATTATATAACAGCAACTTATACCTGTCTGATgg                                   <  2:1090528/150‑1 (MQ=255)
                                                                                                                  ttAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATgg                                   <  1:259698/150‑1 (MQ=255)
                                                                                                                      gCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGAGACAGTGCACTGGCAACCTCACCCCGTTTTTAATCAGCCTATTCCATGAAAGAAGAGCAAGTTATACCTGTGTGATGgcgc                                >  2:307854/1‑150 (MQ=255)
                                                                                                                       cAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcg                               >  1:215481/1‑150 (MQ=255)
                                                                                                                        aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcgc                              >  2:1625246/1‑150 (MQ=255)
                                                                                                                        aaTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgcgcgc                              <  2:400408/150‑1 (MQ=255)
                                                                                                                         aTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGct                             >  1:1816057/1‑150 (MQ=255)
                                                                                                                          tctcTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGctc                            <  1:377948/150‑1 (MQ=255)
                                                                                                                            tctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTg                          <  1:532757/150‑1 (MQ=255)
                                                                                                                            tctTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTCCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTg                          <  2:109885/150‑1 (MQ=255)
                                                                                                                             ctTTCTGTCATGAATCCATTGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGc                         <  2:626738/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)
                                                                                                                                  tGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACCCCGTTTTTAATCACCCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGc                    >  2:43904/1‑150 (MQ=255)
                                                                                                                                   gTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAACCTCACCCCGTTTTTACTCAACCGATACCATGTAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                   >  1:1101902/1‑150 (MQ=255)
                                                                                                                                     cctccaTCAATGGCAGTCACCATTAGTATGGTCACTGCCATGGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                 >  1:29656/6‑150 (MQ=37)
                                                                                                                                     cATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGt                 >  2:1450866/1‑150 (MQ=255)
                                                                                                                                     cATGAATCCATGGCAGTGACCATACTAATGGTGACGGCCATAGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTGTTTAATCAACCAAGACCATTAAATAACAGCAACTTATCCCTGTCTGCTGGCCCGCTCTGAGAAGCGGt                 >  1:920351/1‑150 (MQ=255)
                                                                                                                                      ctccaTCAATGGCAGTCACCATTAGTATGGTCACTGCCATGGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTa                >  1:1748928/5‑150 (MQ=37)
                                                                                                                                       tccaTCAATGGCAGTCACCATTAGTATGGTCAATGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAAATTATACCGGTCTGATGGCGCGCTCTGCGAAGCGGTaa               >  2:1119174/4‑150 (MQ=38)
                                                                                                                                       tGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTaa               >  2:22638/1‑150 (MQ=255)
                                                                                                                                         aaTCCATGGCAGTGGCCATACTAATGGTGACTGCCATTGATGGATGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcg             <  2:1214417/150‑1 (MQ=255)
                                                                                                                                           tcaaTGGCAGTCACCATTAGTATGGTCACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcgcg           >  2:1211218/4‑150 (MQ=39)
                                                                                                                                           tCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAAcgcg           <  2:1798888/150‑1 (MQ=255)
                                                                                                                                            ccATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGt          >  2:1273615/1‑150 (MQ=255)
                                                                                                                                                ggCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAg      >  2:899184/1‑150 (MQ=255)
                                                                                                                                                 gCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAgg     >  1:935669/1‑150 (MQ=255)
                                                                                                                                                  cAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGt    >  2:569865/1‑150 (MQ=255)
                                                                                                                                                    gTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgc  >  2:1471087/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGC  >  NC_000913/4414085‑4414381

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