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→ATGopgEOPG 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

GATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAG  >  NC_000913/852075‑852368
                                                                                                                                                     |                                                                                                                                                 
gATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCAt                                                                                                                                                 >  2:376694/1‑151 (MQ=255)
   aGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTa                                                                                                                                               <  1:631712/150‑1 (MQ=255)
   aGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTAc                                                                                                                                              >  1:276044/1‑151 (MQ=255)
    agcAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTAc                                                                                                                                              <  2:556345/148‑1 (MQ=255)
       aGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGcc                                                                                                                                           >  1:1109170/1‑150 (MQ=255)
       aGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGcc                                                                                                                                           >  1:1109192/1‑150 (MQ=255)
       aGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCg                                                                                                                                          >  2:451623/1‑151 (MQ=255)
        gTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATGCCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGa                                                                                                                                         <  1:134308/151‑1 (MQ=255)
                    ttAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAaaata                                                                                                                              >  1:407931/1‑150 (MQ=255)
                        aTACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAataaataa                                                                                                                         <  2:888259/151‑1 (MQ=255)
                          aCGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGc                                                                                                                       >  1:99925/1‑151 (MQ=255)
                            gCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCg                                                                                                                     <  2:469674/151‑1 (MQ=255)
                             ccACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCg                                                                                                                     >  2:660462/1‑150 (MQ=255)
                             ccACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCg                                                                                                                     >  2:527525/1‑150 (MQ=255)
                                ccccAGGCCAGTTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGAcc                                                                                                                  <  1:187569/150‑1 (MQ=255)
                                ccccAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGAcc                                                                                                                  <  1:1072/150‑1 (MQ=255)
                                   caggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGt                                                                                                              <  1:552126/151‑1 (MQ=255)
                                    aggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGt                                                                                                              >  1:380725/1‑150 (MQ=255)
                                    aggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTa                                                                                                             <  1:533273/151‑1 (MQ=255)
                                     ggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGt                                                                                                              >  2:346318/1‑149 (MQ=255)
                                     ggccaggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAg                                                                                                            <  2:503944/151‑1 (MQ=255)
                                         aggTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGt                                                                                                        >  2:7126/1‑151 (MQ=255)
                                             tttGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTAcc                                                                                                     <  2:441082/150‑1 (MQ=255)
                                              ttGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCa                                                                                                    <  1:156502/150‑1 (MQ=255)
                                               tGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCAc                                                                                                   <  1:704540/150‑1 (MQ=255)
                                               tGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACg                                                                                                  >  1:84538/1‑151 (MQ=255)
                                               tGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACg                                                                                                  >  1:57157/1‑151 (MQ=255)
                                                gCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACgg                                                                                                 <  1:330657/151‑1 (MQ=255)
                                                   cGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGaa                                                                                               >  1:341519/1‑150 (MQ=255)
                                                          aacaaAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTc                                                                                       <  1:831585/151‑1 (MQ=255)
                                                             aaaaaaTGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAgg                                                                                    >  2:1035686/1‑151 (MQ=255)
                                                               aaaaTGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGAt                                                                                  <  2:730039/151‑1 (MQ=255)
                                                                     tCGCAACACTGAATACCAGACACAAGCTGTCCAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCCGCCTCATGCCAATTCCGCCTGCTGCAAAAATAAAAAAGCCGACCTGGCAGCTGTACCACGaaaaaacctcaagaatttcg                                                                              >  1:1034276/1‑136 (MQ=255)
                                                                       gCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                          >  2:801137/1‑151 (MQ=255)
                                                                         cacaCTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                          <  2:43642/149‑1 (MQ=255)
                                                                         cacaCTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCg                                                                          >  1:43642/1‑149 (MQ=255)
                                                                            acTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCa                                                                                                    <  2:552420/120‑1 (MQ=255)
                                                                            acTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCa                                                                                                    >  1:552417/1‑120 (MQ=255)
                                                                              tGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTc                                                                   >  1:754886/1‑151 (MQ=255)
                                                                                aaTACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCa                                                                 <  1:170807/151‑1 (MQ=255)
                                                                                 aTACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCat                                                                >  1:954770/1‑151 (MQ=255)
                                                                                   aCCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCatat                                                              <  2:18633/151‑1 (MQ=255)
                                                                                     cAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtg                                                            <  1:542049/151‑1 (MQ=255)
                                                                                        acacAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtggtg                                                         >  1:168280/1‑151 (MQ=255)
                                                                                           caGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAg                                                      <  2:886120/151‑1 (MQ=255)
                                                                                            aGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGt                                                     >  1:981484/1‑151 (MQ=255)
                                                                                                 gTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCa                                                  <  1:387582/149‑1 (MQ=255)
                                                                                                       ggTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGc                                          >  1:282606/1‑151 (MQ=255)
                                                                                                        gTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCa                                         <  1:41551/151‑1 (MQ=255)
                                                                                                         tATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCa                                         >  1:825108/1‑150 (MQ=255)
                                                                                                             ccAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGt                                     >  1:832085/1‑150 (MQ=255)
                                                                                                               aGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTAtt                                  >  2:1017131/1‑151 (MQ=255)
                                                                                                                gcgcgcTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTAtt                                  >  2:1092187/1‑150 (MQ=255)
                                                                                                                 cgcgcTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTg                                 <  1:79750/150‑1 (MQ=255)
                                                                                                                   cgcTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGaaa                              <  2:247898/151‑1 (MQ=255)
                                                                                                                    gcTTTTTCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAAt                             <  2:291850/151‑1 (MQ=255)
                                                                                                                      tttttCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAAtt                            >  2:386911/1‑150 (MQ=255)
                                                                                                                          tCATTCTCTTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                       >  1:808553/1‑151 (MQ=255)
                                                                                                                                 ctTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGt                 >  1:506862/1‑150 (MQ=255)
                                                                                                                                 ctTTTTTTCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGt                 >  1:403657/1‑150 (MQ=255)
                                                                                                                                   ttttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGt                 <  1:895668/148‑1 (MQ=255)
                                                                                                                                    tttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             >  1:27643/1‑151 (MQ=255)
                                                                                                                                    tttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             >  1:330279/1‑151 (MQ=255)
                                                                                                                                     ttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             <  1:147786/150‑1 (MQ=255)
                                                                                                                                     ttttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGc             <  2:5831/150‑1 (MQ=255)
                                                                                                                                      tttCACCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGcc            <  1:293658/150‑1 (MQ=255)
                                                                                                                                       ttCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCt           <  2:1012444/150‑1 (MQ=255)
                                                                                                                                        tCGCCTGATTGCCATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTg          >  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)
                                                                                                                                                     |                                                                                                                                                 
GATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAG  >  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: