Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,896 T→A 6.9% intergenic (+40/‑43) 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,8960TA 6.9% 211.4 / 7.5 88intergenic (+40/‑43)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base A (3/3);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-02

ATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAATG  >  NC_000913/2993747‑2994043
                                                                                                                                                     |                                                                                                                                                   
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:300694/150‑1 (MQ=255)
                  cATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                                                                                   >  2:980829/1‑150 (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                                                                                           <  2:935325/150‑1 (MQ=255)
                                                          tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                                                                           <  1:606327/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                >  1:686968/1‑150 (MQ=255)
                                                                                                                                     cGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATg                <  2:69196/150‑1 (MQ=255)
                                                                                                                                        taTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAg             <  2:663797/149‑1 (MQ=255)
                                                                                                                                           aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAg          <  1:730625/150‑1 (MQ=255)
                                                                                                                                           aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAg          <  2:300926/150‑1 (MQ=255)
                                                                                                                                           aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAg          <  1:896433/150‑1 (MQ=255)
                                                                                                                                            cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAgg         <  2:986420/150‑1 (MQ=255)
                                                                                                                                               ttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTACCCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAAc      >  2:1104742/1‑150 (MQ=255)
                                                                                                                                               ttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAAc      <  1:1278882/150‑1 (MQ=255)
                                                                                                                                                tataATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAGATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACa     >  1:358841/1‑150 (MQ=255)
                                                                                                                                                   aataGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAATg  >  2:671728/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATTTATTACAAAAATGTCTTTCGCCAGCAACTCTTGCTCTAATGAAGCAGGAACAATG  >  NC_000913/2993747‑2994043

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