Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,347,044 A→G 6.3% V623V (GTT→GTC)  rnb ← RNase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,347,0440AG 6.3% 289.3 / 7.2 95V623V (GTT→GTC) rnbRNase II
Reads supporting (aligned to +/- strand):  ref base A (38/51);  new base G (4/2);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATGCCG  >  NC_000913/1346896‑1347191
                                                                                                                                                    |                                                                                                                                                   
gTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATaa                                                                                                                                                    >  2:819853/1‑150 (MQ=255)
  tAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACg                                                                                                                                                  >  1:1421909/1‑150 (MQ=255)
      aaaTCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGATTTCGGCAATAGTGACGTCAATAACgtcag                                                                                                                                              <  1:804492/150‑1 (MQ=255)
         tCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACgtcagtca                                                                                                                                           <  1:233797/150‑1 (MQ=255)
          cTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTTATTGACgtcactatt                                                                                                                                          <  1:589727/150‑6 (MQ=255)
           tGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTTATTGACgtcactattg                                                                                                                                         <  2:1508108/150‑7 (MQ=255)
            gCCTGCATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACtt                                                                                                                                        <  1:183570/150‑1 (MQ=255)
                 gATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTaa                                                                                                                                   >  1:1616392/1‑150 (MQ=255)
                   tataACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAAc                                                                                                                                 >  1:821651/1‑150 (MQ=255)
                    ataACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAAcc                                                                                                                                <  1:144663/150‑1 (MQ=255)
                        cGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGttt                                                                                                                            >  2:1455172/1‑150 (MQ=255)
                         gAAGTTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTc                                                                                                                           <  2:1445784/150‑1 (MQ=255)
                         gAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTc                                                                                                                           <  2:1567249/150‑1 (MQ=255)
                          aaGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACGGTTTCa                                                                                                                          <  2:354871/150‑1 (MQ=255)
                          aaGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCa                                                                                                                          <  1:1522506/150‑1 (MQ=255)
                             gTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCt                                                                                                                       <  1:210067/150‑1 (MQ=255)
                                 agCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTaa                                                                                                                   <  1:1245387/150‑1 (MQ=255)
                                    ggggAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAAttt                                                                                                                <  2:718769/150‑1 (MQ=255)
                                     gggAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTg                                                                                                               <  2:121947/150‑1 (MQ=255)
                                     gggAAATAAACGGCCCATCCATGAGGAATGGGCCATGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATTCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTTTAAACCGTTTCACCTTTAAtttt                                                                                                               >  2:462478/1‑149 (MQ=255)
                                       gAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTa                                                                                                             <  2:1593523/150‑1 (MQ=255)
                                        cctttcaCGGCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTAc                                                                                                            >  1:1185619/7‑150 (MQ=255)
                                        aaaTAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCAcc                                                                                                                        >  1:1561820/1‑138 (MQ=255)
                                        aaaTAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCAcc                                                                                                                        <  2:1561820/138‑1 (MQ=255)
                                        aaaTAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTAc                                                                                                            <  2:901817/150‑1 (MQ=255)
                                          tttcaCGGCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAg                                                                                                          >  1:1618143/5‑150 (MQ=255)
                                           ttcaCGGCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGt                                                                                                         >  2:1455736/4‑150 (MQ=255)
                                           tAAACGGCCCATCCATGAGGAATGGGCAGTGAAAGGAGATTACGCGACCCGGAGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAAAAACGTCAGTCACTTTGTAAACCGTTTCACATTTAATTTGTAAAGt                                                                                                         >  2:614630/1‑150 (MQ=255)
                                            tcaCGGCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTg                                                                                                        >  2:1280590/3‑150 (MQ=255)
                                             caCGGCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGc                                                                                                       >  2:933947/2‑150 (MQ=255)
                                               cGGCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCg                                                                                                     >  1:1085870/1‑150 (MQ=255)
                                                ggCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGt                                                                                                    >  2:304217/1‑150 (MQ=255)
                                                ggCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGt                                                                                                    >  2:555090/1‑150 (MQ=255)
                                                ggCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGt                                                                                                    <  2:1616392/150‑1 (MQ=255)
                                                ggCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGt                                                                                                    <  2:821651/150‑1 (MQ=255)
                                                      tccatGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCt                                                                                              >  1:213717/1‑150 (MQ=255)
                                                      tccatGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCt                                                                                              <  1:1532333/150‑1 (MQ=255)
                                                       ccatGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTg                                                                                             <  2:985034/150‑1 (MQ=255)
                                                        catGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTgg                                                                                            <  2:1356624/150‑1 (MQ=255)
                                                         atGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGc                                                                                           >  1:615410/1‑150 (MQ=255)
                                                         atGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGc                                                                                           <  1:801068/150‑1 (MQ=255)
                                                         atGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGc                                                                                           <  2:746941/150‑1 (MQ=255)
                                                          tGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCt                                                                                          <  2:101338/150‑1 (MQ=255)
                                                           gAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTg                                                                                         >  1:505382/1‑150 (MQ=255)
                                                            aGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCTTTTTCCTGGATGc                                                                                        <  2:29109/150‑1 (MQ=255)
                                                             ggAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCa                                                                                       >  2:437444/1‑150 (MQ=255)
                                                              gAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCaa                                                                                      <  1:1368787/150‑1 (MQ=255)
                                                               aaTGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCaaa                                                                                     <  2:362715/150‑1 (MQ=255)
                                                                  gggCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCa                                                                                  >  1:1343576/1‑150 (MQ=255)
                                                                   ggCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAg                                                                                 <  1:414168/150‑1 (MQ=255)
                                                                    gCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGt                                                                                >  2:1619606/1‑150 (MQ=255)
                                                                     ccGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGtt                                                                               <  1:1226511/150‑1 (MQ=255)
                                                                       gTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCa                                                                             <  1:819853/150‑1 (MQ=255)
                                                                            aGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATcgcg                                                                        <  1:164221/150‑1 (MQ=255)
                                                                             ggAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCgcgc                                                                       >  2:1626287/1‑150 (MQ=255)
                                                                              gagaTTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGca                                                                      <  2:615410/150‑1 (MQ=255)
                                                                                gaTTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGAAGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCATTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGcaca                                                                    >  2:925886/1‑150 (MQ=255)
                                                                                   tACGCGACCGGGCGCGCAATAATGGTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAGCCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAGCCAGTTCATCGCGCACAGCg                                                                 >  2:1056722/1‑150 (MQ=255)
                                                                                        gACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTaa                                                            <  2:356448/150‑1 (MQ=255)
                                                                                         aCCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAg                                                           <  1:597585/150‑1 (MQ=255)
                                                                                         aCCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAg                                                           <  1:1455736/150‑1 (MQ=255)
                                                                                          ccGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGa                                                          >  2:1259872/1‑150 (MQ=255)
                                                                                               cgcgcAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGt                                                     >  1:423397/1‑150 (MQ=255)
                                                                                                 cgcAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGc                                                   <  2:183961/150‑1 (MQ=255)
                                                                                                     ataatGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGa                                               <  1:555090/150‑1 (MQ=255)
                                                                                                     ataatGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGa                                               <  2:1185619/150‑1 (MQ=255)
                                                                                                       aatGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAAt                                             <  1:1251709/150‑1 (MQ=255)
                                                                                                           ctgcGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAg                                         >  1:839027/1‑150 (MQ=255)
                                                                                                             gcGGGTTTCGATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCTTTTCACCTTTAATTTGTTCAGTGCCGTTTTACTTGCTGCAAACCAGTTCATCGCGCACATCGTGTAAGAAAGTTGACGGAATAAATGc                                       <  2:976541/150‑1 (MQ=255)
                                                                                                             gcGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGc                                       <  2:1272723/150‑1 (MQ=255)
                                                                                                                ggTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGAt                                    <  2:1343576/150‑1 (MQ=255)
                                                                                                                     ccATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAgcgc                               >  1:387298/1‑150 (MQ=255)
                                                                                                                      cATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGcc                              >  2:1493763/1‑150 (MQ=255)
                                                                                                                       aTGCGGACTTCTGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTAGAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCg                             <  1:315702/150‑1 (MQ=255)
                                                                                                                         gCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGtt                           <  1:355045/150‑1 (MQ=255)
                                                                                                                          cGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTa                          >  1:860181/1‑150 (MQ=255)
                                                                                                                             aCTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGTGATAGCGCCGTTATCa                       <  2:1421909/150‑1 (MQ=255)
                                                                                                                               ttCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAAc                     <  1:1026379/150‑1 (MQ=255)
                                                                                                                                tCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAAcc                    >  1:208044/1‑150 (MQ=255)
                                                                                                                                tCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAAcc                    <  2:213717/150‑1 (MQ=255)
                                                                                                                                 cGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCa                   >  1:1102603/1‑150 (MQ=255)
                                                                                                                                 cGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTAAAGTGCCGTTTTCCTGGCTCCAAACCAGTTCATAGCGCACAGCGTCTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCa                   >  2:841284/1‑150 (MQ=255)
                                                                                                                                     gtgactGACGTTATTGACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACg               >  1:1577433/6‑150 (MQ=255)
                                                                                                                                      aTAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGa              <  2:1306185/150‑1 (MQ=255)
                                                                                                                                       gactGACGTTATTGACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGaa             >  2:595666/4‑150 (MQ=255)
                                                                                                                                       gactGACGTTATTGACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGaa             >  1:414012/4‑150 (MQ=255)
                                                                                                                                       gactGACGTTATTGACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGaa             >  1:1321115/4‑150 (MQ=255)
                                                                                                                                        aGTGATGTCAATAACGTCAGTGACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGGTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAATGCGATAGCGCCGTTATCAACCAAACGAac            <  1:462478/150‑1 (MQ=255)
                                                                                                                                         gTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAaca           >  2:362954/1‑150 (MQ=255)
                                                                                                                                           gACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACg         >  2:562991/1‑150 (MQ=255)
                                                                                                                                              gTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCAt      >  1:312293/1‑150 (MQ=255)
                                                                                                                                              gTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCAt      <  1:437444/150‑1 (MQ=255)
                                                                                                                                               tCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATg     <  2:234277/150‑1 (MQ=255)
                                                                                                                                                cAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAACCCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATgc    >  1:899262/1‑150 (MQ=255)
                                                                                                                                                cAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATgc    >  1:786030/1‑150 (MQ=255)
                                                                                                                                                  aTAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATgccg  <  1:724015/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATGCCG  >  NC_000913/1346896‑1347191

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