Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 869,018 G→A G500E (GGG→GAG)  gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,0180GA98.7% 285.0 / ‑4.5 74G500E (GGG→GAG) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (41/32);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

TGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGG  >  NC_000913/868874‑869163
                                                                                                                                                |                                                                                                                                                    
tGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGc                                                                                                                                                    >  1:289319/1‑147 (MQ=255)
   tACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTa                                                                                                                                              <  2:438045/150‑1 (MQ=255)
    aCCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTAt                                                                                                                                             <  1:602069/150‑1 (MQ=255)
      cTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCa                                                                                                                                           >  2:328624/1‑150 (MQ=255)
           cATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACttgt                                                                                                                                     >  2:753352/1‑151 (MQ=255)
              gCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACttgttg                                                                                                                                   <  2:641318/150‑1 (MQ=255)
              gCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGc                                                                                                                                  >  1:121342/1‑151 (MQ=255)
              gCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGc                                                                                                                                  <  1:112588/151‑1 (MQ=255)
              gCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGc                                                                                                                                  <  1:133715/151‑1 (MQ=255)
               ccTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGc                                                                                                                                  >  1:88249/1‑150 (MQ=255)
                 tGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCg                                                                                                                               >  2:713173/1‑151 (MQ=255)
                  ggCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCg                                                                                                                               >  1:6291/1‑150 (MQ=255)
                     gTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGAtct                                                                                                                           <  2:386351/151‑1 (MQ=255)
                          ccGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGc                                                                                                                      <  1:77682/151‑1 (MQ=255)
                            gCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGt                                                                                                                    <  2:121508/151‑1 (MQ=255)
                                gAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGAtt                                                                                                                >  1:529639/1‑151 (MQ=255)
                                gAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGAtt                                                                                                                >  1:1540/1‑151 (MQ=255)
                                gAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGAtt                                                                                                                >  1:48355/1‑151 (MQ=255)
                                 aGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGc                                                                                                                                  >  1:514196/1‑132 (MQ=255)
                                 aGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGc                                                                                                                                  <  2:514172/132‑1 (MQ=255)
                                 aGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGAttt                                                                                                               <  1:352611/151‑1 (MQ=255)
                                       ccGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCAt                                                                                                         >  2:349529/1‑151 (MQ=255)
                                       ccGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCAt                                                                                                         >  2:260465/1‑151 (MQ=255)
                                       ccGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCAt                                                                                                         >  2:189343/1‑151 (MQ=255)
                                        cggcggTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCAtt                                                                                                        <  2:48665/151‑1 (MQ=255)
                                           cggTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCg                                                                                                     <  1:656190/151‑1 (MQ=255)
                                               cagcgccagcgcATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTAt                                                                                                  >  2:168613/1‑150 (MQ=255)
                                                agcgccagcgcATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCt                                                                                                >  1:518054/1‑151 (MQ=255)
                                                 gcgccagcgcATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTg                                                                                               <  2:140527/151‑1 (MQ=255)
                                                         gcATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTAtctc                                                                                       >  1:278094/1‑151 (MQ=255)
                                                         gcATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTAtctc                                                                                       <  1:581498/151‑1 (MQ=255)
                                                          cATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTcc                                                                                      >  2:654677/1‑151 (MQ=255)
                                                               gCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCa                                                                                    <  2:658370/148‑1 (MQ=255)
                                                               gCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCAc                                                                                   <  1:747252/149‑1 (MQ=255)
                                                                 aTTGCTCGCGCGTTGGCATTGAATCCAAAAGTGACCATTGCCGACGAAGCCGTTTCGGGGCTGGTTGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGCCTTTGCGTATCTGTTTTTCTcc                                                                                      >  1:558880/1‑144 (MQ=255)
                                                                  ttGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGatat                                                                              >  1:233654/1‑151 (MQ=255)
                                                                         cgcgTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCggcg                                                                       >  2:445878/1‑149 (MQ=255)
                                                                          gcgTTGGCATTGAATCCAAAAGTGATCAGTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTAGTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGCGATTTCGGGATTGCGTATCTGTTTATCTCCCACGATATGGCgggcg                                                                      >  1:36946/1‑148 (MQ=255)
                                                                            gTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCggtgg                                                                      >  1:174715/1‑149 (MQ=255)
                                                                                 cATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCgg                                                                         >  2:294104/1‑141 (MQ=255)
                                                                                    tGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTGGAGAgg                                                              >  2:135185/1‑149 (MQ=255)
                                                                                    tGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGa                                                             <  1:580214/150‑1 (MQ=255)
                                                                                      aaTCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGAt                                                            <  2:70161/149‑1 (MQ=255)
                                                                                           aaaaGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTc                                                       >  1:700488/1‑149 (MQ=255)
                                                                                            aaaGTGATCATTGCCGACGAATCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGGGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTagag                                                                 <  2:648696/138‑1 (MQ=255)
                                                                                            aaaGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCgggaggtg                                                                   >  1:648735/1‑130 (MQ=255)
                                                                                            aaaGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTc                                                       >  1:615837/1‑148 (MQ=255)
                                                                                            aaaGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGGAGAGCggggtaattatc                                                    >  2:586993/1‑141 (MQ=255)
                                                                                              aGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCgt                                                  <  1:572131/151‑1 (MQ=255)
                                                                                                tGATCATTGCCGACGAAGCCGTTTCGGAGCTGGATGTTTCTATTCGCGGGCAGATTATCAACGTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCggttggtgagg                                                                >  2:671506/1‑127 (MQ=255)
                                                                                                    cATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCg                                            <  1:587029/151‑1 (MQ=255)
                                                                                                    cATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGTGTGGCg                                            >  1:591677/1‑151 (MQ=255)
                                                                                                     gttGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGGGGTGCGGGGTAGTCATCGcgtg                                               >  2:32932/2‑147 (MQ=255)
                                                                                                     aTTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCgg                                           <  1:314796/151‑1 (MQ=255)
                                                                                                     aTTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCgg                                           <  2:485867/151‑1 (MQ=255)
                                                                                                     aTTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCgg                                           >  1:633842/1‑151 (MQ=255)
                                                                                                     aTTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCgg                                           <  2:265824/151‑1 (MQ=255)
                                                                                                       tGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGt                                          <  2:234011/150‑1 (MQ=255)
                                                                                                          cgacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGAGGGTTTTTTATCGTGTGGCgg                                           >  2:645426/1‑146 (MQ=255)
                                                                                                           gacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGt                                     <  2:679391/151‑1 (MQ=255)
                                                                                                           gacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGATTAGTCATCGTGTGGCgggg                                         >  1:167251/1‑145 (MQ=255)
                                                                                                               aaGCCGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTAtct                                <  2:408922/151‑1 (MQ=255)
                                                                                                                  ccGTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGGGGGGTTGTCTTCGTGTGGCgggggtggt                                    >  2:158717/1‑138 (MQ=255)
                                                                                                                    gTTTCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGc                         <  2:169888/151‑1 (MQ=255)
                                                                                                                       tCGGCGCTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCa                        <  1:625209/149‑1 (MQ=255)
                                                                                                                          gcgcTGGATGTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGGGGTGTATTTCgggggaa                      >  1:537575/1‑144 (MQ=255)
                                                                                                                                   gTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATTGGGTGGGGGGGGTGGATTTCGGGGGAAATTTTGGAAAt            >  1:94829/1‑150 (MQ=255)
                                                                                                                                   gTTTCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGGGTGGATTAGTAATCGTGTGGCGGGGGTGGATTTTCGGGAAAAtggttttaa             >  1:675963/1‑145 (MQ=255)
                                                                                                                                     ttCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTgg         <  2:319237/150‑1 (MQ=255)
                                                                                                                                     ttCTATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTgg         <  1:438330/150‑1 (MQ=255)
                                                                                                                                        tATTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGGATTTCGGGGCAAATTTTTTAAATTgtccc      >  1:728997/1‑147 (MQ=255)
                                                                                                                                         aTTCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCAc    <  1:262326/151‑1 (MQ=255)
                                                                                                                                          ttCGCGAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGTTTTTCCGGGGAAATTTTTGAAATTGGTccc     >  1:344584/1‑149 (MQ=255)
                                                                                                                                             gcgAGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTTGTCATCGTTTGGCGGGGGTGGATTTCCGGGGAAAATTTTTAAATTGTTCCCCgc  >  2:606650/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGG  >  NC_000913/868874‑869163

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