Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,559:1 +A 55.6% pseudogene (18/663 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,5591.A55.6% 157.7 / 140.5 81pseudogene (18/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base . (21/15);  new base A (28/17);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACC  >  NC_000913/3951460‑3951654
                                                                                                      |                                                                                                
aaCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgc                                                                                                   <  2:173025/100‑1 (MQ=255)
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgca                                                                                                  >  2:1406785/1‑99 (MQ=255)
  cAAGCTGCGTCAGGTACAGGTGGTATTACAATGTGATTTAAATGCTCTGTTACCAGCATTTCATCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgcag                                                                                                 <  2:1900709/100‑3 (MQ=255)
  cAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgcag                                                                                                 <  1:664951/100‑3 (MQ=255)
  cAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgcag                                                                                                 <  1:2028301/100‑3 (MQ=255)
  cAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgcag                                                                                                 <  2:1488099/100‑3 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgcagcg                                                                                               <  1:695651/100‑5 (MQ=255)
    aGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTgcagcg                                                                                               >  2:346691/1‑96 (MQ=255)
     gCTGCGTCAGGCACATGTGGCATTACAGGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATATCGATGACTGGCAGCAACACTgcg                                                                                                 >  2:161405/1‑98 (MQ=255)
           tCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCg                                                                                       >  2:275032/1‑100 (MQ=255)
           tCAGGCACATGTGCCATTACAAGGTGATTTACATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGc                                                                                        >  1:387666/1‑100 (MQ=255)
            cAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGt                                                                                      >  2:409644/1‑100 (MQ=255)
            cAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGt                                                                                      >  2:2161595/1‑100 (MQ=255)
             aGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtg                                                                                     <  2:1663311/100‑1 (MQ=255)
                  caTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtgatga                                                                                 <  1:971708/100‑1 (MQ=255)
                   aTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGaa                                                                                >  1:1945485/1‑100 (MQ=255)
                   aTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGaa                                                                                <  2:1851107/100‑1 (MQ=255)
                     gtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACa                                                                              >  2:1328704/1‑100 (MQ=255)
                     gtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACa                                                                              >  1:295539/1‑100 (MQ=255)
                       ggCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAt                                                                             <  2:796429/99‑1 (MQ=255)
                        gCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCCATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTc                                                                           >  2:52149/1‑100 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTcc                                                                          <  1:1737184/100‑1 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCt                                                                         >  1:866096/1‑100 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCt                                                                         >  2:873588/1‑100 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTgg                                                                       >  2:2038444/1‑100 (MQ=255)
                              cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                     >  2:847830/1‑100 (MQ=255)
                              cAAGGTGAGTTAAATGCTCTGTTACCAGCATTACAGCAGCCGGTAA‑‑ATCAATGACTGGCAGCAACACGGCGCGCAGCTGCGTGATGAACATTCCTGGCGt                                                                    <  2:1068834/100‑1 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTa                                                                  >  2:2236311/1‑100 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTCGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                 >  2:734366/1‑100 (MQ=255)
                                    gATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                               >  2:1774849/1‑100 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAc                                                              <  2:1945472/100‑1 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                             >  2:823123/1‑100 (MQ=255)
                                      tttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                             <  1:1869315/100‑1 (MQ=255)
                                       ttAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCa                                                            <  1:1623963/100‑1 (MQ=255)
                                        tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCAt                                                           >  1:1489877/1‑100 (MQ=255)
                                        tAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATc                                                          >  1:56091/1‑100 (MQ=255)
                                          aaTGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATcc                                                         <  2:130862/100‑1 (MQ=255)
                                           aTGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATccc                                                        >  2:1443385/1‑100 (MQ=255)
                                           aTGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCg                                                       <  1:1647535/100‑1 (MQ=255)
                                            tGCTCTGTTACCAGCAGTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATccc                                                        <  2:805958/99‑1 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                      <  2:627719/100‑1 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGa                                                   >  1:1824229/1‑100 (MQ=255)
                                                    tACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCt                                               >  1:2077657/1‑100 (MQ=255)
                                                     aCCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTa                                              >  1:294737/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                             >  1:88764/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                             <  2:188877/100‑1 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                             >  1:1671696/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                             >  1:958853/1‑100 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTAt                                             >  1:1188663/1‑99 (MQ=255)
                                                      ccAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATc                                            >  2:1085631/1‑100 (MQ=255)
                                                        aGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTa                                          <  1:173933/100‑1 (MQ=255)
                                                          cATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                         >  1:578209/1‑100 (MQ=255)
                                                          cATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                         <  2:323339/100‑1 (MQ=255)
                                                           aTTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgc                                       <  2:866095/100‑1 (MQ=255)
                                                            ttACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgcg                                      <  1:1622445/100‑1 (MQ=255)
                                                            ttACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcgcg                                      >  2:491973/1‑100 (MQ=255)
                                                              aCAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACgcgc                                     >  2:11757/1‑100 (MQ=255)
                                                              aCAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                    >  1:1160713/1‑100 (MQ=255)
                                                              aCAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                    <  1:409643/100‑1 (MQ=255)
                                                               cagcagCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                    >  2:504303/1‑100 (MQ=255)
                                                                 gcagcCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                  >  2:536984/1‑100 (MQ=255)
                                                                   agcCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttg                                >  2:530354/1‑100 (MQ=255)
                                                                     ccGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                             >  2:483921/1‑100 (MQ=255)
                                                                     ccGTCAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                             <  1:427658/100‑1 (MQ=255)
                                                                      ngTTAA‑‑ATCAATGACTGGCAGCACCACTGCAGCGCAGCTGGGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCCCCgttgtcg                             >  2:203717/2‑98 (MQ=255)
                                                                      cGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTCCGCGCCgttgttg                             >  1:1300695/1‑100 (MQ=255)
                                                                      cGTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                             <  1:121003/100‑1 (MQ=255)
                                                                       gTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                            <  2:685795/100‑1 (MQ=255)
                                                                       gTTAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                            <  1:1415076/100‑1 (MQ=255)
                                                                       gTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                           >  2:2197381/1‑100 (MQ=255)
                                                                        ttAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                          <  1:791423/100‑1 (MQ=255)
                                                                        ttAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                          <  1:1076612/100‑1 (MQ=255)
                                                                        ttAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                          <  2:56091/100‑1 (MQ=255)
                                                                         tAA‑‑ATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                          <  1:86281/100‑1 (MQ=255)
                                                                         tAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                         >  2:1413823/1‑100 (MQ=255)
                                                                              aTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                     <  1:646135/100‑1 (MQ=255)
                                                                              aTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                     >  1:1792475/1‑100 (MQ=255)
                                                                               tCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                     >  1:475360/1‑100 (MQ=255)
                                                                               tCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCCATCTACGCGCCGTTGTTGTTAAaac                     <  1:373477/100‑1 (MQ=255)
                                                                                cAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                   >  1:783889/1‑100 (MQ=255)
                                                                                     aCTGGTAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGTCGTTACGACCATCCCGGTGACGCTATCTACGCTCCGTTGTTTTaacaacaacggcg              <  2:1062785/100‑5 (MQ=255)
                                                                                     aCTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGt               >  2:2247213/1‑100 (MQ=255)
                                                                                      cTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc              >  1:348686/1‑100 (MQ=255)
                                                                                      cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg             >  2:1405182/1‑100 (MQ=255)
                                                                                      cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg             >  2:20112/1‑100 (MQ=255)
                                                                                           agcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg        <  2:1641203/100‑1 (MQ=255)
                                                                                            gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGt       >  1:1935843/1‑100 (MQ=255)
                                                                                                cacTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGATGTTAAAACAACTGTCGAATCGt       >  2:1725264/1‑97 (MQ=255)
                                                                                                 acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc  <  1:1254366/100‑1 (MQ=255)
                                                                                                 acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc  >  2:1729535/1‑100 (MQ=255)
                                                                                                      |                                                                                                
AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAA‑‑ATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACC  >  NC_000913/3951460‑3951654

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