Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 6.8% intergenic (+33/‑50) ygeI → / → pbl uncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,8890TA 6.8% 239.5 / 5.5 87intergenic (+33/‑50)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (38/43);  new base A (3/3);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGG  >  NC_000913/2993740‑2994036
                                                                                                                                                     |                                                                                                                                                   
gCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAGTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                                                                                                                                                     >  2:1254744/1‑150 (MQ=255)
    cTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCATTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                                                                                                                                                 <  1:177492/150‑1 (MQ=255)
       aTACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                                                                              <  1:376551/150‑1 (MQ=255)
          cAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt                                                                                                                                           <  2:185185/150‑1 (MQ=255)
            gATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt                                                                                                                                         >  1:741213/1‑150 (MQ=255)
             aTCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAta                                                                                                                                        >  2:477168/1‑149 (MQ=255)
              tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                                                                       >  1:567041/1‑150 (MQ=255)
               aaaaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc                                                                                                                                      <  1:347590/149‑1 (MQ=255)
                aaaTACTAGCTTTTATAGTTTCAACTTGGATGGCGCCTTTTCAAAATGATTTGTTTTCAGATGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctcn                                                                                                                                     <  2:329196/150‑2 (MQ=255)
                aaaTACAAGCATTTATAGTTTCAACATTGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCg                                                                                                                                     <  1:738060/150‑1 (MQ=255)
                aaaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCg                                                                                                                                     <  1:1104256/150‑1 (MQ=255)
                 aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                                                                    <  1:699240/150‑1 (MQ=255)
                  aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt                                                                                                                                   <  1:300776/150‑1 (MQ=255)
                     cAAGCATTTATAGGTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                                                                <  1:946967/150‑4 (MQ=255)
                     aaaGCATTTATAGTTTCAAAATTGATGTCTACTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGATATACTTTTATAATAGTATCTCGTTCCt                                                                                                                                <  1:226497/149‑1 (MQ=255)
                       aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                                                                              >  1:300926/1‑150 (MQ=255)
                         cATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                                                                            >  2:980829/1‑150 (MQ=255)
                         cATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                                                                            <  2:300694/150‑1 (MQ=255)
                          aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                                                                           >  1:844718/1‑150 (MQ=255)
                             tataGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                                                                        >  2:1025545/1‑150 (MQ=255)
                               taGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAAt                                                                                                                      <  1:259450/150‑1 (MQ=255)
                                 gTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                                                                    >  2:282478/1‑150 (MQ=255)
                                        aatggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacaca                                                                                                             <  2:263996/149‑1 (MQ=255)
                                           ggatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                                                                          >  2:268267/1‑150 (MQ=255)
                                            gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                                                                         <  1:320946/150‑1 (MQ=255)
                                              tggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaa                                                                                                       >  1:528503/1‑150 (MQ=255)
                                                  gcCTTATCAAAATGATATGTATTCAGAAGATAATCCTATCACACCTTACTATAAAATTGAGGGGTTGACGATTGCCCAAATAAATCATAACGAGATACTTTTATAATAGTAGCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                                                                   >  2:688646/1‑150 (MQ=255)
                                                     ttttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgc                                                                                                >  1:829630/1‑150 (MQ=255)
                                                     ttttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgc                                                                                                <  1:857419/150‑1 (MQ=255)
                                                      tttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcg                                                                                               >  2:587848/1‑150 (MQ=255)
                                                        tCAAAATGATATGTATTCAGAATATAAAACTATCTCACCTTACTATAAAATTTAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                                                                             <  1:706599/150‑1 (MQ=255)
                                                        tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                                                                             <  2:522932/150‑1 (MQ=255)
                                                        tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACATTACTATAAAAATGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCGCACACAGGAAAAAATATTCGCt                                                                                             >  2:1200242/1‑150 (MQ=255)
                                                         cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTc                                                                                            >  1:1080808/1‑150 (MQ=255)
                                                         cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCGTACTATAAAATTGAGTGGTAGACGATTGCCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCATCACACATGAAAAAATATGCtctc                                                                                            >  1:286142/1‑150 (MQ=255)
                                                         cAAAATGATATGAATTCAGAAGATAATCCAATCTCACCTTACTATAAAATTGAGTGGTACACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATATCGTTCCTTTAATTCAATCCTCACACAAGAAAAAATATGCGCTc                                                                                            >  2:274528/1‑150 (MQ=255)
                                                           aaaTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAc                                                                                          >  2:882978/1‑150 (MQ=255)
                                                                 tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                                                                    <  1:606327/150‑1 (MQ=255)
                                                                 tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                                                                    <  2:935325/150‑1 (MQ=255)
                                                                 tatGTATTCAGAAGATAATCCTATCTCAACTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                                                                    <  1:634446/150‑1 (MQ=255)
                                                                    ttATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                                                                 <  2:380844/149‑1 (MQ=255)
                                                                       ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                                                                              <  1:980829/150‑1 (MQ=255)
                                                                          agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                                                                           >  1:857892/1‑150 (MQ=255)
                                                                           gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                                                                          >  1:522037/1‑150 (MQ=255)
                                                                            aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAACTCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                                                                         <  2:33234/150‑1 (MQ=255)
                                                                            aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                                                                         <  1:134423/150‑1 (MQ=255)
                                                                                     cTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc                                                                >  2:48018/1‑150 (MQ=255)
                                                                                        tctcACCTTACTATTAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGTTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                             <  1:468439/150‑1 (MQ=255)
                                                                                        tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                             <  1:1098713/150‑1 (MQ=255)
                                                                                        tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg                                                             <  1:320637/150‑1 (MQ=255)
                                                                                         ctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt                                                            >  1:465617/1‑150 (MQ=255)
                                                                                           cACCTTACTATAAAATTGTGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTAAACCATCTTCACGttt                                                          <  1:716789/150‑1 (MQ=255)
                                                                                           cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt                                                          <  2:896659/150‑1 (MQ=255)
                                                                                             acTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg                                                        <  1:859643/149‑1 (MQ=255)
                                                                                              cTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTCACCATCTTCACGTTTTGt                                                       >  2:1080691/1‑150 (MQ=255)
                                                                                                  cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat                                                   >  1:1251115/1‑150 (MQ=255)
                                                                                                  cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGAGCTCCCCAAAAGTTAAATTCTTATCCATCTTCCCGTTTTGtattt                                                   >  1:1000984/1‑148 (MQ=255)
                                                                                                   tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt                                                  >  2:243227/1‑150 (MQ=255)
                                                                                                    ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCATTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACAATCTTAAAGTTTTGTAtattt                                                 >  2:353625/1‑150 (MQ=255)
                                                                                                       aaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatttatt                                              >  2:1184392/1‑150 (MQ=255)
                                                                                                         aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTAc                                            >  1:269951/1‑150 (MQ=255)
                                                                                                            gAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACaaa                                         >  1:747817/1‑150 (MQ=255)
                                                                                                                ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAAtt                                     >  1:627649/1‑149 (MQ=255)
                                                                                                                 gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGt                                    <  1:544415/150‑1 (MQ=255)
                                                                                                                     aCGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCttt                                <  2:287433/150‑1 (MQ=255)
                                                                                                                      cGATTGACCAAATAAATCATAACGAGATACGTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACAGGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTCTCCATCTTCACGTTTTGTATATTTAGTACAAAAATGGCTGTc                               >  1:940336/1‑150 (MQ=255)
                                                                                                                      cGATTGACCAAATAAATCATAACGAGATAATTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCCCACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTAAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTc                               >  2:752955/1‑150 (MQ=255)
                                                                                                                         ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGcc                            <  1:983883/150‑1 (MQ=255)
                                                                                                                          tGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGATCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCa                           >  2:840199/1‑150 (MQ=255)
                                                                                                                           gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAg                          <  2:1199746/150‑1 (MQ=255)
                                                                                                                           gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCACACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGAATATTTATTACAAAAATGTCTTTCCCCAg                          >  2:723635/1‑150 (MQ=255)
                                                                                                                                aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCCATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAAct                     >  1:655947/1‑150 (MQ=255)
                                                                                                                                  taaatCATAACGAGATACATTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatttatt                                              >  1:443877/1‑123 (MQ=255)
                                                                                                                                  taaatCATAACGAGATACATTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatttatt                                              <  2:443877/123‑1 (MQ=255)
                                                                                                                                  taaatAATAAAGAGATAATTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACAAATGAAAAAATATGCGCTAACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATTTCTTTCGCCAGCAActct                   <  2:818189/150‑1 (MQ=255)
                                                                                                                                   aaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCtt                  <  1:700158/150‑1 (MQ=255)
                                                                                                                                     atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                <  1:1261548/150‑1 (MQ=255)
                                                                                                                                     atCATAACGAGATACTTTTATAATAGTATCTCGTTCATTTCATTCAATCCTCACACATGAAAAAATATGCTCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGc                <  1:558170/150‑1 (MQ=255)
                                                                                                                                      aaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGct               >  1:707465/5‑150 (MQ=255)
                                                                                                                                       aggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATTTCTTTCGCCAGCAACTCTTGcta              >  2:373480/4‑149 (MQ=255)
                                                                                                                                        ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCACTCCTCACACATGAAATAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGctct             >  1:1172112/3‑150 (MQ=255)
                                                                                                                                            cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATg         <  2:69196/150‑1 (MQ=255)
                                                                                                                                            cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATg         >  1:686968/1‑150 (MQ=255)
                                                                                                                                               taTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAg      <  2:663797/149‑1 (MQ=255)
                                                                                                                                                  aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAg   <  2:300926/150‑1 (MQ=255)
                                                                                                                                                  aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAg   <  1:730625/150‑1 (MQ=255)
                                                                                                                                                  aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAg   <  1:896433/150‑1 (MQ=255)
                                                                                                                                                   cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAgg  <  2:986420/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGG  >  NC_000913/2993740‑2994036

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: