Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 852,223 T→C I125M (ATA→ATG)  opgE ← OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913852,2230TC98.7% 287.4 / NA 74I125M (ATA→ATG) opgEOPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (36/37);  minor base A (1/0);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAG  >  NC_000913/852075‑852368
                                                                                                                                                     |                                                                                                                                                 
gATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCAt                                                                                                                                                 >  2:376694/1‑151 (MQ=255)
   aGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTa                                                                                                                                               <  1:631712/150‑1 (MQ=255)
   aGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTAc                                                                                                                                              >  1:276044/1‑151 (MQ=255)
    agcAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTAc                                                                                                                                              <  2:556345/148‑1 (MQ=255)
       aGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGcc                                                                                                                                           >  1:1109170/1‑150 (MQ=255)
       aGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGcc                                                                                                                                           >  1:1109192/1‑150 (MQ=255)
       aGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCg                                                                                                                                          >  2:451623/1‑151 (MQ=255)
        gTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATGCCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGa                                                                                                                                         <  1:134308/151‑1 (MQ=255)
                    ttAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAaaata                                                                                                                              >  1:407931/1‑150 (MQ=255)
                        aTACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAataaataa                                                                                                                         <  2:888259/151‑1 (MQ=255)
                          aCGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGc                                                                                                                       >  1:99925/1‑151 (MQ=255)
                            gCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCg                                                                                                                     <  2:469674/151‑1 (MQ=255)
                             ccACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCg                                                                                                                     >  2:660462/1‑150 (MQ=255)
                             ccACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCg                                                                                                                     >  2:527525/1‑150 (MQ=255)
                                ccccAGGCCAGTTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGAcc                                                                                                                  <  1:187569/150‑1 (MQ=255)
                                ccccAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGAcc                                                                                                                  <  1:1072/150‑1 (MQ=255)
                                   caggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGt                                                                                                              <  1:552126/151‑1 (MQ=255)
                                    aggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGt                                                                                                              >  1:380725/1‑150 (MQ=255)
                                    aggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTa                                                                                                             <  1:533273/151‑1 (MQ=255)
                                     ggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGt                                                                                                              >  2:346318/1‑149 (MQ=255)
                                     ggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAg                                                                                                            <  2:503944/151‑1 (MQ=255)
                                         aggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGt                                                                                                        >  2:7126/1‑151 (MQ=255)
                                             tttGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTAcc                                                                                                     <  2:441082/150‑1 (MQ=255)
                                              ttGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCa                                                                                                    <  1:156502/150‑1 (MQ=255)
                                               tGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCAc                                                                                                   <  1:704540/150‑1 (MQ=255)
                                               tGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACg                                                                                                  >  1:84538/1‑151 (MQ=255)
                                               tGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACg                                                                                                  >  1:57157/1‑151 (MQ=255)
                                                gCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACgg                                                                                                 <  1:330657/151‑1 (MQ=255)
                                                   cGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGaa                                                                                               >  1:341519/1‑150 (MQ=255)
                                                          aacaaAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTc                                                                                       <  1:831585/151‑1 (MQ=255)
                                                             aaaaaaTGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAgg                                                                                    >  2:1035686/1‑151 (MQ=255)
                                                               aaaaTGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGAt                                                                                  <  2:730039/151‑1 (MQ=255)
                                                                     tCGCAACACTGAATACCAGACACAAGCTGTCCAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCCGCCTCATGCCAATTCCGCCTGCTGCAAAAATAAAAAAGCCGACCTGGCAGCTGTACCACGaaaaaacctcaagaatttcg                                                                              >  1:1034276/1‑136 (MQ=255)
                                                                       gCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                          >  2:801137/1‑151 (MQ=255)
                                                                         cacaCTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                          <  2:43642/149‑1 (MQ=255)
                                                                         cacaCTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                          >  1:43642/1‑149 (MQ=255)
                                                                            acTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCa                                                                                                    <  2:552420/120‑1 (MQ=255)
                                                                            acTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCa                                                                                                    >  1:552417/1‑120 (MQ=255)
                                                                              tGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTc                                                                   >  1:754886/1‑151 (MQ=255)
                                                                                aaTACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCa                                                                 <  1:170807/151‑1 (MQ=255)
                                                                                 aTACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCat                                                                >  1:954770/1‑151 (MQ=255)
                                                                                   aCCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCatat                                                              <  2:18633/151‑1 (MQ=255)
                                                                                     cAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtg                                                            <  1:542049/151‑1 (MQ=255)
                                                                                        acacAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtggtg                                                         >  1:168280/1‑151 (MQ=255)
                                                                                           caGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAg                                                      <  2:886120/151‑1 (MQ=255)
                                                                                            aGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGt                                                     >  1:981484/1‑151 (MQ=255)
                                                                                                 gTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCa                                                  <  1:387582/149‑1 (MQ=255)
                                                                                                       ggTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGc                                          >  1:282606/1‑151 (MQ=255)
                                                                                                        gTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCa                                         <  1:41551/151‑1 (MQ=255)
                                                                                                         tATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCa                                         >  1:825108/1‑150 (MQ=255)
                                                                                                             ccAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGt                                     >  1:832085/1‑150 (MQ=255)
                                                                                                               aGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTAtt                                  >  2:1017131/1‑151 (MQ=255)
                                                                                                                gcgcgcTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTAtt                                  >  2:1092187/1‑150 (MQ=255)
                                                                                                                 cgcgcTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTg                                 <  1:79750/150‑1 (MQ=255)
                                                                                                                   cgcTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGaaa                              <  2:247898/151‑1 (MQ=255)
                                                                                                                    gcTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAAt                             <  2:291850/151‑1 (MQ=255)
                                                                                                                      tttttCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAAtt                            >  2:386911/1‑150 (MQ=255)
                                                                                                                          tCATTCTCTTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                       >  1:808553/1‑151 (MQ=255)
                                                                                                                                 ctTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGt                 >  1:506862/1‑150 (MQ=255)
                                                                                                                                 ctTTTTTT‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGt                 >  1:403657/1‑150 (MQ=255)
                                                                                                                                   tttttt‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGt                 <  1:895668/148‑1 (MQ=255)
                                                                                                                                    ttttt‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             >  1:27643/1‑151 (MQ=255)
                                                                                                                                    ttttt‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             >  1:330279/1‑151 (MQ=255)
                                                                                                                                     tttt‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             <  1:147786/150‑1 (MQ=255)
                                                                                                                                     tttt‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             <  2:5831/150‑1 (MQ=255)
                                                                                                                                      ttt‑CACCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGcc            <  1:293658/150‑1 (MQ=255)
                                                                                                                                       tt‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCt           <  2:1012444/150‑1 (MQ=255)
                                                                                                                                        t‑CGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTg          >  1:395226/1‑150 (MQ=255)
                                                                                                                                          cGCCTGTGTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTg          >  1:597810/1‑149 (MQ=255)
                                                                                                                                          cGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAg        <  2:742376/151‑1 (MQ=255)
                                                                                                                                           gCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGc       <  1:98202/151‑1 (MQ=255)
                                                                                                                                               gATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCaaaa   <  2:657931/151‑1 (MQ=255)
                                                                                                                                               gATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCaa     <  1:645606/149‑1 (MQ=255)
                                                                                                                                                aTTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAg  <  2:950806/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAG  >  NC_000913/852075‑852368

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑