Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,953,910 A→G N320S (AAC→AGC)  ilvA → threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,953,9100AG100.0% 822.5 / NA 257N320S (AAC→AGC) ilvAthreonine deaminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (134/123);  total (134/123)

Only 100 of 274 total aligned reads displayed.

CGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAG  >  NC_000913/3953820‑3953999
                                                                                          |                                                                                         
 gctggcgctggcGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCCTATTGTTTCCGGTGCCAACGGGAGCtt                                                                                        >  2:1957332/1‑93 (MQ=255)
  ctggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAAAGTGAGCTTc                                                                                       >  1:4847377/1‑93 (MQ=255)
   tggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTcc                                                                                      <  1:1429020/93‑1 (MQ=255)
    ggctctggcgGGAATGAAAAAATATATCGCCCTGCACAACAGTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCa                                                                                     <  1:3365434/93‑1 (MQ=255)
    ggcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCa                                                                                     <  1:2275772/93‑1 (MQ=255)
     gcgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCAc                                                                                    <  1:231564/93‑1 (MQ=255)
      cgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACg                                                                                   <  2:6535496/93‑1 (MQ=255)
      cgctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACg                                                                                   >  1:707769/1‑93 (MQ=255)
        ctggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCCACGTGAGCTTCCACGGc                                                                                 >  1:2677118/1‑93 (MQ=255)
         tggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGcc                                                                                >  2:3691164/1‑93 (MQ=255)
          ggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCt                                                                               >  1:4390620/1‑93 (MQ=255)
          ggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCt                                                                               >  1:1713566/1‑93 (MQ=255)
          ggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCt                                                                               >  1:4527203/1‑93 (MQ=255)
          ggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCt                                                                               <  1:6578079/93‑1 (MQ=255)
          ggcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCCTCCACGGCCt                                                                               >  1:430219/1‑93 (MQ=255)
           gcgGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTg                                                                              >  1:5862586/1‑93 (MQ=255)
             gggAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTgcg                                                                            <  1:3619536/93‑1 (MQ=255)
              ggAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTgcgc                                                                           <  1:29831/93‑1 (MQ=255)
                aaTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTa                                                                         >  1:5185201/1‑93 (MQ=255)
                aaTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTa                                                                         >  2:4787916/1‑93 (MQ=255)
                aaTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTa                                                                         >  2:5531078/1‑93 (MQ=255)
                 aTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCCACGTGAGCTTCCACGGCCTGCGCTAc                                                                        >  2:1479526/1‑93 (MQ=255)
                 aTGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTAc                                                                        >  1:2863257/1‑93 (MQ=255)
                    aaaaaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGtc                                                                     >  2:3595256/1‑93 (MQ=255)
                    aaaaaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGtc                                                                     <  2:6281192/93‑1 (MQ=255)
                     aaaaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGtct                                                                    >  2:3915088/1‑93 (MQ=255)
                      aaaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGtctc                                                                   <  1:1006519/93‑1 (MQ=255)
                       aaaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCa                                                                  <  1:5494299/93‑1 (MQ=255)
                        aaTATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAg                                                                 >  1:5291918/1‑93 (MQ=255)
                          tatatCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGaa                                                               <  1:2072374/93‑1 (MQ=255)
                             atCGCCCTGCACAACATTCGGGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGc                                                            <  2:2962915/93‑1 (MQ=255)
                               cGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGCACGCTg                                                          >  1:2921369/1‑93 (MQ=255)
                                gCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGc                                                         <  2:6717241/93‑1 (MQ=255)
                                  ccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGa                                                       >  2:4260944/1‑93 (MQ=255)
                                  ccTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGa                                                       >  1:1968888/1‑93 (MQ=255)
                                    tGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAAc                                                     >  2:5812848/1‑93 (MQ=255)
                                      cacaACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTg                                                   <  2:5850194/93‑1 (MQ=255)
                                      cacaACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTg                                                   <  1:1613218/93‑1 (MQ=255)
                                      cacaACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTg                                                   <  1:4741380/93‑1 (MQ=255)
                                        caacaTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCCACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTggg                                                 >  2:1901202/1‑93 (MQ=255)
                                          acaTTCGCGGCGAACGGCTGGCGCATATTCTTTCAGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCg                                               >  1:1093823/1‑93 (MQ=255)
                                            aTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGTGCTACGTCTCAGAACGCTGCGAACTGGGCGaa                                             <  1:3738629/93‑1 (MQ=255)
                                            aTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGaa                                             <  2:4716748/93‑1 (MQ=255)
                                            aTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGaa                                             <  1:2806190/93‑1 (MQ=255)
                                            aTTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGaa                                             <  1:2709645/93‑1 (MQ=255)
                                             ttCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAAc                                            >  2:5734929/1‑93 (MQ=255)
                                               cgcgGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAg                                          <  1:4427149/93‑1 (MQ=255)
                                               cgcgGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAg                                          >  2:3345737/1‑93 (MQ=255)
                                               cgcgGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAg                                          <  2:5247393/93‑1 (MQ=255)
                                                gcggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGc                                         <  1:3860722/93‑1 (MQ=255)
                                                gcggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTCGGAGAACAGc                                         >  1:1831139/1‑93 (MQ=255)
                                                 cggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCg                                        >  2:770873/1‑93 (MQ=255)
                                                  ggcgAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGGGAACAGCgg                                       >  1:5649300/1‑92 (MQ=255)
                                                     gAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGaa                                    <  2:257399/93‑1 (MQ=255)
                                                      aaCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAg                                   >  1:3922308/1‑93 (MQ=255)
                                                       aCGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGc                                  >  2:4031833/1‑93 (MQ=255)
                                                         ggCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgt                                <  1:603108/93‑1 (MQ=255)
                                                         ggCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgt                                <  2:6101728/93‑1 (MQ=255)
                                                         ggCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgt                                >  1:3982715/1‑93 (MQ=255)
                                                         ggCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgt                                <  1:5382251/93‑1 (MQ=255)
                                                          gCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgtt                               <  1:1299323/93‑1 (MQ=255)
                                                          gCTGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgtt                               <  2:5325620/93‑1 (MQ=255)
                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             >  1:2091563/1‑93 (MQ=255)
                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             <  1:1323002/93‑1 (MQ=255)
                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             <  1:2583459/93‑1 (MQ=255)
                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             <  2:3042453/93‑1 (MQ=255)
                                                            tGGCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGATACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgt                             <  1:5369148/93‑1 (MQ=255)
                                                             ggCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgtt                            <  1:6649234/93‑1 (MQ=255)
                                                             ggCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgtt                            >  1:3286939/1‑93 (MQ=255)
                                                             ggCGCATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCgttgtt                            <  1:5962671/93‑1 (MQ=255)
                                                              gcgcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGttgttg                           >  2:1715461/1‑93 (MQ=255)
                                                              gcgcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGttgttg                           <  1:856825/93‑1 (MQ=255)
                                                              gcgcATATTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGttgttg                           <  1:4040277/93‑1 (MQ=255)
                                                                  atatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCgg                       >  1:556820/1‑93 (MQ=255)
                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCgga                      >  1:4050099/1‑92 (MQ=255)
                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGt                      <  1:4964330/93‑1 (MQ=255)
                                                                   tatTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGt                      <  2:5005846/93‑1 (MQ=255)
                                                                    atTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTg                     >  2:4873733/1‑93 (MQ=255)
                                                                    atTCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTg                     <  1:1365195/93‑1 (MQ=255)
                                                                      tCTTTCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGAc                   >  2:5376801/1‑93 (MQ=255)
                                                                        tttCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCa                 >  1:3185711/1‑93 (MQ=255)
                                                                         ttCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCAt                <  2:3460143/93‑1 (MQ=255)
                                                                          tCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCAtt               >  2:6407652/1‑93 (MQ=255)
                                                                          tCCGGTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCAtt               <  2:1023626/93‑1 (MQ=255)
                                                                             ggTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCg            <  2:6646851/93‑1 (MQ=255)
                                                                             ggTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCg            <  2:2700508/93‑1 (MQ=255)
                                                                             ggTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCg            >  2:5560961/1‑93 (MQ=255)
                                                                             ggTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCg            >  1:341375/1‑93 (MQ=255)
                                                                             ggTGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCg            <  2:4512373/93‑1 (MQ=255)
                                                                               tGCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGga          >  2:6179239/1‑93 (MQ=255)
                                                                                gCCAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaa         >  1:4601257/1‑93 (MQ=255)
                                                                                 ccAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaag        <  1:5694836/93‑1 (MQ=255)
                                                                                 ccAACGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGgaag        >  2:323040/1‑93 (MQ=255)
                                                                                    aCGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaa     >  2:2412930/1‑93 (MQ=255)
                                                                                    aCGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaa     >  2:4201607/1‑93 (MQ=255)
                                                                                    aCGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaa     >  2:18398/1‑93 (MQ=255)
                                                                                     cGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaa    <  1:5130481/93‑1 (MQ=255)
                                                                                     cGTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaa    <  2:2530820/93‑1 (MQ=255)
                                                                                      gTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaaa   <  2:6350288/93‑1 (MQ=255)
                                                                                      gTGAGCTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGaaaaa   <  1:6242911/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCGGCGAACGGCTGGCGCATATTCTTTCCGGTGCCAACGTGAACTTCCACGGCCTGCGCTACGTCTCAGAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAG  >  NC_000913/3953820‑3953999

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