Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 888,802 (A)7→6 coding (241/1686 nt) b0847 ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913888,7960A.99.5% 904.9 / ‑4.0 198coding (247/1686 nt)b0847putative transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (100/97);  minor base G (0/1);  total (100/98)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

CAGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAATGT  >  NC_000913/888704‑888884
                                                                                            |                                                                                        
 aGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaa                                                                                         <  2:4104744/93‑1 (MQ=255)
  gCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCCCCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaa                                                                                        <  2:4690960/93‑1 (MQ=255)
  gCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaa                                                                                        <  2:2257444/93‑1 (MQ=255)
  gCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaa                                                                                        >  1:4766460/1‑93 (MQ=255)
    ttACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaa                                                                                      >  2:3245924/1‑93 (MQ=255)
     tACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaa                                                                                     <  1:5631956/93‑1 (MQ=255)
     tACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaa                                                                                     <  2:5664339/93‑1 (MQ=255)
      aCCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaag                                                                                    >  2:6172523/1‑92 (MQ=255)
      aCCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaag                                                                                    <  2:4439982/93‑2 (MQ=255)
       ccTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaagt                                                                                   >  1:5123443/1‑91 (MQ=255)
        cTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGaaaaaagtt                                                                                  <  1:2052702/93‑4 (MQ=255)
          aaCCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCg                                                                               <  2:2371287/93‑1 (MQ=255)
           aCCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCgg                                                                              >  2:2505880/1‑93 (MQ=255)
           aCCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCgg                                                                              <  1:467828/93‑1 (MQ=255)
           aCCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCgg                                                                              >  2:3431422/1‑93 (MQ=255)
           aCCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCgg                                                                              >  1:1311959/1‑93 (MQ=255)
            cccTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGt                                                                             <  2:1658590/93‑1 (MQ=255)
             ccTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTc                                                                            >  2:5415454/1‑93 (MQ=255)
             ccTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTc                                                                            >  2:708123/1‑93 (MQ=255)
             ccTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTc                                                                            >  1:2847421/1‑93 (MQ=255)
             ccTAAGGCGATCACCAGCGCACTGCCAACCATAACCAGTGCAAAAATAAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTc                                                                            >  1:3010437/1‑93 (MQ=255)
              ccaaGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTcc                                                                           >  2:3971125/3‑93 (MQ=255)
              cTAAGGCGATCACCAGCGCACTGCCAACCATCGCCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTcc                                                                           <  1:3579346/93‑1 (MQ=255)
               tAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTTGGTCCg                                                                          <  1:4955362/93‑1 (MQ=255)
               tAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCg                                                                          >  1:1392274/1‑93 (MQ=255)
                aaGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCgg                                                                         >  1:2908793/1‑93 (MQ=255)
                 aGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGc                                                                        >  2:336035/1‑93 (MQ=255)
                  ggCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCt                                                                       >  1:3486068/1‑93 (MQ=255)
                    cGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTc                                                                     <  2:4212331/93‑1 (MQ=255)
                    cGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTc                                                                     >  2:4397111/1‑93 (MQ=255)
                    cGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTc                                                                     <  1:5503125/93‑1 (MQ=255)
                    cGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTc                                                                     >  1:669442/1‑93 (MQ=255)
                     gATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCg                                                                    >  1:3425872/1‑93 (MQ=255)
                      aTCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGa                                                                   <  2:1709004/93‑1 (MQ=255)
                      aTCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGa                                                                   >  1:2077074/1‑93 (MQ=255)
                       tCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAc                                                                  >  1:3836276/1‑93 (MQ=255)
                       tCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAc                                                                  >  2:5132947/1‑93 (MQ=255)
                         aCCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAccc                                                                >  1:3892266/1‑93 (MQ=255)
                         aCCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAccc                                                                <  2:5211563/93‑1 (MQ=255)
                          ccAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGAcccc                                                               >  1:696445/1‑93 (MQ=255)
                            aGCGCACTGCCAGCCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGa                                                             <  1:3331296/93‑1 (MQ=255)
                            aGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGa                                                             >  2:556918/1‑93 (MQ=255)
                                cACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCa                                                         <  1:4303809/93‑1 (MQ=255)
                                 aCTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAg                                                        >  1:176054/1‑93 (MQ=255)
                                  cTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGa                                                       >  2:3735243/1‑93 (MQ=255)
                                  cTGCCAACCATAACCAGTGCTAACATTAGGTAATTGTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGa                                                       <  2:1276857/93‑1 (MQ=255)
                                   tGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGaa                                                      <  1:3223984/93‑1 (MQ=255)
                                     ccAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGaaaa                                                    <  2:3827775/93‑1 (MQ=255)
                                     ccAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGaaaa                                                    <  2:4977274/93‑1 (MQ=255)
                                       aaCCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATa                                                  >  1:2529742/1‑93 (MQ=255)
                                            tCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAg                                             <  1:1000025/93‑1 (MQ=255)
                                            tCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAg                                             <  1:709961/93‑1 (MQ=255)
                                            tCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAg                                             >  2:4906420/1‑93 (MQ=255)
                                              aCCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCa                                           <  1:4881735/93‑1 (MQ=255)
                                              aCCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCa                                           <  2:4564184/93‑1 (MQ=255)
                                               ccAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCAt                                          >  1:4465963/1‑93 (MQ=255)
                                                cAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATa                                         <  2:4338387/93‑1 (MQ=255)
                                                 agcgCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATaa                                        <  1:2278631/90‑1 (MQ=255)
                                                  gTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATaaa                                       <  1:1131573/93‑1 (MQ=255)
                                                  gTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATaaa                                       <  1:2330444/93‑1 (MQ=255)
                                                     cTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGcc                                    <  2:3876679/93‑1 (MQ=255)
                                                       aaCATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCaa                                  >  1:5423323/1‑93 (MQ=255)
                                                        aCATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAg                                 >  1:2641984/1‑93 (MQ=255)
                                                         cATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGa                                <  2:5828501/93‑1 (MQ=255)
                                                         cATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGa                                <  1:4810515/93‑1 (MQ=255)
                                                          aTTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGACAAGAt                               >  1:5297555/1‑93 (MQ=255)
                                                           ttAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAAGAGCCTAAAGCCAAGAta                              >  2:4013750/1‑92 (MQ=255)
                                                           ttAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGAtt                              >  2:2166242/1‑93 (MQ=255)
                                                           ttAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGAtt                              >  2:1100519/1‑93 (MQ=255)
                                                             aGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGTTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTaa                            <  2:2922485/93‑1 (MQ=255)
                                                             aGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTaa                            <  1:5952765/93‑1 (MQ=255)
                                                             aGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTaa                            >  2:4860285/1‑93 (MQ=255)
                                                             aGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTaa                            <  1:2738215/93‑1 (MQ=255)
                                                              ggTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAg                           >  2:4702820/1‑93 (MQ=255)
                                                              ggTAATTTTTCCCATCGCGAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgc                          <  1:3714367/93‑1 (MQ=255)
                                                               gTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgc                          >  1:6230051/1‑93 (MQ=255)
                                                               gTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgc                          <  2:3067935/93‑1 (MQ=255)
                                                               gTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGATCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgc                          >  1:1519855/1‑93 (MQ=255)
                                                                tAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgcg                         >  1:5921437/1‑93 (MQ=255)
                                                                tAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAgcg                         <  1:6256539/93‑1 (MQ=255)
                                                                 aaTTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGCTGCAGAAAATAAACAGCATAAAGCAAAGATTAAgcga                        >  1:877137/1‑92 (MQ=255)
                                                                  aTTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCa                       >  2:3475942/1‑93 (MQ=255)
                                                                  aTTTTTCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCa                       <  1:2404478/93‑1 (MQ=255)
                                                                    ttttCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATc                     <  2:4274633/93‑1 (MQ=255)
                                                                     tttCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCg                    >  1:3779418/1‑93 (MQ=255)
                                                                     tttCCCATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCg                    <  2:1313311/93‑1 (MQ=255)
                                                                         ccATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGgtgt                >  1:3082441/1‑93 (MQ=255)
                                                                         ccATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGgtgt                <  1:2141751/93‑1 (MQ=255)
                                                                          cATCGCGACAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGtt               >  2:6132134/1‑93 (MQ=255)
                                                                          cATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGtt               <  2:2737758/93‑1 (MQ=255)
                                                                          cATCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGtt               <  2:66054/93‑1 (MQ=255)
                                                                           aTCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTa              >  1:4059992/1‑93 (MQ=255)
                                                                            tCGCGAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTaa             <  1:4613727/93‑1 (MQ=255)
                                                                                gAAAAAAAATGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCt         <  2:4252094/93‑1 (MQ=255)
                                                                                 aaaaaaaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTg        <  1:4107241/93‑1 (MQ=255)
                                                                                  aaaaaaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGa       <  1:1694882/93‑1 (MQ=255)
                                                                                    aaaaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGaaa     <  2:1696659/93‑1 (MQ=255)
                                                                                      aaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAAtg   >  2:841623/1‑93 (MQ=255)
                                                                                      aaaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAAtg   >  2:1936719/1‑93 (MQ=255)
                                                                                       aaTGGAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAAtgt  >  2:6146913/1‑93 (MQ=255)
                                                                                            |                                                                                        
CAGCTTACCTAACCCTAAGGCGATCACCAGCGCACTGCCAACCATCACCAGTGCTAACATTAGGTAATTTTTCCCATCGCGAAAAAAAATGGAAAAAAAGTTCGGTCCGGCTTCGACCCCGACGCAGAAAATAAACAGCATAAAGCCAAGATTAAGCGCATCGGTGTTAATGCTGAAATGT  >  NC_000913/888704‑888884

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