Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,871,984 A→G E129E (GAA→GAG yeaK → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,871,9840AG98.6% 427.6 / ‑3.3 139E129E (GAA→GAGyeaKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (73/64);  total (74/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 140 total aligned reads displayed.

GATGAACTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGA  >  NC_000913/1871895‑1872073
                                                                                         |                                                                                         
 aTGAACTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGAt                                                                                       >  2:3215458/1‑93 (MQ=255)
 aTGAACTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGAt                                                                                       <  2:2011445/93‑1 (MQ=255)
 aTGAACGCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGAt                                                                                       <  2:48459/93‑1 (MQ=255)
  tGAACTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGAtt                                                                                      >  2:2564527/1‑93 (MQ=255)
  tGAACTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGAtt                                                                                      >  1:3262092/1‑93 (MQ=255)
  tGAACTCACAGGCTGGTGCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGAtt                                                                                      <  2:45153/93‑1 (MQ=255)
   gAACTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGAttt                                                                                     <  1:3382811/93‑1 (MQ=255)
    aaCTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTtg                                                                                    <  1:4059298/93‑1 (MQ=255)
     aCTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTCCTCTTCGAGCGATTtga                                                                                   >  2:1298473/1‑93 (MQ=255)
       tCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTtgatg                                                                                 >  1:4539448/1‑93 (MQ=255)
        cacaGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTtgatga                                                                                <  2:4065530/93‑1 (MQ=255)
        cacaGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTtgatga                                                                                >  2:215516/1‑93 (MQ=255)
          caGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGaaa                                                                              <  1:3024761/93‑1 (MQ=255)
             gCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCg                                                                           >  1:3250520/1‑93 (MQ=255)
              cTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGc                                                                          >  1:2200551/1‑93 (MQ=255)
              cTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGc                                                                          >  2:1849000/1‑93 (MQ=255)
               tgtgTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGct                                                                         >  1:839741/1‑92 (MQ=255)
                gtgtCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCt                                                                        <  1:3944876/93‑1 (MQ=255)
                gtgtCTGCGGCGCAACCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCt                                                                        <  2:3699813/93‑1 (MQ=255)
                  gtCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTc                                                                      <  2:4138699/93‑1 (MQ=255)
                   tCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCa                                                                     <  1:200615/93‑1 (MQ=255)
                     ttCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAAt                                                                   >  2:4078841/1‑93 (MQ=255)
                     ttCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAAt                                                                   >  2:3483729/1‑93 (MQ=255)
                     ttCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAAt                                                                   <  2:2969445/93‑1 (MQ=255)
                     ttCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAAt                                                                   <  1:2238869/93‑1 (MQ=255)
                        ggCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCa                                                                <  1:1538989/93‑1 (MQ=255)
                        ggCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCa                                                                <  1:3713850/93‑1 (MQ=255)
                         gcgcGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAg                                                               <  1:3089804/93‑1 (MQ=255)
                            cgATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCa                                                            <  2:237767/93‑1 (MQ=255)
                            cgATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCa                                                            >  2:483310/1‑93 (MQ=255)
                             gATCCCCCCTTTCAGCTTCCATCCAACACTCAACCTGGTTGCCGTCCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCAt                                                           <  1:1792597/93‑1 (MQ=255)
                             gATCCCCCCTTTCAGCTTCAATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCAt                                                           >  2:772597/1‑93 (MQ=255)
                              aTCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATg                                                          >  2:449915/1‑93 (MQ=255)
                               tCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGc                                                         <  2:3107388/93‑1 (MQ=255)
                               tCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGc                                                         >  1:487719/1‑93 (MQ=255)
                               tCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGc                                                         >  1:2675871/1‑93 (MQ=255)
                                ccccccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCt                                                        <  2:2825192/93‑1 (MQ=255)
                                 cccccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTg                                                       >  2:1276720/1‑93 (MQ=255)
                                 cccccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTg                                                       >  1:2389235/1‑93 (MQ=255)
                                  ccccTTTCCGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTgg                                                      >  2:4584069/1‑93 (MQ=255)
                                  ccccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTgg                                                      >  2:2125048/1‑93 (MQ=255)
                                  ccccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTgg                                                      >  1:3478084/1‑93 (MQ=255)
                                   cccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGa                                                     <  2:1533533/93‑1 (MQ=255)
                                    ccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGAt                                                    <  1:4018599/93‑1 (MQ=255)
                                    ccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGAt                                                    <  2:2097919/93‑1 (MQ=255)
                                    ccTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGAt                                                    >  2:751942/1‑93 (MQ=255)
                                      tttCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATaa                                                  <  2:1925133/93‑1 (MQ=255)
                                       ttCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATaaa                                                 <  2:2171068/93‑1 (MQ=255)
                                       ttCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATaaa                                                 <  1:4482115/93‑1 (MQ=255)
                                         cAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATc                                               <  2:2786534/93‑1 (MQ=255)
                                          aGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATcc                                              <  1:3415622/93‑1 (MQ=255)
                                            cTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGt                                            >  1:1047547/1‑93 (MQ=255)
                                             ttCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGtt                                           >  1:1261887/1‑93 (MQ=255)
                                             ttCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGtt                                           <  1:1509124/93‑1 (MQ=255)
                                             ttCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGtt                                           >  2:43178/1‑93 (MQ=255)
                                             ttCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGtt                                           >  2:3187056/1‑93 (MQ=255)
                                                catccaAAACTCAAACTGGTTTCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTAtt                                        <  1:4190230/93‑1 (MQ=255)
                                                catccaAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTAtt                                        >  1:4134634/1‑93 (MQ=255)
                                                 atccaAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTAttt                                       <  1:979058/93‑1 (MQ=255)
                                                 atccaAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTAttt                                       >  2:4361613/1‑93 (MQ=255)
                                                  tccaAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGCTCAATCCGTTAtttt                                      >  1:4642079/1‑93 (MQ=255)
                                                  tccaAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTAtttt                                      <  2:1360578/93‑1 (MQ=255)
                                                  tccaAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTAtttt                                      >  2:3408896/1‑93 (MQ=255)
                                                  tccaAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTAtttt                                      <  2:2541728/93‑1 (MQ=255)
                                                     aaaaCTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGaa                                   >  2:3641027/1‑93 (MQ=255)
                                                       aaCTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGaaaa                                 <  1:3683585/93‑1 (MQ=255)
                                                        aCTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAAc                                <  1:986875/93‑1 (MQ=255)
                                                         cTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAAcc                               >  1:4378302/1‑93 (MQ=255)
                                                         cTCAAACTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAAcc                               <  2:1088270/93‑1 (MQ=255)
                                                          tCAAACTGGTTGCCGACCCTCTACTCTTCCAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACcg                              >  1:2044863/1‑93 (MQ=255)
                                                           cAAACTGGTTGCCGACCCTGTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACcgc                             <  1:1190251/93‑1 (MQ=255)
                                                            aaaCTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACcgcg                            >  2:4014526/1‑93 (MQ=255)
                                                            aaaCTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACcgcg                            >  2:1835584/1‑93 (MQ=255)
                                                             aaCTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCgg                           >  2:2309668/1‑93 (MQ=255)
                                                             aaCTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCgg                           >  1:2212310/1‑93 (MQ=255)
                                                              aCTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGa                          >  2:4167668/1‑93 (MQ=255)
                                                               cTGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGat                         >  1:612375/1‑93 (MQ=255)
                                                                tGGTTGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGatt                        <  1:3961585/93‑1 (MQ=255)
                                                                  ggtGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATACGTTATTTTGAAAACCGCGGattat                      >  1:4529112/3‑93 (MQ=255)
                                                                   ttGCCGTCCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGattatt                     <  1:2808979/93‑1 (MQ=255)
                                                                   ttGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGattatt                     >  2:1630402/1‑93 (MQ=255)
                                                                    tGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTAttt                    <  1:1561295/93‑1 (MQ=255)
                                                                    tGCCGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCAGTTATTTTGAAAACCGCGGATTAttt                    >  1:4058306/1‑93 (MQ=255)
                                                                    tGCCGACCCTCTACTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTAttt                    <  1:1320329/93‑1 (MQ=255)
                                                                      ccGACCCTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgc                  <  2:2814403/93‑1 (MQ=255)
                                                                          cccTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcat              >  2:1070401/1‑93 (MQ=255)
                                                                           ccTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcatt             <  1:2857503/93‑1 (MQ=255)
                                                                           ccTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcatt             >  2:2731716/1‑93 (MQ=255)
                                                                           ccTCTACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcatt             <  2:4302867/93‑1 (MQ=255)
                                                                            ctctACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattg            <  2:1963339/93‑1 (MQ=255)
                                                                            ctctACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattg            >  1:4362620/1‑93 (MQ=255)
                                                                             tctACTCTTTGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattgg           >  2:1246955/1‑92 (MQ=255)
                                                                             tctACTCTTCGAGCGATTTGATGAAATCGCCTTCAATTCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattgc           >  1:4125675/1‑93 (MQ=255)
                                                                             tctACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattgc           >  2:1589673/1‑93 (MQ=255)
                                                                             tctACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattgc           <  2:2560126/93‑1 (MQ=255)
                                                                              ctACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattgcg          >  1:1858659/1‑93 (MQ=255)
                                                                               tACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattgcgc         >  1:1023639/1‑93 (MQ=255)
                                                                               tACTCTTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATttgcgcattgcgc         <  2:3312907/93‑1 (MQ=255)
                                                                                  tctTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAAc      >  1:3478563/1‑93 (MQ=255)
                                                                                   ctTCGAGCGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAAcc     >  1:2017686/1‑93 (MQ=255)
                                                                                         |                                                                                         
GATGAACTCACAGGCTGTGTCTTCGGCGCGATCCCCCCTTTCAGCTTCCATCCAAAACTCAAACTGGTTGCCGACCCTCTACTCTTCGAACGATTTGATGAAATCGCCTTCAATGCAGGCATGCTGGATAAATCCGTTATTTTGAAAACCGCGGATTATTTGCGCATTGCGCAACCAGA  >  NC_000913/1871895‑1872073

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: