Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,916,382 A→G V368V (GTT→GTC atpA ← membrane‑bound ATP synthase, F1 sector, alpha‑subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,916,3820AG99.6% 1550.5 / ‑3.7 486V368V (GTT→GTCatpAmembrane‑bound ATP synthase, F1 sector, alpha‑subunit
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (228/256);  total (228/258)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 504 total aligned reads displayed.

TGCCAGAGCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTAACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGC  >  NC_000913/3916292‑3916471
                                                                                          |                                                                                         
tgCCAGAGCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGAc                                                                                         <  1:3895852/93‑1 (MQ=255)
   cAGAGCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACTCGGGTTGACCGc                                                                                      >  2:1666949/1‑93 (MQ=255)
   cAGAGCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGAccac                                                                                      >  1:4284457/1‑91 (MQ=255)
    agagCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCa                                                                                     >  1:934258/1‑93 (MQ=255)
    agagCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCa                                                                                     <  1:732462/93‑1 (MQ=255)
       gCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGa                                                                                  >  1:6808955/1‑93 (MQ=255)
         ggTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACg                                                                                <  2:2846569/93‑1 (MQ=255)
          gTACGGATACCACCTGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGa                                                                               <  2:7713663/93‑1 (MQ=255)
          gTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGa                                                                               >  2:7794887/1‑93 (MQ=255)
           tACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGaa                                                                              <  2:9636440/93‑1 (MQ=255)
           tACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGaa                                                                              >  1:6688454/1‑93 (MQ=255)
              ggATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATAc                                                                           >  1:8342068/1‑93 (MQ=255)
               gATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATAcc                                                                          >  1:1650959/1‑93 (MQ=255)
               gATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATAcc                                                                          <  1:700039/93‑1 (MQ=255)
               gATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATAcc                                                                          <  2:5885971/93‑1 (MQ=255)
                aTACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCg                                                                         >  1:8566658/1‑93 (MQ=255)
                aTACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCg                                                                         <  1:2636187/93‑1 (MQ=255)
                aTACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCg                                                                         >  2:8583660/1‑93 (MQ=255)
                aTACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCg                                                                         <  2:729706/93‑1 (MQ=255)
                 tACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCgg                                                                        <  1:9659033/93‑1 (MQ=255)
                 tACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCgg                                                                        <  2:5665351/93‑1 (MQ=255)
                   ccaccGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCg                                                                      <  1:3415890/93‑1 (MQ=255)
                    caccGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGt                                                                     <  2:8160326/93‑1 (MQ=255)
                    caccGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGt                                                                     >  1:4645384/1‑93 (MQ=255)
                     accGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGtt                                                                    >  1:8821911/1‑93 (MQ=255)
                      ccGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGATATACTCGGGTTTACCGCAGGACGAATACCGGCGTTg                                                                   >  2:9024957/1‑93 (MQ=255)
                      ccGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTg                                                                   >  1:452622/1‑93 (MQ=255)
                         gACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAAc                                                                <  2:3858625/93‑1 (MQ=255)
                         gACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAAc                                                                >  2:2007624/1‑93 (MQ=255)
                         gACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAAc                                                                >  1:5645794/1‑93 (MQ=255)
                         gACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAAc                                                                >  1:6624271/1‑93 (MQ=255)
                          aCAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACa                                                               <  1:6323631/93‑1 (MQ=255)
                           cAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAg                                                              >  2:4475109/1‑93 (MQ=255)
                           cAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAg                                                              <  2:942774/93‑1 (MQ=255)
                           cAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAg                                                              >  1:3514957/1‑93 (MQ=255)
                            aGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAgg                                                             <  1:2383710/93‑1 (MQ=255)
                             gTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggt                                                            <  2:5025172/93‑1 (MQ=255)
                             gTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggt                                                            <  2:3158487/93‑1 (MQ=255)
                             gTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGAAGGACGAATACCGGAGTTGACCAggt                                                            >  1:7496674/1‑93 (MQ=255)
                              ttttttCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggtt                                                           <  1:9865946/93‑1 (MQ=255)
                              ttttttCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggtt                                                           <  1:1312932/93‑1 (MQ=255)
                              ttttttCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggtt                                                           <  1:3757208/93‑1 (MQ=255)
                               tttttCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACCggttg                                                          >  1:2861827/1‑93 (MQ=255)
                                ttttCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggttgg                                                         >  2:5232104/1‑93 (MQ=255)
                                 tttCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggttggt                                                        <  1:5220968/93‑1 (MQ=255)
                                 tttCATGATATTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAggttggt                                                        <  2:1837242/93‑1 (MQ=255)
                                   tCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATGCCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGttt                                                      <  2:3707062/93‑1 (MQ=255)
                                   tCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGttt                                                      <  1:6165297/93‑1 (MQ=255)
                                    cATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTc                                                     <  2:8020010/93‑1 (MQ=255)
                                     aTGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTcc                                                    <  2:1485342/93‑1 (MQ=255)
                                       gATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAg                                                  >  2:2346662/1‑93 (MQ=255)
                                       gATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAg                                                  <  2:3980238/93‑1 (MQ=255)
                                        aTCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAgg                                                 >  1:1160160/1‑93 (MQ=255)
                                          cTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGaa                                               >  1:5480766/1‑93 (MQ=255)
                                          cTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGaa                                               <  2:2939396/93‑1 (MQ=255)
                                           ttGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAg                                              >  1:1976847/1‑93 (MQ=255)
                                            tGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGa                                             <  2:10002655/93‑1 (MQ=255)
                                             ggTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCTGGTTGACCGCAGGACGAATACCGTCGTTGAACAGGTTGGTTTCCAGGAAGAt                                            <  2:9785587/93‑1 (MQ=255)
                                             ggTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGAt                                            <  1:756450/93‑1 (MQ=255)
                                              gTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATc                                           <  1:9764484/93‑1 (MQ=255)
                                              gTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATc                                           <  1:9457191/93‑1 (MQ=255)
                                                cTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTg                                         <  2:7542843/93‑1 (MQ=255)
                                                 tgtgCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTgc                                        >  1:7847551/1‑92 (MQ=255)
                                                 tgtgCTGCACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGa                                        >  2:6231631/1‑93 (MQ=255)
                                                    gctgcACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCa                                     <  1:2330198/93‑1 (MQ=255)
                                                      tgcACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATc                                   >  2:6172349/1‑93 (MQ=255)
                                                      tgcACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACAGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATc                                   >  1:8061570/1‑93 (MQ=255)
                                                       gcACCACCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCg                                  <  1:8387647/93‑1 (MQ=255)
                                                        caccacCAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCgg                                 >  2:4629198/1‑93 (MQ=255)
                                                         accaccAACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCTGt                                >  2:5446108/1‑93 (MQ=255)
                                                          ccaccaACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTa                               >  1:8943475/1‑93 (MQ=255)
                                                          ccaccaACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTa                               <  1:7122661/93‑1 (MQ=255)
                                                           caccaACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTaa                              <  2:8870339/93‑1 (MQ=255)
                                                           caccaACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTaa                              <  2:4812409/93‑1 (MQ=255)
                                                            accaACACGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAAt                             <  1:6448547/93‑1 (MQ=255)
                                                                 cacGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGgaga                        >  2:4318482/1‑93 (MQ=255)
                                                                  acGGGATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGAt                       <  2:5978414/93‑1 (MQ=255)
                                                                    gggATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTAGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTa                     >  1:601178/1‑93 (MQ=255)
                                                                     ggATACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTAc                    >  2:6744195/1‑93 (MQ=255)
                                                                       aTACGGAAATCCCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGt                  <  1:5277576/93‑1 (MQ=255)
                                                                       aTACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGt                  <  2:3073692/93‑1 (MQ=255)
                                                                        tACGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGtt                 >  1:3140558/1‑93 (MQ=255)
                                                                         aCGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTg                <  2:2554370/93‑1 (MQ=255)
                                                                          cGGAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTgg               >  2:7179746/1‑93 (MQ=255)
                                                                           ggAAATACCCGGGTTGACCGCAGGACGAATACCGGGGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGt              >  1:1798654/1‑93 (MQ=255)
                                                                           ggAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGt              <  2:4110076/93‑1 (MQ=255)
                                                                            gAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTc             <  2:6547368/93‑1 (MQ=255)
                                                                            gAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTc             >  1:793429/1‑93 (MQ=255)
                                                                            gAAATACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTc             <  2:3574278/93‑1 (MQ=255)
                                                                              aaTACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCgg           <  2:2179133/93‑1 (MQ=255)
                                                                               aTACCCGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGga          >  1:1003940/1‑93 (MQ=255)
                                                                                  cccGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGgaacg       <  1:4864059/93‑1 (MQ=255)
                                                                                  cccGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGgaacg       >  1:9051024/1‑93 (MQ=255)
                                                                                   ccGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGgaacga      >  2:7009023/1‑93 (MQ=255)
                                                                                    cGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGgaacgaa     >  2:83669/1‑93 (MQ=255)
                                                                                    cGGGTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGgaacgaa     >  2:3292495/1‑93 (MQ=255)
                                                                                     gggTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGGAATGGAGATTACGTTGGTCGgaacgaac    >  1:2161705/1‑93 (MQ=255)
                                                                                      ggTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGaacgaacg   >  2:7417402/1‑93 (MQ=255)
                                                                                       gTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGc  >  1:5008229/1‑93 (MQ=255)
                                                                                       gTTGACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGc  <  2:1310707/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGCCAGAGCGGTACGGATACCACCGGACAGTTTTTTCATGATCTTGGTCTGTGCTGCACCACCAACACGGGATACGGAAATACCCGGGTTAACCGCAGGACGAATACCGGCGTTGAACAGGTTGGTTTCCAGGAAGATCTGACCATCGGTAATGGAGATTACGTTGGTCGGAACGAACGC  >  NC_000913/3916292‑3916471

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