Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,254,807 C→T D554N (GAT→AAT)  plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,254,8070CT98.6% 276.3 / NA 73D554N (GAT→AAT) plsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (39/32);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGGCAGCAATATTATTG  >  NC_000913/4254661‑4254952
                                                                                                                                                  |                                                                                                                                                 
ttGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGAGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAgg                                                                                                                                               >  1:735591/1‑151 (MQ=255)
  gCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGt                                                                                                                                              <  2:623348/150‑1 (MQ=255)
  gCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGt                                                                                                                                              >  1:742927/1‑150 (MQ=255)
       cgcgGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCagt                                                                                                                                        >  2:70185/1‑151 (MQ=255)
        gcgGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCagt                                                                                                                                        <  1:741642/150‑1 (MQ=255)
         cgGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCagttg                                                                                                                                      >  1:751358/1‑151 (MQ=255)
                atgatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCt                                                                                                                               <  1:674659/151‑1 (MQ=255)
                  gatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCg                                                                                                                             >  1:4212/1‑151 (MQ=255)
                  gatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCg                                                                                                                             <  2:111599/151‑1 (MQ=255)
                  gatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCg                                                                                                                             <  2:742882/151‑1 (MQ=255)
                         tCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACt                                                                                                                      >  1:338614/1‑151 (MQ=255)
                          cGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACt                                                                                                                      <  1:555129/150‑1 (MQ=255)
                             cGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTc                                                                                                                  <  2:309366/151‑1 (MQ=255)
                               aTTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTcgc                                                                                                                >  1:164886/1‑151 (MQ=255)
                                ttGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTcgcg                                                                                                               >  2:550539/1‑151 (MQ=255)
                                  gtgtCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCgg                                                                                                              <  2:737552/150‑1 (MQ=255)
                                    gtCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtg                                                                                                           >  1:707965/1‑151 (MQ=255)
                                    gtCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgt                                                                                                            <  1:206490/150‑1 (MQ=255)
                                        ttCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtgagtg                                                                                                       >  1:200816/1‑151 (MQ=255)
                                         tctcGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGgtgagtg                                                                                                       <  2:316825/150‑1 (MQ=255)
                                           tcGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAgc                                                                                                    >  1:587699/1‑151 (MQ=255)
                                              aCTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCt                                                                                                 <  2:736199/151‑1 (MQ=255)
                                                 tCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTg                                                                                                <  1:720748/149‑1 (MQ=255)
                                                      cctGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGAAGAGTCCGTGGAGTAGTGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGAt                                                                                         >  2:530207/3‑151 (MQ=255)
                                                      cTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGAt                                                                                         >  1:299841/1‑151 (MQ=255)
                                                          ttCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGCGATGc                                                                                       >  2:441533/1‑149 (MQ=255)
                                                           tCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCa                                                                                     <  1:483316/150‑1 (MQ=255)
                                                             aTTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAg                                                                                    >  2:204248/1‑149 (MQ=255)
                                                              tttGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGt                                                                                   >  1:528278/1‑149 (MQ=255)
                                                               ttGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAgc                                                                                >  1:65269/1‑151 (MQ=255)
                                                                tGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAgc                                                                                <  2:666597/150‑1 (MQ=255)
                                                                  cAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAgcgc                                                                              <  1:504625/150‑1 (MQ=255)
                                                                  cAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAgcgc                                                                              >  1:376304/1‑150 (MQ=255)
                                                                   aGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAgcgc                                                                              >  2:155506/1‑149 (MQ=255)
                                                                    gcgcGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGt                                                                           <  1:654863/151‑1 (MQ=255)
                                                                       cgTGATCGATAAGCTCGCTGGCGCTGGTTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACa                                                                        <  1:466506/151‑1 (MQ=255)
                                                                         tGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGCTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGGAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAg                                                                       >  2:5012/1‑150 (MQ=255)
                                                                            tCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGcac                                                                    >  2:180251/1‑150 (MQ=255)
                                                                             cGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGcaca                                                                   <  2:280076/150‑1 (MQ=255)
                                                                               aTAAGCTCGCTGGCGCTGGCTGAAGGAACGGGCGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATGCAGGTAGCAGTTGAGTTGCTCGGGTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGGCAGTAGCGCAGTACAGCACAGGt                                                                >  2:246883/1‑151 (MQ=255)
                                                                                 aaGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTc                                                              >  1:368744/1‑151 (MQ=255)
                                                                                  aGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTc                                                              >  1:213846/1‑150 (MQ=255)
                                                                                    cTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATg                                                           >  1:138147/1‑151 (MQ=255)
                                                                                      cgctggcgctggCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGc                                                         >  1:77430/1‑151 (MQ=255)
                                                                                         tggcgctggcTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGtt                                                      >  1:669755/1‑151 (MQ=255)
                                                                                          ggcgctggcTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGt                                                       <  2:200526/149‑1 (MQ=255)
                                                                                            cgctggcTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGc                                                   >  1:539964/1‑151 (MQ=255)
                                                                                            cgctggcTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGc                                                   <  2:449468/151‑1 (MQ=255)
                                                                                            cgctggcTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGc                                                   >  1:409278/1‑151 (MQ=255)
                                                                                              ctggctgAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCcg                                                 >  2:254156/1‑151 (MQ=255)
                                                                                                ggctgAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCcgc                                                <  2:339017/150‑1 (MQ=255)
                                                                                                   tgAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCt                                            >  2:95429/1‑151 (MQ=255)
                                                                                                       ggAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCgt                                        >  1:401618/1‑151 (MQ=255)
                                                                                                         aaCGGTAGAGTCCGTGGAGTAGGGCACGTTGCGGATCAGATCCAGGTAGCAGTTGCGTTCCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTGCAGCACAGGTTCATTGCGTTTGCCGCGCCTGCgtt                                       <  2:320550/150‑1 (MQ=255)
                                                                                                          aCGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCgttgt                                     <  2:338600/151‑1 (MQ=255)
                                                                                                             gTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTaa                                  <  2:305987/151‑1 (MQ=255)
                                                                                                              tAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAAt                                 >  1:751325/1‑151 (MQ=255)
                                                                                                                    ccGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCt                                            <  2:151643/134‑1 (MQ=255)
                                                                                                                    ccGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCt                                            >  1:151650/1‑134 (MQ=255)
                                                                                                                     cGTGGAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTAcc                          <  2:725021/151‑1 (MQ=255)
                                                                                                                        ggAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATc                       <  2:752922/151‑1 (MQ=255)
                                                                                                                         gAGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATc                       >  1:177286/1‑150 (MQ=255)
                                                                                                                          aGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATc                       >  1:143576/1‑149 (MQ=255)
                                                                                                                          aGTAGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAg                     <  2:600193/151‑1 (MQ=255)
                                                                                                                             aGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATc                       >  1:185202/1‑146 (MQ=255)
                                                                                                                             aGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATc                       >  1:760192/1‑146 (MQ=255)
                                                                                                                             aGGGCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGAt                   <  2:30713/150‑1 (MQ=255)
                                                                                                                                gCACGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGgc               <  1:191853/151‑1 (MQ=255)
                                                                                                                                   cGTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGgcagc            >  1:688386/1‑151 (MQ=255)
                                                                                                                                    gTTGCGCATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGgcagca           <  1:622640/151‑1 (MQ=255)
                                                                                                                                       gcgcATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGGCAGCAata        <  1:574951/151‑1 (MQ=255)
                                                                                                                                        cgcATCAGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGGCAGCAatat       <  1:584838/151‑1 (MQ=255)
                                                                                                                                              aGATTCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGGCAGCAATATTATTg  >  2:407052/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAGGAACGGTAGAGTCCGTGGAGTAGGGCACGTTGCGCATCAGATCCAGGTAGCAGTTGAGTTGCTCGGTTAACTGCTCGCGGGTGAGTGAGCGCTGACGTGATGCCAGTAGCGCAGTACAGCACAGGTTCATGGCGTTTGCCGCGCCTGCGTTGTTAATGCGTACCATCAGATCGGCAGCAATATTATTG  >  NC_000913/4254661‑4254952

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