Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,521,346 A→T 5.6% intergenic (+18/+247) lysV → / ← xapR tRNA‑Lys/transcriptional activator of xapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,521,3460AT 5.6% 235.5 / 4.3 90intergenic (+18/+247)lysV/xapRtRNA‑Lys/transcriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base A (34/50);  new base T (3/2);  total (37/53)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

GGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/2521211‑2521493
                                                                                                                                       |                                                                                                                                                   
ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                                                                                                       >  1:197266/1‑150 (MQ=255)
ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                                                                                                       <  2:539421/150‑1 (MQ=255)
ggaggGGGGCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                                                                                                       <  1:1259478/150‑1 (MQ=255)
 gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAc                                                                                                                                      <  1:1067801/150‑1 (MQ=255)
 gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAc                                                                                                                                      <  1:1233966/150‑1 (MQ=255)
  tgggggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCTTTTTTAca                                                                                                                                     <  2:619658/149‑2 (MQ=255)
      ggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCTTTTTTAcattgg                                                                                                                                 <  1:1042416/150‑6 (MQ=255)
        tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGc                                                                                                                               <  1:672577/150‑1 (MQ=255)
        tCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGc                                                                                                                               <  1:170709/150‑1 (MQ=255)
         cggcggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCg                                                                                                                              <  2:733922/150‑1 (MQ=255)
         cggcggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCg                                                                                                                              <  2:1421211/150‑1 (MQ=255)
            cggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                                                                                                           >  2:1221772/1‑150 (MQ=255)
            cggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTATTACTATCTGCGATa                                                                                                                           >  2:1012539/1‑150 (MQ=255)
                tctaTCACGTCATAACCCACCTCTTCGGGTCGTTTTCTCTGTTGGTAGAGCAGTTGACTTTTAATCAATTGTTCGCAGGTTCGAATCATGCACGAACCACCAATGTAAAAAAGCTCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCa                                                                                                                       <  1:1535268/147‑1 (MQ=14)
                tCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCa                                                                                                                       <  1:522779/150‑1 (MQ=255)
                 cGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaa                                                                                                                      >  2:1174526/1‑150 (MQ=255)
                  gATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                                                                                                     >  1:963011/1‑150 (MQ=255)
                  gATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                                                                                                     >  1:1664804/1‑150 (MQ=255)
                   aTCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAg                                                                                                                    >  1:178584/1‑150 (MQ=255)
                     cccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAt                                                                                                                  <  1:587894/150‑1 (MQ=255)
                       cGTCATCACCCACCTCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTc                                                                                                                <  2:66435/150‑1 (MQ=255)
                         tcatcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGa                                                                                                              <  2:575511/150‑1 (MQ=255)
                           atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTATAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAc                                                                                                            <  2:1450925/150‑1 (MQ=255)
                            tcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAcc                                                                                                           <  1:1458286/150‑1 (MQ=255)
                               cccaccACTTGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgc                                                                                                        <  1:809437/150‑1 (MQ=255)
                               cccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgc                                                                                                        <  1:786474/150‑1 (MQ=40)
                                 caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGATTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCTATACTCAAAGATTCGAACCTgcag                                                                                                      <  1:840836/150‑1 (MQ=255)
                                 caccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                      <  1:214213/150‑1 (MQ=39)
                                  accacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagc                                                                                                     <  2:1409655/150‑1 (MQ=38)
                                   ccacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagca                                                                                                    <  2:50245/150‑1 (MQ=38)
                                      cTTCGGGGCGTTACCGCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTGCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt                                                                                                 <  1:1264031/150‑1 (MQ=11)
                                        tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                               <  1:127502/150‑1 (MQ=33)
                                          gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  2:790602/150‑1 (MQ=32)
                                          gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                             <  1:469848/150‑1 (MQ=32)
                                              cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGGttgagttg                                                                                         >  2:53674/1‑150 (MQ=17)
                                              cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg                                                                                         <  2:536267/150‑1 (MQ=32)
                                                     tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                  >  1:1222029/1‑150 (MQ=32)
                                                      cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAATCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                                                                 >  2:354584/1‑150 (MQ=17)
                                                        gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                               >  1:1873962/1‑150 (MQ=32)
                                                          tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                             <  2:1663062/150‑1 (MQ=32)
                                                           ggTAGAGCAGTTGACTTTTAATCAATTGCTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                            <  1:194143/150‑1 (MQ=17)
                                                            gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                           >  2:1799498/1‑150 (MQ=32)
                                                                cgCAGTTGACTTTTACTCACTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTCAAATCGCGCCCTCAAGGCGCTTTTTTACTAGCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTGGAGCGCAGCGAAACAACgg                                                                       <  1:46557/149‑1 (MQ=11)
                                                                 gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                      >  1:174082/1‑150 (MQ=32)
                                                                   aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGCCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAGTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                    >  1:1146227/1‑150 (MQ=17)
                                                                   aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                    >  1:1332627/1‑150 (MQ=32)
                                                                        aCTTTTAATCAATTGGTCGCAGGTTCGAAGCCTGCACGACCCCCCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                                                               >  1:1099861/1‑150 (MQ=17)
                                                                         cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                                                              <  1:773434/150‑1 (MQ=32)
                                                                          ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                                                             >  1:1402517/1‑150 (MQ=32)
                                                                            ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                           <  1:1039987/150‑1 (MQ=32)
                                                                            ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                           >  1:689757/1‑150 (MQ=32)
                                                                             tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                                          <  2:595746/150‑1 (MQ=32)
                                                                             tAATCAATTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGAAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCTATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                                          <  2:1702427/150‑1 (MQ=12)
                                                                                 cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAGGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACg                                                      >  1:1128011/1‑150 (MQ=17)
                                                                                   aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTACAACAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCACAGATTCGCACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc                                                    <  1:292069/150‑1 (MQ=255)
                                                                                    ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc                                                   <  2:526768/150‑1 (MQ=32)
                                                                                      ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAAGGCCCg                                                 >  2:691657/1‑150 (MQ=17)
                                                                                       gTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGa                                                >  2:576593/1‑150 (MQ=32)
                                                                                       gTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCGGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCAGa                                                >  1:946079/1‑150 (MQ=17)
                                                                                           cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                            <  2:19306/150‑1 (MQ=32)
                                                                                           cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAATCAAAGATTCGAACATGCAGCAGGTTTGAGTTGAGAGCAGCGAAACAACGGAGCCGCTAGCGGCGACGGCCCGAAggg                                            >  2:1107784/1‑150 (MQ=255)
                                                                                            aGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc                                           <  2:98491/150‑1 (MQ=32)
                                                                                               ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                        <  2:1222029/150‑1 (MQ=32)
                                                                                                 cGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcg                                      <  2:1861515/150‑1 (MQ=32)
                                                                                                     tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                  >  1:671642/1‑150 (MQ=32)
                                                                                                     cccTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                  <  1:506334/149‑1 (MQ=17)
                                                                                                     accTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                  <  2:342492/149‑1 (MQ=17)
                                                                                                     accTGCACGACCCACCAATGAAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAAACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                  <  1:743872/149‑1 (MQ=17)
                                                                                                      ccTGCACGACCCACCAATTTAAAAACGCGCCATAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcg                                 <  2:531797/150‑1 (MQ=255)
                                                                                                        cgCACGACCCACCAATGTAAAAACGCGCCCTACAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAg                               <  1:81144/149‑1 (MQ=17)
                                                                                                               cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCt                        >  2:266949/1‑150 (MQ=32)
                                                                                                                ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACCCAAAGATTCGAACCTGCAGCAGGTTGGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGCCGCCGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTg                       >  1:758101/1‑150 (MQ=255)
                                                                                                                  cagatagTAAAAAAGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCATGCa                     >  2:1738882/7‑150 (MQ=11)
                                                                                                                  accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     <  2:854651/150‑1 (MQ=32)
                                                                                                                   aaaaTTTAAAAAAGCGCATTAAAGGCGTTTTTTTATTATTTGCGATACTAAAATTTTGTACGCTGAAGTTTGTTTTTGTTGTGTGATGGTAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc                    <  2:1760964/148‑1 (MQ=2)
                                                                                                                      aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAc                 >  2:1096953/1‑150 (MQ=32)
                                                                                                                        gTAAAAAAGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc               >  1:1417776/1‑150 (MQ=12)
                                                                                                                        gTAAAAAAGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcccc              >  1:456818/1‑149 (MQ=11)
                                                                                                                         tAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGTTACTCAACGATTCGACCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca              <  1:688001/150‑1 (MQ=12)
                                                                                                                         tAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGccccc              >  2:2009747/1‑149 (MQ=17)
                                                                                                                          aaaaaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac             >  1:1313404/1‑150 (MQ=32)
                                                                                                                            aaacGCGCCCTAACGGCGATTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           <  2:197266/150‑1 (MQ=11)
                                                                                                                             aaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGTTACTCATAGATTCGAACCTGAAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGTCGCGAAGGGCGAGCGATGCGAGTCATCCTGCCCGGCCCGCCaa          >  1:130483/1‑150 (MQ=11)
                                                                                                                             aaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa          >  1:1174071/1‑150 (MQ=32)
                                                                                                                              aaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt         <  1:614686/150‑1 (MQ=32)
                                                                                                                               cgcgcCCTAAAGGCGCTTTTTTACTATCTGCGATAATCAAAGATTCGAACCTGCAGCAGGTTTTAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCTAAGCGAGTCATCCTGCACGACCCACCAATg        <  1:162503/149‑1 (MQ=12)
                                                                                                                                 cgcCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTa      >  2:162138/1‑150 (MQ=32)
                                                                                                                                   cccTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaa    >  2:1027760/1‑150 (MQ=32)
                                                                                                                                   cccTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaa    >  1:498837/1‑150 (MQ=32)
                                                                                                                                     cTAAACGCGCTTTTTTACTATCTGCGATACTCAAAGCTTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa  <  2:560761/150‑1 (MQ=17)
                                                                                                                                       |                                                                                                                                                   
GGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/2521211‑2521493

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