Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,845,598 T→C G545G (GGT→GGC adeD → adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,845,5980TC100.0% 1244.2 / NA 294G545G (GGT→GGCadeDadenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (130/164);  total (130/164)

Only 100 of 302 total aligned reads displayed.

CATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGA  >  NC_000913/3845500‑3845695
                                                                                                  |                                                                                                 
cATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGccc                                                                                                 >  1:5678538/1‑101 (MQ=255)
cATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGccc                                                                                                 >  2:5678537/1‑101 (MQ=255)
cATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGccc                                                                                                 <  2:5678107/101‑1 (MQ=255)
   cTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGtt                                                                                              <  1:5678003/101‑1 (MQ=255)
    tGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTa                                                                                             >  1:5678541/1‑101 (MQ=255)
      ccGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTAcc                                                                                           <  2:5678275/101‑1 (MQ=255)
       cGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTAccc                                                                                          >  1:5678546/1‑101 (MQ=255)
       cGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTAccc                                                                                          <  2:5678134/101‑1 (MQ=255)
       cGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTAccc                                                                                          <  2:5678053/101‑1 (MQ=255)
        gTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCg                                                                                         >  2:5678548/1‑101 (MQ=255)
        gTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCg                                                                                         <  1:5678090/101‑1 (MQ=255)
          tACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGAt                                                                                       <  2:5678191/101‑1 (MQ=255)
            cccATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGa                                                                                     <  1:5678018/101‑1 (MQ=255)
               aTTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGa                                                                                     <  1:5678307/98‑1 (MQ=255)
               aTTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGcc                                                                                  >  2:5678552/1‑101 (MQ=255)
                 tGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGt                                                                                >  2:5678557/1‑101 (MQ=255)
                   ccGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGttt                                                                              >  1:5678567/1‑101 (MQ=255)
                   ccGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGttt                                                                              >  1:5678563/1‑101 (MQ=255)
                   ccGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGttt                                                                              >  2:5678561/1‑101 (MQ=255)
                   ccGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGttt                                                                              <  1:5677758/101‑1 (MQ=255)
                     gggCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTAt                                                                            >  2:5678571/1‑101 (MQ=255)
                     gggCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTAt                                                                            >  1:5678570/1‑101 (MQ=255)
                       gCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTc                                                                          <  1:5678072/101‑1 (MQ=255)
                        cTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCa                                                                         <  1:5678195/101‑1 (MQ=255)
                         tgatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAg                                                                        <  1:5678137/101‑1 (MQ=255)
                         tgatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAg                                                                        >  1:5678575/1‑101 (MQ=255)
                         tgatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAg                                                                        >  2:5678578/1‑101 (MQ=255)
                         tgatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAg                                                                        <  1:5678043/101‑1 (MQ=255)
                            tgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATg                                                                     <  2:5678302/101‑1 (MQ=255)
                             gaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATgg                                                                    <  2:5678117/101‑1 (MQ=255)
                             gaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATgg                                                                    >  1:5678582/1‑101 (MQ=255)
                              aGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGc                                                                   <  1:5678325/101‑1 (MQ=255)
                                cACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGt                                                                 <  1:5678215/101‑1 (MQ=255)
                                 aCCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGtt                                                                >  2:5678586/1‑101 (MQ=255)
                                 aCCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGtt                                                                >  1:5678587/1‑101 (MQ=255)
                                  ccGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGttt                                                               <  2:5678044/101‑1 (MQ=255)
                                  ccGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGttt                                                               <  2:5678062/101‑1 (MQ=255)
                                    gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGtttct                                                             >  1:5678591/1‑101 (MQ=255)
                                    gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGtttct                                                             <  2:5678125/101‑1 (MQ=255)
                                     acacGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGtttctt                                                            <  2:5678415/101‑1 (MQ=255)
                                     acacGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGtttctt                                                            <  2:5678173/101‑1 (MQ=255)
                                     acacGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGtttctt                                                            >  2:5678593/1‑101 (MQ=255)
                                         ggCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTtctc                                                        <  2:5678159/101‑1 (MQ=255)
                                          gcgcAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTctct                                                       <  2:5678216/101‑1 (MQ=255)
                                            gcAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGc                                                     <  2:5678019/101‑1 (MQ=255)
                                            gcAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGc                                                     >  1:5678597/1‑101 (MQ=255)
                                            gcAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGc                                                     <  1:5677994/101‑1 (MQ=255)
                                             cAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGcc                                                    >  1:5678600/1‑101 (MQ=255)
                                              aGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCa                                                   <  2:5678208/101‑1 (MQ=255)
                                                tCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGt                                                 >  1:5678603/1‑101 (MQ=255)
                                                 cGCTGGCGGAACAAAGTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTg                                                      <  1:5678183/95‑1 (MQ=255)
                                                  gCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGa                                               >  1:5678606/1‑101 (MQ=255)
                                                    tGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATc                                             >  1:5678611/1‑101 (MQ=255)
                                                    tGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATc                                             >  2:5678610/1‑101 (MQ=255)
                                                     ggCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATcc                                            <  1:5678102/101‑1 (MQ=255)
                                                       cGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATcccc                                          <  2:5678030/101‑1 (MQ=255)
                                                       cGGAACAAATTGACCCCTTGAAACCCGCCGCCCGTTAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGt                                                                 <  2:5678225/78‑1 (MQ=255)
                                                         gAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCcgc                                        >  1:5678616/1‑101 (MQ=255)
                                                         gAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCcgc                                        >  2:5678619/1‑101 (MQ=255)
                                                          aaCAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCcgcg                                       >  1:5678620/1‑101 (MQ=255)
                                                           aCAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCgcgc                                      <  2:5678251/101‑1 (MQ=255)
                                                            cAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCt                                     <  1:5678268/101‑1 (MQ=255)
                                                              aaTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTaa                                   >  1:5678624/1‑101 (MQ=255)
                                                              aaTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTaa                                   <  1:5678286/101‑1 (MQ=255)
                                                                ttGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTaaaa                                 <  1:5678233/101‑1 (MQ=255)
                                                                  gACGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACt                               <  1:5678253/101‑1 (MQ=255)
                                                                    cGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTaa                             <  1:5678054/101‑1 (MQ=255)
                                                                    cGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTaa                             <  1:5678158/101‑1 (MQ=255)
                                                                    cGCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTaa                             >  1:5678628/1‑101 (MQ=255)
                                                                     gCCTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAAc                            >  2:5678629/1‑101 (MQ=255)
                                                                      ccTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAAcc                           >  1:5678632/1‑101 (MQ=255)
                                                                       cTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCa                          <  2:5678178/101‑1 (MQ=255)
                                                                       cTTGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCa                          <  2:5677981/101‑1 (MQ=255)
                                                                         tGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAAc                                <  2:5677918/93‑1 (MQ=255)
                                                                         tGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGt                        >  1:5678636/1‑101 (MQ=255)
                                                                         tGAAAGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGt                        >  2:5678638/1‑101 (MQ=255)
                                                                           aaaGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCa                      >  1:5678642/1‑101 (MQ=255)
                                                                           aaaGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCa                      <  2:5678241/101‑1 (MQ=255)
                                                                            aaGCCGCCGCCCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAg                     <  1:5678065/101‑1 (MQ=255)
                                                                               ccgccgccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAgggg                  >  2:5678647/1‑101 (MQ=255)
                                                                                 gccgccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCt                <  1:5678128/101‑1 (MQ=255)
                                                                                 gccgccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCt                <  2:5678290/101‑1 (MQ=255)
                                                                                 gccgccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCt                >  1:5678650/1‑101 (MQ=255)
                                                                                  ccgccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTa               <  1:5678402/101‑1 (MQ=255)
                                                                                   cgccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAt              >  2:5678652/1‑101 (MQ=255)
                                                                                    gccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAtt             >  1:5678655/1‑101 (MQ=255)
                                                                                    gccCGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAtt             <  1:5678223/101‑1 (MQ=255)
                                                                                     cccGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAttt            <  1:5677960/101‑1 (MQ=255)
                                                                                     cccGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAttt            >  1:5678659/1‑101 (MQ=255)
                                                                                     cccGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAttt            <  1:5678383/101‑1 (MQ=255)
                                                                                      ccGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTg           <  2:5678369/101‑1 (MQ=255)
                                                                                       cGTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGa          <  2:5678340/101‑1 (MQ=255)
                                                                                        gTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGAt         <  2:5678151/101‑1 (MQ=255)
                                                                                        gTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGAt         >  1:5678662/1‑101 (MQ=255)
                                                                                        gTGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGAt         >  2:5678663/1‑101 (MQ=255)
                                                                                         tGAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATg        <  1:5678343/101‑1 (MQ=255)
                                                                                          gAATGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATgg       <  1:5678272/101‑1 (MQ=255)
                                                                                           aaTGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGc      <  2:5678472/101‑1 (MQ=255)
                                                                                            aTGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCg     >  1:5678668/1‑101 (MQ=255)
                                                                                             tGCGGCCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCga    >  2:5678669/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGA  >  NC_000913/3845500‑3845695

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