Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,900 G→A G113G (GGC→GGT yieO ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,9000GA99.6% 899.7 / ‑5.2 239G113G (GGC→GGTyieOputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (129/109);  total (130/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

Only 100 of 253 total aligned reads displayed.

CGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACT  >  NC_000913/3937810‑3937989
                                                                                          |                                                                                         
cgacATGGCGACAAAATTCAATACTGGAAGAAGTTCATGACGAGGATAAGCGCTCATTAAGGCCAGCCGAGCAACAGGCATCATCATTGcacc                                                                                         <  1:5351332/90‑1 (MQ=255)
cGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcacc                                                                                         >  1:1181342/1‑93 (MQ=255)
cGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcacc                                                                                         >  2:1283616/1‑93 (MQ=255)
  gCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcaccgc                                                                                       <  2:3930925/93‑1 (MQ=255)
    aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCt                                                                                     <  2:5076053/93‑1 (MQ=255)
    aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCt                                                                                     >  2:3823528/1‑93 (MQ=255)
    aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCt                                                                                     <  2:6272350/93‑1 (MQ=255)
      ggCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTAt                                                                                   >  2:4690427/1‑93 (MQ=255)
      ggCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTAt                                                                                   >  2:5617063/1‑93 (MQ=255)
       gCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATc                                                                                  >  2:2074022/1‑93 (MQ=255)
        cGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATcc                                                                                 <  1:1623243/93‑1 (MQ=255)
        cGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATcc                                                                                 >  2:2764021/1‑93 (MQ=255)
         gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATccc                                                                                >  1:660149/1‑93 (MQ=255)
         gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATccc                                                                                >  1:5900302/1‑93 (MQ=255)
          aCAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATcccc                                                                               <  2:2166362/93‑1 (MQ=255)
           cAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCt                                                                              <  2:6198993/93‑1 (MQ=255)
            aaaaTTCAATACTGGAAGAAGTTCATTACGAGGAGAAGCGCGCAGTAAGGCCAGCGGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTg                                                                             <  2:2490310/93‑1 (MQ=255)
             aaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGa                                                                            >  2:5002202/1‑93 (MQ=255)
              aaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGaa                                                                           >  1:2534218/1‑93 (MQ=255)
               aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAAt                                                                          <  1:6096366/93‑1 (MQ=255)
               aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAAt                                                                          <  2:2912564/93‑1 (MQ=255)
                ttCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATa                                                                         >  2:274723/1‑93 (MQ=255)
                 tCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATaa                                                                        >  1:1341527/1‑93 (MQ=255)
                  cAATACTGGAAGCAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAc                                                                       >  1:6022861/1‑93 (MQ=255)
                  cAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAc                                                                       >  2:6436136/1‑93 (MQ=255)
                  cAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAc                                                                       >  1:2667561/1‑93 (MQ=255)
                    aTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAccc                                                                     >  2:1396738/1‑93 (MQ=255)
                     tACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCg                                                                    <  1:5523317/93‑1 (MQ=255)
                      aCTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCgg                                                                   >  2:4799422/1‑93 (MQ=255)
                      aCTGGAAGAAGTTCATTACGAGGATAAGCGCGAAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCgg                                                                   <  1:6175108/93‑1 (MQ=255)
                       cTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGa                                                                  <  2:3258957/93‑1 (MQ=255)
                        tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGaa                                                                 <  2:6143886/93‑1 (MQ=255)
                          gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGa                                                               >  1:2713041/1‑93 (MQ=255)
                          gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGa                                                               >  1:728713/1‑93 (MQ=255)
                          gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGa                                                               <  1:3913654/93‑1 (MQ=255)
                             gaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACaa                                                            >  1:4048700/1‑93 (MQ=255)
                               atttCATTACGAGGATAAGCTCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAAcc                                                          <  1:1456762/91‑1 (MQ=255)
                               agTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAAcc                                                          <  2:1197593/93‑1 (MQ=255)
                               agTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAAcc                                                          >  2:3168700/1‑93 (MQ=255)
                                gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCa                                                         <  2:3299826/93‑1 (MQ=255)
                                gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCa                                                         <  2:3463194/93‑1 (MQ=255)
                                 ttCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAg                                                        >  1:805420/1‑93 (MQ=255)
                                  tCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGc                                                       <  2:4038712/93‑1 (MQ=255)
                                   cATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATGCTCATTGCACCGGCTATCCCCTGAATAACCCGCAAGACAACCAGCt                                                      >  1:4371259/1‑93 (MQ=255)
                                   cATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCt                                                      >  1:3045426/1‑93 (MQ=255)
                                   cATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCt                                                      >  1:4858042/1‑93 (MQ=255)
                                   cATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCt                                                      >  2:5560923/1‑93 (MQ=255)
                                    aTTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtg                                                     >  2:5332576/1‑93 (MQ=255)
                                     ttACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtgt                                                    <  2:734240/93‑1 (MQ=255)
                                     ttACGAGGAGAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtgt                                                    <  2:3600571/93‑1 (MQ=255)
                                      tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtgtg                                                   <  1:1872373/93‑1 (MQ=255)
                                         gAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTa                                                <  2:451113/93‑1 (MQ=255)
                                         gAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTa                                                >  1:1935675/1‑93 (MQ=255)
                                           ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              >  1:6130165/1‑93 (MQ=255)
                                           ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              >  1:1114578/1‑93 (MQ=255)
                                           ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              <  2:433408/93‑1 (MQ=255)
                                           ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              <  1:5165620/93‑1 (MQ=255)
                                            gATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCg                                             >  2:4103486/1‑93 (MQ=255)
                                            gATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCg                                             <  1:2230645/93‑1 (MQ=255)
                                            gATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCg                                             >  2:5900483/1‑93 (MQ=255)
                                            gATAAGCGCGCAGCAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCg                                             >  1:1652898/1‑93 (MQ=255)
                                             aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGa                                            >  2:4490186/1‑93 (MQ=255)
                                             aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGa                                            <  2:2381446/93‑1 (MQ=255)
                                              tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGaa                                           >  1:1376258/1‑93 (MQ=255)
                                              tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGaa                                           >  1:6440893/1‑93 (MQ=255)
                                               aaGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTTTCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAt                                          <  1:3556608/93‑1 (MQ=255)
                                               aaGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAt                                          <  1:1770218/93‑1 (MQ=255)
                                               aaGCGCGCAGTAAGCCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAt                                          <  2:1777027/93‑1 (MQ=255)
                                                aGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAtt                                         >  1:1433946/1‑93 (MQ=255)
                                                 gcgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTa                                        <  2:5791378/93‑1 (MQ=255)
                                                   gcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGa                                      >  2:1480874/1‑93 (MQ=255)
                                                     gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGaaa                                    <  2:2539889/93‑1 (MQ=255)
                                                     gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGaaa                                    <  1:931290/93‑1 (MQ=255)
                                                      cAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAg                                   <  1:4540173/93‑1 (MQ=255)
                                                      cAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAg                                   >  1:504768/1‑93 (MQ=255)
                                                        gTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTg                                 >  1:338672/1‑93 (MQ=255)
                                                          aaGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcg                               <  1:5663197/93‑1 (MQ=255)
                                                           aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              <  1:4195324/93‑1 (MQ=255)
                                                           aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              <  2:5817254/93‑1 (MQ=255)
                                                           aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              >  2:5548762/1‑93 (MQ=255)
                                                             gCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAg                            >  2:370743/1‑93 (MQ=255)
                                                             gCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAg                            <  1:1224047/93‑1 (MQ=255)
                                                               cAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGc                          <  1:1025625/93‑1 (MQ=255)
                                                                aGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGca                         >  2:5416024/1‑92 (MQ=255)
                                                                 gCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCa                        <  2:5186674/93‑1 (MQ=255)
                                                                 gCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCa                        >  1:6266329/1‑93 (MQ=255)
                                                                   cGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                      >  2:5580990/1‑93 (MQ=255)
                                                                   cGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                      <  1:2401080/93‑1 (MQ=255)
                                                                    gAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCagag                     >  2:5664363/1‑93 (MQ=255)
                                                                    gAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCagag                     >  1:2902472/1‑93 (MQ=255)
                                                                    gAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCagag                     >  2:3013729/1‑93 (MQ=255)
                                                                       cAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAc                  <  1:4408870/93‑1 (MQ=255)
                                                                       cAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAc                  >  2:1862883/1‑93 (MQ=255)
                                                                            ggCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  2:892613/1‑93 (MQ=255)
                                                                             gCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtg            <  2:1671389/93‑1 (MQ=255)
                                                                                tcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGa         <  1:4257487/93‑1 (MQ=255)
                                                                                   tcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACa      >  2:5198917/1‑93 (MQ=255)
                                                                                   tcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACa      <  2:3872217/93‑1 (MQ=255)
                                                                                      ttGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGAc   >  1:3216531/1‑93 (MQ=255)
                                                                                       tGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACt  <  2:4284698/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACT  >  NC_000913/3937810‑3937989

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