Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,663,474 A→T 15.6% A11A (GCA→GCT suhB → inositol monophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,663,4740AT15.6% 494.8 / 71.9 192A11A (GCA→GCTsuhBinositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base A (74/88);  new base T (13/17);  total (87/105)
Fisher's exact test for biased strand distribution p-value = 8.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 211 total aligned reads displayed.

ATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAG  >  NC_000913/2663377‑2663571
                                                                                                 |                                                                                                 
aTGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagc                                                                                                 <  2:542916/100‑1 (MQ=255)
    cGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGCGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgc                                                                                             >  2:961532/1‑100 (MQ=255)
    cGATCTCGCGGCTATGCTCTGCGCCGTTCTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgc                                                                                             <  1:1008689/100‑1 (MQ=255)
       tctcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAg                                                                                          <  1:3474230/100‑1 (MQ=255)
         tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                        <  1:872462/100‑1 (MQ=255)
         tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                        <  1:1763186/100‑1 (MQ=255)
         tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                        <  1:243229/100‑1 (MQ=255)
          cGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGcn                                                                                       <  2:3149796/100‑2 (MQ=255)
          cGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCg                                                                                       >  1:3749934/1‑100 (MQ=255)
          cGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCg                                                                                       >  2:3636387/1‑100 (MQ=255)
           gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                      <  2:3059658/100‑1 (MQ=255)
             tgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatg                                                                                    <  2:194322/100‑4 (MQ=38)
              gctATCCTCTTCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGACAGACCGATGCATCCTATGCTGAACATCGCCGTCCGCGCAGCGCGCAAGGCg                                                                                       <  2:829780/96‑1 (MQ=255)
              gctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGta                                                                                   >  2:2883470/1‑100 (MQ=255)
                tataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaat                                                                                 >  1:851060/1‑100 (MQ=255)
                 ataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatt                                                                                >  1:942293/1‑100 (MQ=255)
                  taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgatgttcag                                                                               <  1:3712098/100‑9 (MQ=37)
                  taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaattt                                                                               <  1:2581177/100‑1 (MQ=255)
                    ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaa                                                                             >  1:4054361/1‑100 (MQ=255)
                       tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                          <  2:2842215/100‑1 (MQ=255)
                       tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTg                                                                          <  2:3479804/100‑1 (MQ=255)
                         cgcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGcc                                                                        >  2:1685002/1‑100 (MQ=255)
                          gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCa                                                                       >  1:2877681/1‑100 (MQ=255)
                          gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCa                                                                       <  1:3301389/100‑1 (MQ=255)
                           ccGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaa                                                                      >  2:2975144/1‑100 (MQ=255)
                            cGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaa                                                                     >  1:2449155/1‑100 (MQ=255)
                                  cccGTTGTTGAACATCCACGGAGAGAGGCCGATGCATCGGATGCTGAACGTCGCCGGGCGCGCAGCGCGCAAGGCGGGGAATTTAATTGCCAAAAACTAt                                                               <  2:1875145/100‑1 (MQ=255)
                                   ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATg                                                              <  1:733276/100‑1 (MQ=255)
                                   ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATg                                                              <  1:2855796/100‑1 (MQ=255)
                                   ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATg                                                              >  1:1640925/1‑100 (MQ=255)
                                       tCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAc                                                          <  2:3737389/100‑1 (MQ=255)
                                        cTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcc                                                         >  1:613270/1‑100 (MQ=255)
                                         tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGCAAccc                                                        >  2:3349294/1‑100 (MQ=255)
                                          ttAACGTCCAGTGAGAGAGACCGATGCATCAGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcccc                                                       <  2:1688114/100‑1 (MQ=255)
                                          ttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcccc                                                       >  2:38324/1‑100 (MQ=255)
                                           tAACGTCTAGTGAGAGAGACCGATGCATCCGATTGTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAGACCCCg                                                      <  2:1362634/100‑1 (MQ=255)
                                           tAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCg                                                      >  1:1535523/1‑100 (MQ=255)
                                           tAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCg                                                      <  1:2607672/100‑1 (MQ=255)
                                           tAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCg                                                      <  1:4072304/100‑1 (MQ=255)
                                            aaCATCCGTTGAGAGAGACCGATACATCCTATGCTTAACATCGTCGTGCGCGCAGCGCGCAAGTCTAGTACTTGAATTGCCAAAAACTATGAAAc                                                          <  2:3966879/95‑1 (MQ=255)
                                             aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGa                                                    <  1:3067280/100‑1 (MQ=255)
                                              cATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGAc                                                   >  2:1204703/1‑100 (MQ=255)
                                              cATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGAc                                                   >  1:107749/1‑100 (MQ=255)
                                              cATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGAc                                                   <  1:3108899/100‑1 (MQ=255)
                                                tCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGAc                                                   >  1:1467626/1‑98 (MQ=255)
                                                 ccAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCt                                                <  1:2710489/100‑1 (MQ=255)
                                                  cAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                               >  2:1770296/1‑100 (MQ=255)
                                                   aGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGt                                              >  1:2372363/1‑100 (MQ=255)
                                                    gTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTa                                             <  2:3240025/100‑1 (MQ=255)
                                                    gTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTa                                             <  1:2625272/100‑1 (MQ=255)
                                                     tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAg                                            >  2:2140285/1‑100 (MQ=255)
                                                        gagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAg                                         >  1:4196320/1‑100 (MQ=255)
                                                          gagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCg                                       <  1:2164733/100‑1 (MQ=255)
                                                            gaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                     >  1:2587693/1‑100 (MQ=255)
                                                            gaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                     <  2:74613/100‑1 (MQ=255)
                                                            gaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGGAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                     >  1:2928799/1‑100 (MQ=255)
                                                             aCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGc                                    <  2:1279942/100‑1 (MQ=255)
                                                             aCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGc                                    <  1:248793/100‑1 (MQ=255)
                                                              ccGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGcc                                   <  1:2048980/100‑1 (MQ=255)
                                                              ccGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGcc                                   <  2:757240/100‑1 (MQ=255)
                                                              ccGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGcc                                   >  2:3225057/1‑100 (MQ=255)
                                                                    cATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAg                             >  1:477616/1‑100 (MQ=255)
                                                                    cATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAg                             <  1:1858014/100‑1 (MQ=255)
                                                                     aTCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCAAGAAAgg                            >  1:3220542/1‑100 (MQ=255)
                                                                      tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAgg                            <  2:2389828/99‑1 (MQ=255)
                                                                      tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                           >  1:1727236/1‑100 (MQ=255)
                                                                      tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                           <  2:2215088/100‑1 (MQ=255)
                                                                      tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                           <  2:2727540/100‑1 (MQ=255)
                                                                       ccGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                          >  1:2792735/1‑100 (MQ=255)
                                                                        caatacTGAACATCACCGTGCGCGCAGCGCGCAATGCGGGTAAAAAAAAAAAAAAAAACTATGAAACCCCGGACGCTGTAGGAGCGAGCCAGAAAGGCAg                         <  2:1623468/95‑1 (MQ=255)
                                                                        cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                         >  2:878653/1‑100 (MQ=255)
                                                                         gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                         <  1:129850/99‑1 (MQ=255)
                                                                         gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGt                        >  1:280101/1‑100 (MQ=255)
                                                                         gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGt                        >  1:2539040/1‑100 (MQ=255)
                                                                         gATGCTGAACACCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAg                         <  1:511746/99‑1 (MQ=255)
                                                                          aTGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCGAGAAAGGCAGTa                       <  1:2751137/100‑1 (MQ=255)
                                                                               gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAt                  >  2:2502105/1‑100 (MQ=255)
                                                                                aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAtt                 >  2:3647395/1‑100 (MQ=255)
                                                                                aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAtt                 <  1:4133541/100‑1 (MQ=255)
                                                                                aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAtt                 >  2:1869338/1‑100 (MQ=255)
                                                                                 tacccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAttt                >  2:48437/5‑100 (MQ=38)
                                                                                  cATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTc               >  2:4184426/1‑100 (MQ=255)
                                                                                   cccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCg              >  2:130356/3‑100 (MQ=39)
                                                                                    tCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt             >  1:3788597/1‑100 (MQ=255)
                                                                                    ncGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt             >  1:3573945/2‑100 (MQ=255)
                                                                                    ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt             >  1:1708999/2‑100 (MQ=39)
                                                                                    ccGCCTTGCGCGCTGCGCGAACGGCGGGTAATTTAATTGCGAAAAACTATGAAACCCCGGATGCTGTAGAAGGCAGACAGGAAGGCAGTAACGATTTCGt             >  2:2438105/2‑100 (MQ=255)
                                                                                     cGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTg            >  2:2080630/1‑100 (MQ=255)
                                                                                     cGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTg            <  1:3522656/100‑1 (MQ=255)
                                                                                      gCCTTGCGCGCTGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa           >  2:1157124/1‑100 (MQ=255)
                                                                                        cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAcc         <  1:2900934/100‑1 (MQ=255)
                                                                                        cGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAcc         <  1:625510/100‑1 (MQ=255)
                                                                                         gTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGATGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa        >  2:2676353/1‑100 (MQ=255)
                                                                                         gTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCa        <  2:1752041/100‑1 (MQ=255)
                                                                                            cgcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACg     <  2:4054334/100‑1 (MQ=255)
                                                                                            cgcgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACg     >  1:2685992/1‑100 (MQ=255)
                                                                                              cgcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa   <  1:657000/100‑1 (MQ=255)
                                                                                               gcagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAg  <  1:818688/100‑1 (MQ=255)
                                                                                                cagcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAg  <  2:4148060/99‑1 (MQ=255)
                                                                                                 agcgcgcaAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGCCCAACGt    >  2:3609515/1‑96 (MQ=255)
                                                                                                 |                                                                                                 
ATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAG  >  NC_000913/2663377‑2663571

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