Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,689,482 A→G T303A (ACC→GCC)  manA → mannose‑6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,4820AG90.5% 1152.8 / 85.7 348T303A (ACC→GCC) manAmannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (14/19);  new base G (155/160);  total (169/179)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 349 total aligned reads displayed.

GATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGAAACCACCATTAG  >  NC_000913/1689383‑1689582
                                                                                                   |                                                                                                    
gATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGAc                                                                                                     >  2:3057311/1‑101 (MQ=255)
    aCGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCAGCTAACCAGTTGTTGGCCCAg                                                                                                 >  1:3057321/1‑101 (MQ=255)
         ctgcGTGCGGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGt                                                                                            >  2:3057324/1‑101 (MQ=255)
          tgcgtgcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTg                                                                                           <  1:3056693/101‑1 (MQ=255)
           gcgtgcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGa                                                                                          <  2:3056806/101‑1 (MQ=255)
           gcgtgcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGa                                                                                          <  1:3056836/101‑1 (MQ=255)
            cgtgcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGaa                                                                                         <  2:3056993/101‑1 (MQ=255)
             gtgcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGaaa                                                                                        >  2:3057329/1‑101 (MQ=255)
             gtgcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGaaa                                                                                        >  1:3057328/1‑101 (MQ=255)
             gtgcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGaa                                                                                         >  1:3057333/1‑100 (MQ=255)
               gcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACa                                                                                      <  2:3056934/101‑1 (MQ=255)
               gcgGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACa                                                                                      <  1:3057317/101‑1 (MQ=255)
                 gggTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAg                                                                                    <  2:3057032/101‑1 (MQ=255)
                   gTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGt                                                                                  <  2:3056912/101‑1 (MQ=255)
                    tCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGt                                                                                            <  2:3056867/90‑1 (MQ=255)
                    tCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTg                                                                                 >  1:3057339/1‑101 (MQ=255)
                    tCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTg                                                                                 <  1:3056731/101‑1 (MQ=255)
                      tGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCa                                                                               <  2:3056886/101‑1 (MQ=255)
                      tGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCa                                                                               >  2:3057343/1‑101 (MQ=255)
                       gACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAg                                                                              >  1:3057345/1‑101 (MQ=255)
                         cGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGaa                                                                            >  1:3057350/1‑101 (MQ=255)
                         cGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGaa                                                                            <  1:3056389/101‑1 (MQ=255)
                          gCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAAc                                                                           <  2:3056848/101‑1 (MQ=255)
                          gCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAAc                                                                           >  2:3057352/1‑101 (MQ=255)
                            cTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTg                                                                         <  2:3056959/101‑1 (MQ=255)
                             tAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGaa                                                                                         >  1:3057358/1‑84 (MQ=255)
                              aaaTACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGc                                                                                <  2:3056530/92‑1 (MQ=255)
                              aaaTACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAAc                                                                           >  2:3057359/1‑97 (MQ=255)
                               aaTACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGAc                                                                      >  2:3057364/1‑101 (MQ=255)
                                aTACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACt                                                                     >  1:3057365/1‑101 (MQ=255)
                                  aCATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTc                                                                   <  2:3056823/101‑1 (MQ=255)
                                  aCATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTc                                                                   <  1:3056961/101‑1 (MQ=255)
                                   cATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTcc                                                                  >  1:3057371/1‑101 (MQ=255)
                                    aTTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTccc                                                                 >  2:3057374/1‑101 (MQ=255)
                                    aTTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTccc                                                                 <  1:3056785/101‑1 (MQ=255)
                                    aTTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTccc                                                                 <  2:3056989/101‑1 (MQ=255)
                                     ttGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCg                                                                >  1:3057375/1‑101 (MQ=255)
                                      tGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGa                                                               <  2:3056768/101‑1 (MQ=255)
                                      tGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGa                                                               <  2:3057134/101‑1 (MQ=255)
                                       gATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGAt                                                              <  1:3056770/101‑1 (MQ=255)
                                        atatTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGAtt                                                             <  1:3056880/101‑1 (MQ=255)
                                        atatTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGAtt                                                             >  2:3057380/1‑101 (MQ=255)
                                          atTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTcc                                                           >  1:3057383/1‑101 (MQ=255)
                                          atTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTcc                                                           <  1:3056865/101‑1 (MQ=255)
                                           ttCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCa                                                          >  2:3057386/1‑101 (MQ=255)
                                              cGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCa                                                          >  1:3057392/1‑98 (MQ=255)
                                              cGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTg                                                       >  2:3057390/1‑101 (MQ=255)
                                               ggAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGAc                                                                      <  2:3056749/85‑1 (MQ=255)
                                                 aaCTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgat                                                    >  1:3057399/1‑101 (MQ=255)
                                                  aCTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgatg                                                   >  2:3057401/1‑101 (MQ=255)
                                                  aCTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgatg                                                   <  2:3056976/101‑1 (MQ=255)
                                                   cTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgatga                                                  <  1:3056742/101‑1 (MQ=255)
                                                    tGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGgatgat                                                 >  1:3057404/1‑101 (MQ=255)
                                                     ggTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAtt                                                >  2:3057406/1‑101 (MQ=255)
                                                      gTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAttt                                               <  1:3056640/101‑1 (MQ=255)
                                                      gTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAttt                                               <  1:3056802/101‑1 (MQ=255)
                                                       ttGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAtttt                                              <  1:3057024/101‑1 (MQ=255)
                                                          ccAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGcc                                           <  2:3056919/101‑1 (MQ=255)
                                                           cAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCt                                          >  2:3057413/1‑101 (MQ=255)
                                                             aTGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtc                                        >  2:3057416/1‑101 (MQ=255)
                                                              tgtgAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtct                                       >  1:3057418/1‑101 (MQ=255)
                                                              tgtgAAATTCGAAGCCAAACCGGCTAACCAGTTGTCGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTtct                                       <  1:3056889/101‑1 (MQ=255)
                                                                 gAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGc                                    >  1:3057422/1‑101 (MQ=255)
                                                                 gAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGc                                    <  1:3056810/101‑1 (MQ=255)
                                                                    aTTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGc                                 >  1:3057426/1‑101 (MQ=255)
                                                                    aTTCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGc                                 >  2:3057429/1‑101 (MQ=255)
                                                                     ttCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCa                                >  1:3057432/1‑101 (MQ=255)
                                                                     ttCGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCa                                <  1:3056927/101‑1 (MQ=255)
                                                                       cGAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATg                              <  2:3057076/101‑1 (MQ=255)
                                                                       cGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATg                              <  1:3056933/101‑1 (MQ=255)
                                                                       cGAAGCCAAACCGGATAACCAGTTGTTGGCCCAGCCg                                                                                              >  2:3057434/1‑37 (MQ=255)
                                                                        gAAGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGAtttt                                              <  2:3056903/84‑1 (MQ=255)
                                                                         aaGCCAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGAc                            >  2:3057439/1‑101 (MQ=255)
                                                                            ccAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCtt                         <  2:3056783/101‑1 (MQ=255)
                                                                            ccAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCtt                         <  1:3056896/101‑1 (MQ=255)
                                                                            ccAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCtt                         >  2:3057443/1‑101 (MQ=255)
                                                                             cAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTa                        >  1:3057445/1‑101 (MQ=255)
                                                                             cAAACCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTa                        <  1:3057055/101‑1 (MQ=255)
                                                                              aaaCCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAg                       <  2:3057247/101‑1 (MQ=255)
                                                                               aaCCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGt                      <  2:3056689/101‑1 (MQ=255)
                                                                                aCCGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTg                     >  1:3057449/1‑101 (MQ=255)
                                                                                 ccGGCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGa                    <  2:3056947/101‑1 (MQ=255)
                                                                                   ggCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTg                                                                                 >  1:3057457/1‑38 (MQ=255)
                                                                                   ggCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATa                  >  2:3057454/1‑101 (MQ=255)
                                                                                   ggCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATa                  <  2:3056603/101‑1 (MQ=255)
                                                                                    gCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATaa                 <  1:3057275/101‑1 (MQ=255)
                                                                                    gCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATaa                 <  1:3056977/101‑1 (MQ=255)
                                                                                    gCTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATaa                 <  2:3057169/101‑1 (MQ=255)
                                                                                     cTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATaaa                <  1:3056818/101‑1 (MQ=255)
                                                                                     cTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATaaa                <  1:3057005/101‑1 (MQ=255)
                                                                                     cTAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGACTTTGCCTTCTCGCTGCATGACCTTAGTGATaaa                >  2:3057460/1‑101 (MQ=255)
                                                                                      tAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAg               <  2:3057088/101‑1 (MQ=255)
                                                                                      tAACCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAg               >  1:3057464/1‑101 (MQ=255)
                                                                                       aaCCAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGa              >  1:3057468/1‑101 (MQ=255)
                                                                                         ccAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGaaa            <  1:3057133/101‑1 (MQ=255)
                                                                                         ccAGTTGTTGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGaaa            >  2:3057470/1‑101 (MQ=255)
                                                                                            gttgttGGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGAAacca         >  1:3057475/1‑101 (MQ=255)
                                                                                             ttgttgGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGAAaccac        >  2:3057478/1‑101 (MQ=255)
                                                                                              tgttgGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTg                                                       >  1:3057481/1‑53 (MQ=255)
                                                                                              tgttgGCCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGAAaccacc       <  2:3057046/101‑1 (MQ=255)
                                                                                                   |                                                                                                    
GATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGATATTCCGGAACTGGTTGCCAATGTGAAATTCGAAGCCAAACCGGCTAACCAGTTGTTGACCCAGCCGGTGAAACAAGGTGCAGAACTGGACTTCCCGATTCCAGTGGATGATTTTGCCTTCTCGCTGCATGACCTTAGTGATAAAGAAACCACCATTAG  >  NC_000913/1689383‑1689582

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