Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,347,042 A→T 6.1% I624N (ATT→AAT)  rnb ← RNase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,347,0420AT 6.1% 310.2 / 5.7 97I624N (ATT→AAT) rnbRNase II
Reads supporting (aligned to +/- strand):  ref base A (40/51);  new base T (4/2);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATGCCG  >  NC_000913/1346894‑1347191
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                        <  2:1561820/138‑1 (MQ=255)
                                          aaaTAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCAcc                                                                                                                        >  1:1561820/1‑138 (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: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)
                                                                                                                                                    |                                                                                                                                                     
CAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCGCGCAATAATGCTGCGGGTTTCCATGCGGACTTCGGCAATAGTGACGTCAATAACGTCAGTCACTTTGTAAACCGTTTCACCTTTAATTTGTACAGTGCCGTTTTCCTGGCTGCAAACCAGTTCATCGCGCACAGCGTGTAAGAAAGGTGCCGGAATAAAGGCGATAGCGCCGTTATCAACCAAACGAACACGCATGCCG  >  NC_000913/1346894‑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: