Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,402,284 T→C V7A (GTA→GCA)  purA → adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,402,2840TC80.5% 578.3 / 152.0 292V7A (GTA→GCA) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base C (120/115);  total (150/142)
Fisher's exact test for biased strand distribution p-value = 8.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 318 total aligned reads displayed.

GTGCAAAAAGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAGGGCGGTC  >  NC_000913/4402194‑4402376
                                                                                          |                                                                                            
 tGCAAAAAGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACt                                                                                           >  1:473053/1‑93 (MQ=255)
   cAAAAAGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTgg                                                                                         <  1:4487384/93‑1 (MQ=255)
   cAAAAAGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTgg                                                                                         <  2:2714103/93‑1 (MQ=255)
   cAAAAAGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTgg                                                                                         <  2:1252578/93‑1 (MQ=255)
      aaaGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCa                                                                                      <  1:7383905/93‑1 (MQ=255)
      aaaGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCa                                                                                      >  2:7048614/1‑93 (MQ=255)
      aaaGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCa                                                                                      <  1:1881217/93‑1 (MQ=255)
       aaGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCAc                                                                                     <  1:6385852/93‑1 (MQ=255)
         gTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCAccc                                                                                   >  1:4986191/1‑93 (MQ=255)
          tGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCa                                                                                  >  1:5698049/1‑93 (MQ=255)
            cTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAAt                                                                                <  2:899436/93‑1 (MQ=255)
             tGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATg                                                                               <  1:2023069/93‑1 (MQ=255)
             tGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATg                                                                               <  2:1789998/93‑1 (MQ=255)
              gTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATgg                                                                              >  1:6658607/1‑93 (MQ=255)
               tAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATggg                                                                             <  2:3134833/93‑1 (MQ=255)
               tAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATggg                                                                             >  1:3854746/1‑93 (MQ=255)
                 aCTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGt                                                                           <  2:1694549/93‑1 (MQ=255)
                 aCTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGt                                                                           <  2:7233406/93‑1 (MQ=255)
                 aCTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGt                                                                           <  2:5550914/93‑1 (MQ=255)
                  ctctGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTg                                                                          >  2:478554/1‑93 (MQ=255)
                   tctGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGa                                                                         >  2:5710729/1‑93 (MQ=255)
                   tctGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGa                                                                         >  1:2988900/1‑93 (MQ=255)
                     tGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACg                                                                       >  2:540597/1‑93 (MQ=255)
                     tGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACg                                                                       >  2:2441532/1‑93 (MQ=255)
                     tGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACg                                                                       >  1:160703/1‑93 (MQ=255)
                      gAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGa                                                                      >  1:662461/1‑93 (MQ=255)
                      gAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGa                                                                      <  2:5894184/93‑1 (MQ=255)
                       aaaaaGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGaa                                                                     <  1:873901/93‑1 (MQ=255)
                        aaaaGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGaag                                                                    >  1:6260690/1‑93 (MQ=255)
                           aGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGaaggta                                                                 >  2:3509897/1‑93 (MQ=255)
                            gCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGaaggtaa                                                                <  2:7529385/93‑1 (MQ=255)
                            gCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGaaggtaa                                                                >  2:3187272/1‑93 (MQ=255)
                                 ggTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGaaggtaaaggta                                                           <  2:1890696/93‑1 (MQ=255)
                                 ggTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGaaggtaaaggta                                                           <  1:5320830/93‑1 (MQ=255)
                                  gTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAaggtaaaggtaa                                                          >  2:1129110/1‑93 (MQ=255)
                                  gTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAaggtaaaggtaa                                                          <  2:4352283/93‑1 (MQ=255)
                                  gTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAaggtaaaggtaa                                                          <  2:1357304/93‑1 (MQ=255)
                                    aGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGa                                                        <  1:1788022/93‑1 (MQ=255)
                                    aGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGAAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGa                                                        >  1:3482806/1‑93 (MQ=255)
                                      aaTCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGAtc                                                      >  1:6678643/1‑93 (MQ=255)
                                       aTCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGAtcg                                                     >  2:5623163/1‑93 (MQ=255)
                                        tCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGAtcgt                                                    >  2:426229/1‑93 (MQ=255)
                                        tCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGAtcgt                                                    >  1:5615542/1‑93 (MQ=255)
                                        tCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGAtcgt                                                    >  1:2216815/1‑93 (MQ=255)
                                        tCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGAtcgt                                                    >  1:2898234/1‑93 (MQ=255)
                                        tCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGAtcgt                                                    <  1:3207836/93‑1 (MQ=255)
                                         ccATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAGAGGTAAGAtcgtc                                                   >  2:474496/1‑93 (MQ=255)
                                           aTTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGa                                                 >  2:4122092/1‑93 (MQ=255)
                                           aTTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGa                                                 <  2:3631112/93‑1 (MQ=255)
                                           aTTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGa                                                 >  2:2365134/1‑93 (MQ=255)
                                           aTTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGa                                                 >  1:3584517/1‑93 (MQ=255)
                                            tttttAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAt                                                >  1:4971149/1‑93 (MQ=255)
                                             ttttAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtc                                               <  2:6734434/93‑1 (MQ=255)
                                              tttAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtct                                              >  1:3122626/1‑93 (MQ=255)
                                              tttAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtct                                              <  2:5127639/93‑1 (MQ=255)
                                                tAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtcttc                                            <  2:4465002/93‑1 (MQ=255)
                                                 aaGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtcttct                                           <  1:1409815/93‑1 (MQ=255)
                                                 aaGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtcttct                                           <  1:4785862/93‑1 (MQ=255)
                                                 aaGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtcttct                                           <  1:946775/93‑1 (MQ=255)
                                                 aaGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGAtcttct                                           <  1:2958132/93‑1 (MQ=255)
                                                     aaaCGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTctgact                                       >  2:293743/1‑93 (MQ=255)
                                                       aCGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTctgactga                                     >  1:3663441/1‑93 (MQ=255)
                                                       aCGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTctgactga                                     >  2:7444351/1‑93 (MQ=255)
                                                       aCGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTctgactga                                     >  1:6958840/1‑93 (MQ=255)
                                                        cGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGaa                                    <  1:1028372/93‑1 (MQ=255)
                                                        cGGTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGaa                                    <  1:1319436/93‑1 (MQ=255)
                                                         ggTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAAc                                   >  1:4262077/1‑93 (MQ=255)
                                                         ggTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAAc                                   <  2:6242546/93‑1 (MQ=255)
                                                         ggTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAAc                                   >  2:2620522/1‑93 (MQ=255)
                                                          gTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACg                                  >  1:3943313/1‑93 (MQ=255)
                                                          gTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACg                                  <  2:2199959/93‑1 (MQ=255)
                                                          gTGATTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACg                                  >  2:3841045/1‑93 (MQ=255)
                                                            gATTTTGAAAAATGGGTAACAACCTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACggg                                >  2:6368659/1‑93 (MQ=255)
                                                             aTTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGc                               <  1:2127607/93‑1 (MQ=255)
                                                             aTTTTGAAAAATGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGc                               >  1:1097640/1‑93 (MQ=255)
                                                                   aaaaaTGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAata                         <  1:2614193/93‑1 (MQ=255)
                                                                   aaaaaTGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAata                         <  1:243509/93‑1 (MQ=255)
                                                                    aaaaTGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAatat                        >  1:6858655/1‑93 (MQ=255)
                                                                      aaTGGGTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATAtgt                      >  1:7052793/1‑93 (MQ=255)
                                                                          ggTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTa                  >  1:762941/1‑93 (MQ=255)
                                                                          ggTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTa                  <  2:5183234/93‑1 (MQ=255)
                                                                          ggTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTa                  <  2:6074933/93‑1 (MQ=255)
                                                                           gTAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTAc                 <  1:7192790/93‑1 (MQ=255)
                                                                            tAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACg                >  2:7151529/1‑93 (MQ=255)
                                                                            tAACAACGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACg                <  1:5560403/93‑1 (MQ=255)
                                                                             aacaacGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGc               <  1:148254/93‑1 (MQ=255)
                                                                             aacaacGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGc               >  2:1938974/1‑93 (MQ=255)
                                                                              acaacGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCt              >  2:400398/1‑93 (MQ=255)
                                                                              acaacGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCt              >  1:5425743/1‑93 (MQ=255)
                                                                              acaacGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCt              <  2:4925988/93‑1 (MQ=255)
                                                                               caacGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTa             <  1:4590679/93‑1 (MQ=255)
                                                                                aacGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTAc            <  1:593317/93‑1 (MQ=255)
                                                                                 acGTCGTCGCACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTAcc           >  2:7292178/1‑93 (MQ=255)
                                                                                   gtcgtcgcACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAg         <  1:1162599/93‑1 (MQ=255)
                                                                                   gtcgtcgcACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAg         <  1:501995/93‑1 (MQ=255)
                                                                                   gtcgtcgcACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAg         <  1:4030515/93‑1 (MQ=255)
                                                                                   gtcgtcgcACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAg         <  2:1520090/93‑1 (MQ=255)
                                                                                    tcgtcgcACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAgg        >  2:3229272/1‑93 (MQ=255)
                                                                                      gtcgtACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAGGGc      <  2:5791924/93‑1 (MQ=255)
                                                                                       tcgcACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAGGGCg     >  1:6560442/1‑93 (MQ=255)
                                                                                          |                                                                                            
GTGCAAAAAGTGCTGTAACTCTGAAAAAGCGATGGTAGAATCCATTTTTAAGCAAACGGTGATTTTGAAAAATGGGTAACAACGTCGTCGTACTGGGCACCCAATGGGGTGACGAAGGTAAAGGTAAGATCGTCGATCTTCTGACTGAACGGGCTAAATATGTTGTACGCTACCAGGGCGGTC  >  NC_000913/4402194‑4402376

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