Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,593,3200TG54.7% ‑1.3 / 100.3 64N40H (AAC→CAC) livFbranched‑chain amino acid ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (17/12);  new base G (19/16);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-03

CCTTCCGGGACAATCGCCACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCTAAATAGGCACGG  >  NC_000913/3593175‑3593463
                                                                                                                                                 |                                                                                                                                               
ccTTCCGGGACAATCGCCACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTgcgcgc                                                                                                                                            <  2:181749/151‑1 (MQ=255)
    ccGGGACAATCGCCACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAAt                                                                                                                                        >  2:777954/1‑151 (MQ=255)
      gggACAATCGCCACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCa                                                                                                                                      <  1:578648/151‑1 (MQ=255)
         aCAATCGCCACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCg                                                                                                                                   >  2:46415/1‑151 (MQ=255)
          cAATCGCCACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCg                                                                                                                                   >  2:160942/1‑150 (MQ=255)
                 cACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAAtc                                                                                                                           <  1:674920/151‑1 (MQ=255)
                  aCCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAAtct                                                                                                                          >  2:930367/1‑151 (MQ=255)
                         cgcgCATGATTTTCGCTGTCTGCCAGTCGTTAATGTCTTTATCATCAAACAAAATTCGCCGGCTGGTGGCACGCGGATCTCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGc                                                                                                                       <  1:860527/147‑1 (MQ=255)
                         cgcgCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCtg                                                                                                                   <  2:307993/151‑1 (MQ=255)
                         cgcgCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCtg                                                                                                                   >  1:8060/1‑151 (MQ=255)
                             cATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCtgattg                                                                                                               >  2:687485/1‑151 (MQ=255)
                                gATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCtgattgat                                                                                                             >  2:643475/1‑150 (MQ=255)
                                gATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCtgattgat                                                                                                             >  2:524979/1‑150 (MQ=255)
                                gATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCtgattgat                                                                                                             <  1:385382/150‑1 (MQ=255)
                                    ttCGCTGTCTGCCAGTCGGTAATGTCTTCATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGc                                                                                                         >  1:440872/1‑150 (MQ=255)
                                       gCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGc                                                                                                     >  1:75201/1‑151 (MQ=255)
                                            cTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCAcc                                                                                                >  2:883962/1‑151 (MQ=255)
                                              gCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCt                                                                                               <  1:864884/150‑1 (MQ=255)
                                                cAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCAt                                                                                            >  2:421464/1‑151 (MQ=255)
                                                cAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCAt                                                                                            >  2:607455/1‑151 (MQ=255)
                                                 aGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATg                                                                                           <  2:542739/151‑1 (MQ=255)
                                                  gTCGGTAATGTCTTTATCATCAAACCCAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGc                                                                                          <  1:589353/151‑1 (MQ=255)
                                                      gTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAgcg                                                                                      >  1:97852/1‑151 (MQ=255)
                                                       tAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAgcg                                                                                      <  2:3765/150‑1 (MQ=255)
                                                           gTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTgg                                                                                 >  2:126028/1‑151 (MQ=255)
                                                               ttATCATCAATCACAATTCGCCCGCTGGTGGCCCGCGGATCGCCGCATAACGGGCTGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGa                                                                                <  2:865662/148‑1 (MQ=255)
                                                                   catcaAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGcc                                                                         >  1:118154/1‑151 (MQ=255)
                                                                     tcaAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCg                                                                        <  1:725503/150‑1 (MQ=255)
                                                                       aaaCACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGt                                                                       >  1:35376/1‑149 (MQ=255)
                                                                        aaCACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGGGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGt                                                                       >  1:786935/1‑148 (MQ=255)
                                                                          cacaATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTgg                                                                  >  2:99224/1‑151 (MQ=255)
                                                                           acaATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTgg                                                                  <  1:758322/150‑1 (MQ=255)
                                                                           acaATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTgg                                                                  <  1:114686/150‑1 (MQ=255)
                                                                                  gCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTg                                                            <  1:759564/149‑1 (MQ=255)
                                                                                     cGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACttt                                                       <  1:95530/151‑1 (MQ=255)
                                                                                      gCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTg                                                      >  2:155213/1‑151 (MQ=255)
                                                                                       cTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGAGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTg                                                      <  2:865640/150‑1 (MQ=255)
                                                                                       cTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGt                                                     >  2:66234/1‑151 (MQ=255)
                                                                                        tggtggCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTc                                                    >  1:288579/1‑151 (MQ=255)
                                                                                                   ggATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGacaaca                                         <  2:111699/151‑1 (MQ=255)
                                                                                                   ggATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGacaaca                                         <  1:160939/151‑1 (MQ=255)
                                                                                                   ggATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGacaac                                          >  2:654898/1‑150 (MQ=255)
                                                                                                     aTCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATg                                       <  2:844130/151‑1 (MQ=255)
                                                                                                       cgccgcATAACGTGCCGAGCAAGGTGGTTTTCCTCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGAc                                     >  1:912668/1‑151 (MQ=255)
                                                                                                        gccgcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACt                                    >  1:855271/1‑151 (MQ=255)
                                                                                                         ccgcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACtt                                   >  2:702284/1‑151 (MQ=255)
                                                                                                         ccgcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACtt                                   >  1:895857/1‑151 (MQ=255)
                                                                                                         ccgcATAACGTGCCGAGCAAGGTGGTTTGCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGAc                                     <  2:756879/149‑1 (MQ=255)
                                                                                                          cgcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACttt                                  >  2:808375/1‑151 (MQ=255)
                                                                                                           gcATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACtttt                                 >  1:160557/1‑151 (MQ=255)
                                                                                                            cATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACtttt                                 <  1:539634/150‑1 (MQ=255)
                                                                                                               aaCGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTcc                              >  1:342923/1‑150 (MQ=255)
                                                                                                               aaCGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCa                             <  1:568846/151‑1 (MQ=255)
                                                                                                                    gCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTa                        <  1:412270/151‑1 (MQ=255)
                                                                                                                     ccGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTAt                       >  1:117773/1‑151 (MQ=255)
                                                                                                                      cGAGCAAGGTGGTTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATg                      <  2:647239/151‑1 (MQ=255)
                                                                                                                           aaGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCa                 <  1:26842/151‑1 (MQ=255)
                                                                                                                           aaGGTGGTTTTCCCCACCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCa                 <  1:768615/151‑1 (MQ=255)
                                                                                                                             ggtggtTTTCCCCGCCCCGTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCAc                <  1:47487/150‑1 (MQ=255)
                                                                                                                              gtggtTTTCCCCGCCCCGGGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCCCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCt              >  1:199513/1‑151 (MQ=255)
                                                                                                                                        ccgccccgTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCTAAATAGGCAc    >  1:525451/1‑151 (MQ=255)
                                                                                                                                        ccgccccgTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCTAAATAGGCAc    <  1:541441/151‑1 (MQ=255)
                                                                                                                                        ccgccccgTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCTAAATAGGCAc    >  2:606853/1‑151 (MQ=255)
                                                                                                                                          gccccgTGCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCTAAATAGGCACgg  >  1:220569/1‑151 (MQ=255)
                                                                                                                                           ccccgTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCTAAATAGGCACgg  >  2:323484/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCTTCCGGGACAATCGCCACCGCTTCGCGCATGATTTTCGCTGTCTGCCAGTCGGTAATGTCTTTATCATCAAACACAATTCGCCCGCTGGTGGCACGCGGATCGCCGCATAACGTGCCGAGCAAGGTGGTTTTCCCCGCCCCGTTCGCGCCAATCAGCGTGACAATCTCGCCCTGATTGATATGCAGGCTCACCTCATGCAGCGCCTGGATTTTGCCGTAGTGGGCGCTGACTTTGTCAAAGGACAACATGACTTTTTCCATCTTATGCCTCACCTAAATAGGCACGG  >  NC_000913/3593175‑3593463

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