Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,347,040 C→T 6.3% D625N (GAC→AAC)  rnb ← RNase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,347,0400CT 6.3% 298.1 / 6.4 95D625N (GAC→AAC) rnbRNase II
Reads supporting (aligned to +/- strand):  ref base C (38/51);  new base T (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 = 6.64e-01

AGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATGC  >  NC_000913/1346891‑1347189
                                                                                                                                                     |                                                                                                                                                     
aGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTc                                                                                                                                                       <  2:25438/150‑1 (MQ=255)
   cAGTTAATAAATCTGCCTGGTTATAACGAAGGTAGTGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCGGTGTATGGAGATTACGCGAGCGGGCGCGCAATATTGCTGCGGGTTTCCTTGCGGACTTCGGCAATAGTGACGTCAAt                                                                                                                                                    <  1:1166141/150‑1 (MQ=255)
   cAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAAt                                                                                                                                                    >  1:366897/1‑150 (MQ=255)
     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:555090/1‑150 (MQ=255)
                                                     ggCCCATTCCTCATGGATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGt                                                                                                  >  2:304217/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:801068/150‑1 (MQ=255)
                                                              atGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGc                                                                                         >  1:615410/1‑150 (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                                             <  2:1185619/150‑1 (MQ=255)
                                                                                                          ataatGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGa                                             <  1:555090/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                  <  2:213717/150‑1 (MQ=255)
                                                                                                                                     tCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAAcc                  >  1:208044/1‑150 (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:1321115/4‑150 (MQ=255)
                                                                                                                                            gactGACGTTATTGACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGaa           >  1:414012/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)
                                                                                                                                                     |                                                                                                                                                     
AGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATGC  >  NC_000913/1346891‑1347189

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: