Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,004,502 T→C R49R (CGT→CGC ygeW → putative carbamoyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,004,5020TC99.5% 573.0 / NA 186R49R (CGT→CGCygeWputative carbamoyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (107/78);  minor base G (1/0);  total (108/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

AGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGC  >  NC_000913/3004412‑3004589
                                                                                          |                                                                                       
agcaTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCAc                                                                                       >  1:1205067/1‑93 (MQ=255)
   aTTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCcacccgcacccg                                                                                    >  2:975745/1‑93 (MQ=255)
    ttAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTAATTCCGTATTCCGCGACAACTCCacccgcacccgc                                                                                   >  2:1446524/1‑93 (MQ=255)
    ttAAAAGCACTGCGTGCTAAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCacccgcacccgc                                                                                   >  1:373572/1‑93 (MQ=255)
     tAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCt                                                                                  <  2:2654910/93‑1 (MQ=255)
       aaaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTtc                                                                                >  2:804237/1‑93 (MQ=255)
        aaGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTtct                                                                               >  2:4207109/1‑93 (MQ=255)
             cTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTAt                                                                          <  1:1702666/93‑1 (MQ=255)
             cTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTAt                                                                          <  2:1915475/93‑1 (MQ=255)
              tGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATg                                                                         >  1:3500041/1‑93 (MQ=255)
                cGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCt                                                                       <  2:1207715/93‑1 (MQ=255)
                cGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCt                                                                       <  1:363265/93‑1 (MQ=255)
                  tGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTc                                                                     <  2:577270/93‑1 (MQ=255)
                  tGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTc                                                                     >  1:4049575/1‑93 (MQ=255)
                      gAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCcgcg                                                                 >  1:155615/1‑93 (MQ=255)
                       aaaaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCgcgc                                                                >  2:1294683/1‑93 (MQ=255)
                         aaCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCtt                                                              >  2:1269308/1‑93 (MQ=255)
                          aCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTa                                                             <  1:3940910/93‑1 (MQ=255)
                          aCATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTa                                                             >  1:3863610/1‑93 (MQ=255)
                           cATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTaa                                                            >  2:900051/1‑93 (MQ=255)
                            aTCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAAc                                                           <  2:2198077/93‑1 (MQ=255)
                               tcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                        >  2:2979387/1‑93 (MQ=255)
                               tcAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTg                                                        >  1:2703085/1‑93 (MQ=255)
                                 aaCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCt                                                      >  2:3590818/1‑93 (MQ=255)
                                  aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                     <  1:1047788/93‑1 (MQ=255)
                                  aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                     <  1:2149070/93‑1 (MQ=255)
                                  aCCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCAACGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTc                                                     >  2:2624325/1‑93 (MQ=255)
                                   ccAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCg                                                    <  2:405948/93‑1 (MQ=255)
                                    cAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCgg                                                   >  1:3287191/1‑93 (MQ=255)
                                    cAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCgg                                                   >  1:1874051/1‑93 (MQ=255)
                                    cAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCgg                                                   >  2:412048/1‑93 (MQ=255)
                                     aaaGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGc                                                  >  1:3039478/1‑93 (MQ=255)
                                       aGTCTTTAATAGTGGATTAGGTATTTCCGCATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCt                                                >  2:2032892/1‑93 (MQ=255)
                                        gTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCCCCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCtt                                               <  2:4283893/93‑1 (MQ=255)
                                        gTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCtt                                               <  1:3574400/93‑1 (MQ=255)
                                        gTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCtt                                               >  1:2266062/1‑93 (MQ=255)
                                        gTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCtt                                               >  2:2764260/1‑93 (MQ=255)
                                        gTCTTTAATAGGGGATTAGGTATTTCCGTATTCCGAGACAACTCCACTCGGACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCtt                                               >  2:635764/1‑93 (MQ=255)
                                         tCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTg                                              >  2:223506/1‑93 (MQ=255)
                                         tCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTg                                              >  2:2295941/1‑93 (MQ=255)
                                          cTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGTTTAAACTGCTCGGCCTTGc                                             >  2:3317540/1‑93 (MQ=255)
                                          cTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGc                                             <  2:3236507/93‑1 (MQ=255)
                                          cTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGc                                             >  1:54234/1‑93 (MQ=255)
                                           tttAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGca                                            <  1:259736/93‑1 (MQ=255)
                                           tttAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGca                                            >  2:3496720/1‑93 (MQ=255)
                                            ttAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcac                                           >  2:2577186/1‑93 (MQ=255)
                                             tAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGcaca                                          >  1:1577321/1‑93 (MQ=255)
                                               aTAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaac                                        <  1:2352594/93‑1 (MQ=255)
                                                tAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCacaaca                                       >  2:2059897/1‑93 (MQ=255)
                                                 aGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCAcaacaa                                      <  2:1105055/93‑1 (MQ=255)
                                                  gTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                     >  1:3234770/1‑93 (MQ=255)
                                                  gTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                     >  1:3477068/1‑93 (MQ=255)
                                                  gTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAg                                     <  2:1602982/93‑1 (MQ=255)
                                                    ggATTAGGTATTTCCGTATTCCGCGACAACTCCTCCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAt                                   <  1:4555262/93‑1 (MQ=255)
                                                     gATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtc                                  >  1:1016228/1‑93 (MQ=255)
                                                       ttAGGTATTTCCTTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTAGGCTTCCGCGCTTAACCTGCTCGGCCTTGCAAAACAAGAtctc                                <  1:3094687/93‑1 (MQ=255)
                                                       ttAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtctc                                <  1:1135861/93‑1 (MQ=255)
                                                       ttAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGAtctc                                >  1:705631/1‑93 (MQ=255)
                                                        tAGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcg                               >  2:470811/1‑93 (MQ=255)
                                                         aGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcga                              >  1:4408805/1‑93 (MQ=255)
                                                         aGGTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcga                              <  2:2692699/93‑1 (MQ=255)
                                                          ggTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgac                             >  1:3972732/1‑93 (MQ=255)
                                                          ggTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAAATCTcgac                             <  2:4576322/93‑1 (MQ=255)
                                                           gTATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacg                            >  2:1507089/1‑93 (MQ=255)
                                                            tATTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTcgacga                           <  2:1967335/93‑1 (MQ=255)
                                                             aTTTCCGTATTCCGCGACAACTCCGCCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                          <  1:2659704/93‑1 (MQ=255)
                                                             aTTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                          <  1:4487357/93‑1 (MQ=255)
                                                             aTTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                          >  2:1839657/1‑93 (MQ=255)
                                                             aTTTCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGaa                          >  1:1753638/1‑93 (MQ=255)
                                                              tttCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                         >  1:289394/1‑93 (MQ=255)
                                                              tttCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                         <  2:681430/93‑1 (MQ=255)
                                                              tttCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                         >  1:465775/1‑93 (MQ=255)
                                                              ttgccGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAg                         <  2:4490060/90‑1 (MQ=255)
                                                               ttCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                        >  1:4380844/1‑93 (MQ=255)
                                                               ttCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                        <  2:1681009/93‑1 (MQ=255)
                                                               ttCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAgg                        >  1:500549/1‑93 (MQ=255)
                                                                tCCGTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGc                       <  2:3758488/93‑1 (MQ=255)
                                                                   gTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaaa                    <  1:2121387/93‑1 (MQ=255)
                                                                   gTATTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCaaa                    <  1:3756462/93‑1 (MQ=255)
                                                                     aTTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAAta                  >  1:2104024/1‑92 (MQ=255)
                                                                     aTTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATc                  >  1:4190142/1‑93 (MQ=255)
                                                                     aTTCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATc                  >  2:4355042/1‑93 (MQ=255)
                                                                     aTTCCGCGACAACACCATCCGCACCTGCTTCTCTTATGCTTCCGACCTTACCCTGCTCGGCCTTGCACAAGCAGATCTCGACGAAGGCAACTc                  >  2:1205826/1‑93 (MQ=255)
                                                                      ttCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACACCAAGACCTCGACGAAGGCAAATca                 >  2:2842392/1‑93 (MQ=255)
                                                                       tCCGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATcac                <  2:3082368/93‑1 (MQ=255)
                                                                        ccGCGACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATcaca               <  2:2884124/93‑1 (MQ=255)
                                                                         cgcgACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaa              <  2:3922032/93‑1 (MQ=255)
                                                                         cgcgACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaa              <  1:2562583/93‑1 (MQ=255)
                                                                          gcgACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACaaa             >  2:597630/1‑93 (MQ=255)
                                                                           cgACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAAt            <  1:2540830/93‑1 (MQ=255)
                                                                           cgACAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAAt            <  1:979931/93‑1 (MQ=255)
                                                                              cAACTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGc         >  1:911128/1‑93 (MQ=255)
                                                                              cAACCCCCCCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGc         <  2:2967462/93‑1 (MQ=255)
                                                                               aaCTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCt        >  2:1370918/1‑93 (MQ=255)
                                                                                 cTCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCa      <  2:2487526/93‑1 (MQ=255)
                                                                                  tCCACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCCCTCAc     >  1:611671/1‑93 (MQ=255)
                                                                                   ccACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACg    <  1:3135346/93‑1 (MQ=255)
                                                                                   ccACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACg    <  1:4611390/93‑1 (MQ=255)
                                                                                   ccACCCGCACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACg    >  2:2972996/1‑93 (MQ=255)
                                                                                    cacccgcacccgCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACgg   <  2:1152870/93‑1 (MQ=255)
                                                                                          |                                                                                       
AGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCTCGGCCTTGCACAACAAGATCTCGACGAAGGCAAATCACAAATCGCTCACGGC  >  NC_000913/3004412‑3004589

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