Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,571 T→C 5.3% intergenic (+4/+155) ykgQ → / ← rclC pseudogene, putative dehydrogenase/reactive chlorine species (RCS) stress resistance inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,5710TC 5.3% 438.5 / 10.8 150intergenic (+4/+155)ykgQ/rclCpseudogene, putative dehydrogenase/reactive chlorine species (RCS) stress resistance inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (76/66);  new base C (3/5);  total (79/71)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-02

Only 100 of 154 total aligned reads displayed.

GATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCG  >  NC_000913/317422‑317720
                                                                                                                                                     |                                                                                                                                                     
gATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATAtt                                                                                                                                                       <  1:532041/150‑1 (MQ=255)
  tGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                                                                                     >  1:299094/1‑150 (MQ=255)
   tcccATCAAAAATCGATTTCTCACTGACTGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAAt                                                                                                                                                    <  2:581787/149‑1 (MQ=255)
    cccATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                                                                                   >  1:807508/1‑150 (MQ=255)
    cccATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                                                                                   >  1:653447/1‑150 (MQ=255)
     ccATCAAAAATCGATTTCTCACTGACAGATAAAGTAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                                                                                  <  2:163123/150‑1 (MQ=255)
       aTCAAAAATGGATTTCTGACTTACAGATAAAGAATTGGGACTTTTTAGTTCGATTGATTTTGGATTTTTGGGGGAGTCGGTAAAACATTTTAATACTGATTTTGTTCGTGGTTGGCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                                                                                <  1:1190401/150‑1 (MQ=255)
       aTCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATTTTAATCATg                                                                                                                                                >  2:120703/1‑150 (MQ=255)
       aTCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                                                                                >  2:1123429/1‑150 (MQ=255)
         cAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGTTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTAGGTAAAATACATGATTAATATCAATCATGTAtt                                                                                                                                            <  1:906967/150‑1 (MQ=38)
         cAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                                                                              <  1:408366/150‑1 (MQ=255)
         cAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                                                                              >  1:765657/1‑150 (MQ=255)
          aaaaaTCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAt                                                                                                                                             <  1:1695829/150‑1 (MQ=255)
           aaaaTCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtt                                                                                                                                            <  1:689466/150‑1 (MQ=255)
            aaaTCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                                                                           >  2:192765/1‑150 (MQ=255)
             aaTCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTTATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                                                                          <  2:95027/150‑1 (MQ=255)
              aTCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTAGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                                                                         >  1:102287/1‑150 (MQ=255)
                 gATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCt                                                                                                                                      <  2:337652/150‑1 (MQ=255)
                 gATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGATCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCt                                                                                                                                      <  1:767326/150‑1 (MQ=255)
                  aTTTCTCTCTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGGCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacag                                                                                                                                     <  2:1073677/150‑6 (MQ=39)
                     tctcACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAga                                                                                                                                  >  1:813569/1‑150 (MQ=255)
                        cACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacta                                                                                                                               >  1:53926/1‑148 (MQ=255)
                        cACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacga                                                                                                                               >  2:468347/1‑150 (MQ=255)
                               aGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                                                                        <  1:1558800/150‑1 (MQ=255)
                               aGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                                                                        >  2:1267732/1‑150 (MQ=255)
                                gATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                                                                       >  1:823022/1‑150 (MQ=255)
                                gATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                                                                       >  2:307823/1‑150 (MQ=255)
                                   aaaGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                                                                                    <  2:1538397/150‑1 (MQ=255)
                                   aaaGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                                                                                    >  2:152617/1‑150 (MQ=255)
                                    aaGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                                                                   <  2:223723/150‑1 (MQ=255)
                                    aaGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                                                                   >  1:387243/1‑150 (MQ=255)
                                    aaGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                                                                   <  2:881463/150‑1 (MQ=255)
                                    aaGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                                                                   <  2:1258686/150‑1 (MQ=255)
                                      gAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGATGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTaa                                                                                                                 >  1:1305461/1‑150 (MQ=255)
                                       aaaTGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                                                                                <  1:508364/150‑1 (MQ=255)
                                       aaaTGGCAAAAATTAGTTCGCTTGATTTGGGTTATGTTTGGGAGTCGTTAAAACATTTTAATCATTAATTTTTTCTTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAAGCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                                                                                <  1:121393/150‑1 (MQ=255)
                                        aaTGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGAGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                                                                               >  2:1158751/1‑150 (MQ=255)
                                         aTGGCACAAATTAGTTAGATTTATTTGGTTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTTTTCGTTTTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATTTATTTTACCTAAGACGACAAAAACCTTTTAACAg                                                                                                              <  2:1490575/150‑1 (MQ=255)
                                           gtcacaAATTAGTTCGCTTGATTTGGGTTATTTTGGTGATTCGGTAAATCATTTTAATCCTGAATTTGTTCGTTGTTGTCTTGCTGTAATAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                                                                            <  1:1148207/148‑1 (MQ=255)
                                           ggCTCAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                                                                            <  1:1104118/150‑1 (MQ=255)
                                            gCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGa                                                                                                           >  1:1239279/1‑150 (MQ=255)
                                            gCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGa                                                                                                           >  2:1098321/1‑150 (MQ=255)
                                                       ttCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGAt                                                                                                <  1:468347/150‑1 (MQ=255)
                                                         cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCt                                                                                              <  2:18786/150‑1 (MQ=255)
                                                         cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCt                                                                                              <  2:1379392/150‑1 (MQ=255)
                                                          gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCtt                                                                                             >  2:1136340/1‑150 (MQ=255)
                                                           cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTc                                                                                            >  2:1364704/1‑150 (MQ=255)
                                                           cTTGATTTGGGGTATGTTGGGGAGTCGGCAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGAGGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTc                                                                                            <  2:704795/150‑1 (MQ=255)
                                                               aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                                                                        >  1:22922/1‑150 (MQ=255)
                                                                   gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                                                                                    <  1:1267732/150‑1 (MQ=255)
                                                                     gTTATTTTTGGGAGTCGGTAATACATTTTAATCCTGTTTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTTATCATGAATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                                                                                  >  2:1097800/1‑150 (MQ=255)
                                                                      ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGAATGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                                                                                 <  2:224885/150‑1 (MQ=255)
                                                                         tGTTGGGGAGTCTGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtct                                                                              <  1:692127/150‑1 (MQ=255)
                                                                            tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTATTCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                                                                           >  1:151064/1‑150 (MQ=255)
                                                                             ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                                                                            >  2:39129/1‑84 (MQ=255)
                                                                               ggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATATTCACTATGTCATATATCCGGTa                                                                        >  2:598028/1‑150 (MQ=255)
                                                                                      gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGc                                                                 <  1:1696746/150‑1 (MQ=255)
                                                                                        aaaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGt                                                               >  2:336969/1‑150 (MQ=255)
                                                                                         aaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc                                                              >  1:568559/1‑150 (MQ=255)
                                                                                             aTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGc                                                          >  2:165786/1‑150 (MQ=255)
                                                                                             aTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGc                                                          <  2:161022/150‑1 (MQ=255)
                                                                                                  aaTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcat                                                     >  1:1498780/1‑150 (MQ=255)
                                                                                                    tCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatca                                                   >  1:591678/1‑150 (MQ=255)
                                                                                                     ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatcat                                                  >  2:1552161/1‑150 (MQ=255)
                                                                                                      cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtcatcatc                                                 <  2:1595063/150‑1 (MQ=255)
                                                                                                       tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCAGGTAATATGACGTAATGCTtcatcataa                                                >  1:164609/1‑148 (MQ=255)
                                                                                                        gAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                                                                                    >  1:1584811/1‑97 (MQ=255)
                                                                                                         aaTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGa                                              >  2:1675655/1‑150 (MQ=255)
                                                                                                             tGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgg                                          >  1:1612224/1‑150 (MQ=255)
                                                                                                             tGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTAATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgg                                          >  1:1019437/1‑150 (MQ=255)
                                                                                                              gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAggg                                         <  2:332689/150‑1 (MQ=255)
                                                                                                                tCGTGGTTGTCTTGCTGTAAAAATACATGAATGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCTGTAATATGGCGTCAGGCTTCATCATCAGAAAggggg                                       <  2:424616/150‑1 (MQ=255)
                                                                                                                   tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCt                                    >  2:686646/1‑150 (MQ=255)
                                                                                                                     gTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTgc                                  >  1:1446752/1‑149 (MQ=255)
                                                                                                                      ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGAc                                 <  2:813569/150‑1 (MQ=255)
                                                                                                                         gcTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGcc                              <  1:1566659/149‑1 (MQ=255)
                                                                                                                           ttGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTg                            <  2:1690205/150‑1 (MQ=255)
                                                                                                                           ttGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTg                            <  2:1052845/150‑1 (MQ=255)
                                                                                                                            tGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCATCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGa                           <  2:845773/150‑1 (MQ=255)
                                                                                                                                ttAAAAATACATGATTTATATTAATCATGTATTTTACCTAATACGACAAAAACCTTTTAACAGGGAGTGAAATTGAGCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGACAGGGGGCCTGACGCCTGaaaaa                       <  2:231866/149‑1 (MQ=255)
                                                                                                                                gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATATCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGTGGGCCTGACGCATGaaaaa                       >  2:1313649/1‑150 (MQ=255)
                                                                                                                                 tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAg                      <  2:105239/150‑1 (MQ=255)
                                                                                                                                 tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTAATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAg                      <  2:1362101/150‑1 (MQ=255)
                                                                                                                                  aaaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGt                     <  2:1430545/150‑1 (MQ=255)
                                                                                                                                   aaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTg                    >  2:1460462/1‑150 (MQ=255)
                                                                                                                                    aaaTACATGATTGATATTAATCATGTATTTTAACTAAGACGACAAAAACCTTTTAACAGGGAGTGATAATGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGa                   <  2:253060/150‑1 (MQ=255)
                                                                                                                                      aTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACATGGAGTGATATTGATCTTCAATCTGTCATATATCCGGTAATATGGCGTCAGGCTTAATAAACAGAAAGGGGGCCTGACGCCTGAAAAAGTTaac                 >  1:1095051/1‑150 (MQ=255)
                                                                                                                                        aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaacaa               >  1:359119/1‑150 (MQ=255)
                                                                                                                                        aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaacaa               >  2:877059/1‑150 (MQ=255)
                                                                                                                                         cATGATTAATATCAATCATGTATTTTTACAGCAAGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCAGGTAATATGGCGTCATGCTTCATCATCAGAAAGGGGGCCTGACTCCTGAAAAAGTGAacaaca             >  2:728304/1‑150 (MQ=37)
                                                                                                                                          aTGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAacaaca             <  2:793156/150‑1 (MQ=255)
                                                                                                                                            gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACCGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacaga           >  1:439289/1‑150 (MQ=255)
                                                                                                                                             attgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagac          <  2:1305461/150‑1 (MQ=255)
                                                                                                                                             attgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagac          <  2:1702109/150‑1 (MQ=255)
                                                                                                                                               tgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagacag        <  1:1141512/150‑1 (MQ=255)
                                                                                                                                                gatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAgt       >  1:1359800/1‑150 (MQ=255)
                                                                                                                                                 atatTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAgtg      >  1:1481621/1‑150 (MQ=255)
                                                                                                                                                 atatTAAACATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAgtg      >  1:1486393/1‑150 (MQ=255)
                                                                                                                                                     tAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCg  <  1:307823/150‑1 (MQ=255)
                                                                                                                                                     tAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCg  <  1:1527322/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCG  >  NC_000913/317422‑317720

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: