Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 80,603 (T)5→6 coding (262/1401 nt) leuC ← 3‑isopropylmalate isomerase (dehydratase) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091380,5981.T99.2% 991.6 / ‑1.8 268coding (267/1401 nt)leuC3‑isopropylmalate isomerase (dehydratase) subunit
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (125/140);  total (126/142)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 308 total aligned reads displayed.

GGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATG  >  NC_000913/80506‑80685
                                                                                             |                                                                                       
ggCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAg                                                                                          >  1:6667781/1‑93 (MQ=255)
    aGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGtttt                                                                                     >  1:6720714/1‑93 (MQ=255)
    aGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGtttt                                                                                     <  2:5866402/93‑1 (MQ=255)
    aGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGtttt                                                                                     >  1:1904159/1‑93 (MQ=255)
     ggTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGttttt                                                                                    <  1:6041074/93‑1 (MQ=255)
     ggTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGttttt                                                                                    >  2:713652/1‑93 (MQ=255)
      gTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGtttttt                                                                                   >  1:5551133/1‑92 (MQ=255)
       tGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGttttttg                                                                                  <  1:2844040/93‑3 (MQ=255)
           gCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCa                                                                               >  1:3083806/1‑93 (MQ=255)
           gCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCa                                                                               <  1:717005/93‑1 (MQ=255)
            cccTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAg                                                                              >  1:2360662/1‑93 (MQ=255)
             ccTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGt                                                                             >  2:3410968/1‑93 (MQ=255)
              cTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGtt                                                                            <  1:5761678/93‑1 (MQ=255)
                gTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTcc                                                                          <  2:3497485/93‑1 (MQ=255)
                gTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTcc                                                                          <  2:3775577/93‑1 (MQ=255)
                 ttCCGGCCCCATTGCGGGGACGATCCGCTGATACGGGTGGTTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCt                                                                         <  1:1331616/93‑1 (MQ=255)
                 ttCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCt                                                                         <  2:1012951/93‑1 (MQ=255)
                  tCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTg                                                                        >  2:3833304/1‑93 (MQ=255)
                    cGGCCACATTACGTGCACGACCCCCTGATACGGGTGATTCAGGTCATACACTTCGACGCCAAATTCTTTGCCGTTTTTTGATCAGTTCCTGCa                                                                      <  1:5981830/93‑1 (MQ=255)
                      gCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATc                                                                    <  2:4350910/93‑1 (MQ=255)
                      gCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATc                                                                    >  2:5242656/1‑93 (MQ=255)
                        cccATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTg                                                                  >  1:6278336/1‑93 (MQ=255)
                        cccATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTg                                                                  <  2:2523555/93‑1 (MQ=255)
                         ccATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTgg                                                                 >  2:3929814/1‑93 (MQ=255)
                         ccATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTgg                                                                 >  2:5608300/1‑93 (MQ=255)
                         ccATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTgc                                                                 >  1:4418898/1‑92 (MQ=255)
                          cATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGa                                                                >  2:1278618/1‑93 (MQ=255)
                          cATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGa                                                                <  2:4805548/93‑1 (MQ=255)
                              aCGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATAcg                                                            <  2:5351141/93‑1 (MQ=255)
                                gTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATAcgcg                                                          >  2:5131278/1‑93 (MQ=255)
                                 tGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACgcgc                                                         >  1:6121339/1‑93 (MQ=255)
                                 tGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACgcgc                                                         >  1:1270230/1‑93 (MQ=255)
                                 tGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACgcgc                                                         <  2:6486541/93‑1 (MQ=255)
                                  ggACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGAAGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGcc                                                        >  1:2492382/1‑93 (MQ=255)
                                  agACGATCCCCGGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGcc                                                        <  1:5237074/92‑1 (MQ=255)
                                   gACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGAGCAGTTCCTGCATCTGGATACGCGCCa                                                       >  1:5392063/1‑93 (MQ=255)
                                    gcGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCAt                                                      <  1:4451596/92‑1 (MQ=255)
                                    aCGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCAt                                                      >  2:4185760/1‑93 (MQ=255)
                                    aCGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCAt                                                      <  2:4917603/93‑1 (MQ=255)
                                     cGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCAtt                                                     <  2:630740/93‑1 (MQ=255)
                                      gATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCAttt                                                    <  1:2894975/93‑1 (MQ=255)
                                         ccccTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCAc                                                 >  2:3586631/1‑93 (MQ=255)
                                          cccTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCAcc                                                >  2:4463397/1‑93 (MQ=255)
                                           ccTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCg                                               >  2:3192719/1‑93 (MQ=255)
                                            cTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGc                                              >  1:4550114/1‑93 (MQ=255)
                                             tGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCa                                             >  2:2392441/1‑93 (MQ=255)
                                                   gggTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCAtt                                       <  2:5500333/93‑1 (MQ=255)
                                                    ggTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTa                                      >  1:5591621/1‑93 (MQ=255)
                                                        aTTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGt                                  <  2:6560173/93‑1 (MQ=255)
                                                          tCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCt                                >  2:6195332/1‑93 (MQ=255)
                                                          tCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCt                                <  2:4296594/93‑1 (MQ=255)
                                                           cAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCtt                               >  1:4111056/1‑93 (MQ=255)
                                                           cAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCtt                               <  2:2652321/93‑1 (MQ=255)
                                                            aGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCttt                              >  1:1407287/1‑93 (MQ=255)
                                                            aGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCttt                              <  2:930092/93‑1 (MQ=255)
                                                             ggTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTg                             <  1:4279356/93‑1 (MQ=255)
                                                             ggTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTg                             <  1:5006369/93‑1 (MQ=255)
                                                             agTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTg                             <  2:3261245/92‑1 (MQ=255)
                                                               tCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGt                           >  1:4776442/1‑93 (MQ=255)
                                                                cATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTc                          <  1:2289200/93‑1 (MQ=255)
                                                                cATACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTc                          <  1:6473183/93‑1 (MQ=255)
                                                                 aTACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCt                         <  1:5859677/93‑1 (MQ=255)
                                                                 aTACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCt                         >  2:4614236/1‑93 (MQ=255)
                                                                 aTACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCt                         <  1:1024121/93‑1 (MQ=255)
                                                                  tACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTg                        <  2:2027544/93‑1 (MQ=255)
                                                                  tACAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTg                        >  1:5679098/1‑93 (MQ=255)
                                                                   tcAGTTCGACGCCCAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGCATACGCGCCATTTCGCCGCAGGCATTAATGTCTTTGGTCTgg                       <  1:4034339/92‑1 (MQ=255)
                                                                    cAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTggg                      >  2:1870984/1‑93 (MQ=255)
                                                                    cAGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTggg                      <  2:1656863/93‑1 (MQ=255)
                                                                     aGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGt                     <  1:1185722/93‑1 (MQ=255)
                                                                     aGTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGt                     <  1:334558/93‑1 (MQ=255)
                                                                      gTTCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTa                    >  1:1939229/1‑93 (MQ=255)
                                                                        tCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTaga                  >  2:6355523/1‑93 (MQ=255)
                                                                        tCGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTaga                  <  2:1422423/93‑1 (MQ=255)
                                                                         cGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTagag                 >  2:2160571/1‑93 (MQ=255)
                                                                         cGACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTagag                 >  2:5466240/1‑93 (MQ=255)
                                                                          gACTCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAgaga                >  2:6720411/1‑93 (MQ=255)
                                                                          gACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAgaga                >  2:3697922/1‑93 (MQ=255)
                                                                          gACGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAgaga                <  1:3838071/93‑1 (MQ=255)
                                                                           aCGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGAc               <  1:2617758/93‑1 (MQ=255)
                                                                           aCGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGAc               >  2:5751207/1‑93 (MQ=255)
                                                                            cGCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACg              <  1:6346891/93‑1 (MQ=255)
                                                                             gCCAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGt             <  2:1318706/93‑1 (MQ=255)
                                                                              ccAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGtt            <  1:3524479/93‑1 (MQ=255)
                                                                              ccAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGtt            >  1:3673929/1‑93 (MQ=255)
                                                                               cAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTtg           >  2:2339767/1‑93 (MQ=255)
                                                                               cAAATTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTtg           >  1:1513126/1‑93 (MQ=255)
                                                                                aaaTTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTtgt          >  2:337513/1‑93 (MQ=255)
                                                                                aaaTTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTtgt          >  1:3896364/1‑93 (MQ=255)
                                                                                aaaTTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTtgt          >  1:896356/1‑93 (MQ=255)
                                                                                  aTTCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGa        >  2:603070/1‑93 (MQ=255)
                                                                                   ttCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGAt       <  2:2250543/93‑1 (MQ=255)
                                                                                   ttCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGAt       <  1:1837644/93‑1 (MQ=255)
                                                                                    tCTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATc      <  1:1664965/93‑1 (MQ=255)
                                                                                    tCTTTGAAGTTTTTTGATCAGTTCCTGCATCTGAATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATc      <  1:2999727/93‑1 (MQ=255)
                                                                                     cTTTGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATcc     >  1:3117638/1‑93 (MQ=255)
                                                                                      tttGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCa    >  2:1169086/1‑93 (MQ=255)
                                                                                       ttGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCAt   <  1:1976779/93‑1 (MQ=255)
                                                                                       ttGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCAt   <  2:4044237/93‑1 (MQ=255)
                                                                                        tGCAGTTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATg  <  1:5181020/93‑1 (MQ=255)
                                                                                             |                                                                                       
GGCAAGGTGACGCCCTGTTCCGGCCCCATTACGTGGACGATCCCCTGATACGGGTGATTCAGGTCATACAGTTCGACGCCAAATTCTTTGCAGTTTTTGATCAGTTCCTGCATCTGGATACGCGCCATTTCACCGCAGGCATTAATGTCTTTGGTCTGGGTAGAGACGTTGTGATCCATG  >  NC_000913/80506‑80685

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