Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 914,564 T→C intergenic (‑436/+11) ybjE ← / ← aqpZ putative surface protein/transmembrane water channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913914,5640TC83.3% 593.8 / 118.7 270intergenic (‑436/+11)ybjE/aqpZputative surface protein/transmembrane water channel
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base C (131/94);  total (154/116)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 275 total aligned reads displayed.

TATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCAC  >  NC_000913/914474‑914656
                                                                                          |                                                                                            
tatGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc                                                                                            >  1:638293/1‑93 (MQ=255)
tatGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTc                                                                                            <  2:5089874/93‑1 (MQ=255)
tatGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTc                                                                                            >  2:2903291/1‑93 (MQ=255)
 atGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCctct                                                                                           <  1:4968973/93‑1 (MQ=255)
 atGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCctct                                                                                           >  2:1213161/1‑93 (MQ=255)
  tGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                                                                                          >  2:3337197/1‑93 (MQ=255)
  tGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTg                                                                                          >  2:2176838/1‑93 (MQ=255)
  tGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTg                                                                                          >  2:5450794/1‑93 (MQ=255)
   gTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGc                                                                                         <  2:7217106/93‑1 (MQ=255)
   gTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGc                                                                                         >  1:5446625/1‑93 (MQ=255)
     tGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                                                                                       <  1:3250139/93‑1 (MQ=255)
     tGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGt                                                                                       >  2:6886529/1‑93 (MQ=255)
        ttAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGttaa                                                                                    <  2:3537214/93‑1 (MQ=255)
        ttAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGttaa                                                                                    <  1:3669895/93‑1 (MQ=255)
          aGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGttaatt                                                                                  <  1:2504968/93‑1 (MQ=255)
            cTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGttaattaa                                                                                >  1:6646096/1‑93 (MQ=255)
             tCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTtaattaat                                                                               >  1:3145371/1‑93 (MQ=255)
              cATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATc                                                                              >  1:4864555/1‑93 (MQ=255)
              cATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATc                                                                              <  2:4566222/93‑1 (MQ=255)
               aTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCa                                                                             >  2:4020626/1‑93 (MQ=255)
               aTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCa                                                                             <  2:5355942/93‑1 (MQ=255)
               aTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCa                                                                             <  2:5752281/93‑1 (MQ=255)
               aTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCa                                                                             <  2:5891572/93‑1 (MQ=255)
                tGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc                                                                            >  2:4779540/1‑93 (MQ=255)
                tGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCAc                                                                            >  2:6960961/1‑93 (MQ=255)
                   aaGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCt                                                                         <  2:2948895/93‑1 (MQ=255)
                   aaGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCt                                                                         >  1:1597251/1‑93 (MQ=255)
                    aGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCtt                                                                        <  2:6599663/93‑1 (MQ=255)
                     ggAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt                                                                       <  1:269704/93‑1 (MQ=255)
                     ggAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCttt                                                                       <  1:7101514/93‑1 (MQ=255)
                     ggAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCttt                                                                       >  2:5840664/1‑93 (MQ=255)
                       aGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTc                                                                     <  1:218081/93‑1 (MQ=255)
                        gcgcAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTcc                                                                    >  2:3930152/1‑93 (MQ=255)
                        gcgcAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTcc                                                                    >  1:7393932/1‑93 (MQ=255)
                          gcAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCcag                                                                  >  1:3192645/1‑93 (MQ=255)
                          gcAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCcag                                                                  >  1:3220457/1‑93 (MQ=255)
                           caacaaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCcagc                                                                 >  2:2148495/1‑93 (MQ=255)
                             acaaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCcagcag                                                               <  2:3863984/93‑1 (MQ=255)
                             acaaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCcagcag                                                               >  2:5270425/1‑93 (MQ=255)
                              caaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAgg                                                              >  2:5037373/1‑93 (MQ=255)
                               aaGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAggg                                                             >  2:4537634/1‑93 (MQ=255)
                                aGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt                                                            <  2:375545/93‑1 (MQ=255)
                                aGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt                                                            >  2:7035577/1‑93 (MQ=255)
                                aGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt                                                            >  2:7504702/1‑93 (MQ=255)
                                   tgatgaATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCg                                                         <  1:1250848/93‑1 (MQ=255)
                                   tgatgaATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCg                                                         <  1:7501204/93‑1 (MQ=255)
                                       gaATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTaa                                                     <  1:4132492/93‑1 (MQ=255)
                                        aaTCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTaaa                                                    <  1:4410335/93‑1 (MQ=255)
                                        aaTCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTaaa                                                    >  1:6595893/1‑93 (MQ=255)
                                         aTCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAAt                                                   >  1:1103588/1‑93 (MQ=255)
                                              cATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGac                                              <  2:3578566/93‑1 (MQ=255)
                                              cATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGac                                              <  1:985972/93‑1 (MQ=255)
                                                tataACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGacca                                            >  2:1537414/1‑93 (MQ=255)
                                                    acaacaTCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGa                                        <  2:6386912/93‑1 (MQ=255)
                                                     caacaTCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGAt                                       <  1:899279/93‑1 (MQ=255)
                                                     caacaTCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGAt                                       <  1:3128149/93‑1 (MQ=255)
                                                      aacaTCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATa                                      >  1:4634809/1‑93 (MQ=255)
                                                      aacaTCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATa                                      >  2:5223749/1‑93 (MQ=255)
                                                         aTCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATg                                   >  2:3462180/1‑93 (MQ=255)
                                                         aTCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATg                                   >  1:7078074/1‑93 (MQ=255)
                                                         aTCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATg                                   >  2:3827540/1‑93 (MQ=255)
                                                          tCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgc                                  <  1:516320/93‑1 (MQ=255)
                                                          tCTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgc                                  >  2:3195628/1‑93 (MQ=255)
                                                           cTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgcc                                 <  1:7353568/93‑1 (MQ=255)
                                                           cTTAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgcc                                 <  1:2453026/93‑1 (MQ=255)
                                                            ttAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccg                                >  1:6224053/1‑93 (MQ=255)
                                                            ttAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccg                                >  1:1999694/1‑93 (MQ=255)
                                                            ttAAAAAAAGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATaccg                                >  2:4891920/1‑93 (MQ=255)
                                                             tAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccgc                               >  1:3712084/1‑93 (MQ=255)
                                                               aaaaaaGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGccgccg                             >  1:4036107/1‑93 (MQ=255)
                                                                aaaaaGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGa                            >  1:5612733/1‑93 (MQ=255)
                                                                 aaaaGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGAc                           >  1:5832687/1‑93 (MQ=255)
                                                                  aaaGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACa                          >  2:6215204/1‑93 (MQ=255)
                                                                  aaaGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACa                          >  2:863420/1‑93 (MQ=255)
                                                                   aaGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACaa                         >  1:6103141/1‑93 (MQ=255)
                                                                   aaGGACTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACaa                         >  1:151356/1‑93 (MQ=255)
                                                                    aGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAt                        <  1:7168121/93‑1 (MQ=255)
                                                                    aGGCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAt                        >  2:129037/1‑93 (MQ=255)
                                                                     ggCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAtt                       <  2:247755/93‑1 (MQ=255)
                                                                      gCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTg                      <  1:4064476/93‑1 (MQ=255)
                                                                      gCCTGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTg                      >  1:6935719/1‑93 (MQ=255)
                                                                       ccTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTgg                     >  1:2350650/1‑93 (MQ=255)
                                                                        cTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGc                    <  2:5382415/93‑1 (MQ=255)
                                                                         tGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGca                   >  1:3330716/1‑93 (MQ=255)
                                                                         tGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGca                   >  1:2911887/1‑93 (MQ=255)
                                                                         tGACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGca                   >  2:1531882/1‑93 (MQ=255)
                                                                          gACATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGcac                  >  1:5410077/1‑93 (MQ=255)
                                                                           aCATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGcacc                 >  2:188789/1‑93 (MQ=255)
                                                                           aCATTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGcacc                 <  1:5305551/93‑1 (MQ=255)
                                                                             aTTACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGcaccac               >  1:5565422/1‑93 (MQ=255)
                                                                               tACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCAccc             >  1:965267/1‑93 (MQ=255)
                                                                               tACGCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCAccc             >  1:383112/1‑93 (MQ=255)
                                                                                  gCCAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCaga          >  1:6631905/1‑93 (MQ=255)
                                                                                    cAGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCagaag        <  1:360356/93‑1 (MQ=255)
                                                                                     aGGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCagaaga       <  2:2035800/93‑1 (MQ=255)
                                                                                      ggCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAgaagaa      >  2:2298745/1‑93 (MQ=255)
                                                                                      ggCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAgaagaa      >  1:5019020/1‑93 (MQ=255)
                                                                                       gCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAAc     <  1:1832592/93‑1 (MQ=255)
                                                                                       gCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAAc     >  1:738503/1‑93 (MQ=255)
                                                                                          ttCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACcac  <  2:6320819/93‑1 (MQ=255)
                                                                                          |                                                                                            
TATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCAC  >  NC_000913/914474‑914656

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