Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,449,800 T→C T623A (ACC→GCC)  maoC ← putative aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,8000TC98.8% 520.1 / ‑2.9 163T623A (ACC→GCC) maoCputative aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (72/89);  total (72/91)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

Only 100 of 168 total aligned reads displayed.

AGCCCATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGTATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGACAAACAGAC  >  NC_000913/1449711‑1449889
                                                                                         |                                                                                         
aGCCCATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATc                                                                                        <  2:3116197/93‑1 (MQ=255)
aGCCCATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATc                                                                                        >  1:1274598/1‑93 (MQ=255)
aGCCCATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATc                                                                                        >  2:3872331/1‑93 (MQ=255)
 gCCCATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCg                                                                                       >  2:2067921/1‑93 (MQ=255)
 gCCCATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCg                                                                                       >  1:2121240/1‑93 (MQ=255)
    cATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGTATCGCCt                                                                                    <  2:3990402/93‑1 (MQ=255)
    cATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCt                                                                                    <  2:1087093/93‑1 (MQ=255)
     aTTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTg                                                                                   <  2:591332/93‑1 (MQ=255)
     aTTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTg                                                                                   <  1:3899513/93‑1 (MQ=255)
     aTTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTg                                                                                   <  2:4209614/93‑1 (MQ=255)
      ttCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTgg                                                                                  >  1:4601621/1‑93 (MQ=255)
      ttCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTgg                                                                                  <  1:4240381/93‑1 (MQ=255)
      ttCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTgg                                                                                  <  2:3731289/93‑1 (MQ=255)
      ttCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTgg                                                                                  >  2:3825372/1‑93 (MQ=255)
      ttCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTgg                                                                                  >  1:1472754/1‑93 (MQ=255)
         cgccacACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCtt                                                                               <  1:3287348/91‑1 (MQ=255)
         caccccACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCtt                                                                               <  2:2761775/93‑1 (MQ=255)
         caccacACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCtt                                                                               <  1:4737372/93‑1 (MQ=255)
          accacACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCttt                                                                              >  2:1146310/1‑93 (MQ=255)
          accacACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCttt                                                                              >  2:725779/1‑93 (MQ=255)
          accacACCTGTTGGTTTTTCTTCTGCGCTACGCTGTGTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGTCATCACCTGCCttt                                                                              >  1:4283163/1‑93 (MQ=255)
           ccacACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTa                                                                             <  2:3805715/93‑1 (MQ=255)
            cacaCCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTag                                                                            >  1:3355936/1‑92 (MQ=255)
            cacaCCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTAc                                                                            <  2:4713721/93‑1 (MQ=255)
            cacaCCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTAc                                                                            >  1:4061235/1‑93 (MQ=255)
              cacCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACgg                                                                          >  2:460316/1‑93 (MQ=255)
              cacCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACgg                                                                          >  1:3659754/1‑93 (MQ=255)
                ccTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTCACGGGt                                                                        <  1:3146400/93‑1 (MQ=255)
                 cTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGtt                                                                       <  1:3771873/93‑1 (MQ=255)
                 cTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGtt                                                                       <  1:3002143/93‑1 (MQ=255)
                 cTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGtt                                                                       >  1:2924357/1‑93 (MQ=255)
                 cTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGtt                                                                       <  2:1873113/93‑1 (MQ=255)
                 cTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGtt                                                                       <  2:1811867/93‑1 (MQ=255)
                 cTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGtt                                                                       <  1:3570739/93‑1 (MQ=255)
                  tGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTc                                                                      >  2:2668638/1‑93 (MQ=255)
                   gTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCg                                                                     >  2:4530014/1‑93 (MQ=255)
                    ttGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGa                                                                    <  2:92849/93‑1 (MQ=255)
                     tGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGAt                                                                   >  2:2343038/1‑93 (MQ=255)
                        tttttCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATaaa                                                                <  1:4703988/93‑1 (MQ=255)
                          tttCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAacac                                                              >  1:1004088/1‑93 (MQ=255)
                          tttCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAAc                                                              <  1:3566505/93‑1 (MQ=255)
                            tcttctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGc                                                            <  1:1005152/93‑1 (MQ=255)
                            tcttctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGc                                                            <  1:2135801/93‑1 (MQ=255)
                            tcttctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGc                                                            >  1:958541/1‑93 (MQ=255)
                            tcttcgCCGCTACGCTGTTTTTTCAGCGGCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGc                                                            <  1:381766/93‑1 (MQ=255)
                             cttctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCa                                                           >  1:1566134/1‑93 (MQ=255)
                              ttctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCaa                                                          >  1:2230294/1‑93 (MQ=255)
                              ttctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCaa                                                          <  1:894795/93‑1 (MQ=255)
                              ttctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCaa                                                          <  1:4683249/93‑1 (MQ=255)
                                ctGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATACAACGCAAGc                                                        >  1:4125015/1‑93 (MQ=255)
                                  gcgcTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCtt                                                      <  1:2200027/93‑1 (MQ=255)
                                   cgcTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCttt                                                     <  2:788876/93‑1 (MQ=255)
                                    gcTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTc                                                    >  1:2026315/1‑93 (MQ=255)
                                     cTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTcc                                                   >  2:190831/1‑93 (MQ=255)
                                     cTACGCTGTTTTTTCAGCGTCTTGCGCTTAAAAGTCAGACGCACCTGGCTGACATCGCCTGGCATTCCGGGTTAGATAAAACGCAAACTTTcc                                                   >  2:2210927/1‑93 (MQ=255)
                                      tACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCa                                                  <  2:3298799/93‑1 (MQ=255)
                                       aCGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAg                                                 >  2:1225725/1‑93 (MQ=255)
                                       aCGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAg                                                 <  1:358646/93‑1 (MQ=255)
                                        cGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGc                                                <  1:4508650/93‑1 (MQ=255)
                                         gCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGcc                                               >  1:487419/1‑93 (MQ=255)
                                           tgtttgttCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCg                                             <  1:3634045/93‑1 (MQ=255)
                                             ttttttCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTa                                           <  1:3436227/93‑1 (MQ=255)
                                                 ttCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTa                                       <  2:1471/93‑1 (MQ=255)
                                                 ttCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTa                                       <  1:1808029/93‑1 (MQ=255)
                                                 ttCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTa                                       >  2:3072693/1‑93 (MQ=255)
                                                  tCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAg                                      >  1:3374456/1‑93 (MQ=255)
                                                     gCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCaa                                   <  2:3922691/93‑1 (MQ=255)
                                                     gCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCaa                                   >  2:1640682/1‑93 (MQ=255)
                                                     gCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCaa                                   <  2:3013879/93‑1 (MQ=255)
                                                        tCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTGCGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATga                                >  2:1682319/1‑93 (MQ=255)
                                                        tCTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATga                                >  1:3621645/1‑93 (MQ=255)
                                                         cTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgac                               <  1:1086911/93‑1 (MQ=255)
                                                         cTTGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgac                               <  1:4506458/93‑1 (MQ=255)
                                                          ttGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgacc                              >  1:3124876/1‑93 (MQ=255)
                                                           tGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccg                             <  1:60445/93‑1 (MQ=255)
                                                           tGCGCTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTAGCAATgaccgg                            >  1:1817793/1‑93 (MQ=255)
                                                             cgcTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggc                           >  1:2643360/1‑92 (MQ=255)
                                                             cgcTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccgga                           <  2:3196474/93‑1 (MQ=255)
                                                             cgcTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccgga                           <  1:1472447/93‑1 (MQ=255)
                                                              gcTTACAGGTGAGACTCCCCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggac                          <  1:647227/93‑1 (MQ=255)
                                                              gcTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggac                          <  1:3995490/93‑1 (MQ=255)
                                                              gcTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggac                          >  2:3641791/1‑93 (MQ=255)
                                                              gcTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggac                          >  2:1217095/1‑93 (MQ=255)
                                                               cTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggacc                         <  1:2948002/93‑1 (MQ=255)
                                                               cTTACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggacc                         <  1:416958/93‑1 (MQ=255)
                                                                ttACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATgaccggaccg                        >  2:4424517/1‑93 (MQ=255)
                                                                 tACAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGa                       <  2:1885233/93‑1 (MQ=255)
                                                                  aCAGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTTGTTAGCAATGACCGGACCGac                      <  2:4671348/93‑1 (MQ=255)
                                                                    aGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGacac                    >  1:4440061/1‑93 (MQ=255)
                                                                    aGGTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGacac                    <  2:2054854/93‑1 (MQ=255)
                                                                     ggTGAGACGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACAcc                   >  1:4366266/1‑93 (MQ=255)
                                                                          gaCGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGGCCGACACCGGCAt              >  1:2808450/1‑93 (MQ=255)
                                                                           aCGCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATc             <  1:3758576/93‑1 (MQ=255)
                                                                             gCACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGa           >  1:1021726/1‑93 (MQ=255)
                                                                              cACCTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGAc          <  2:3755746/93‑1 (MQ=255)
                                                                                 cTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGACaaa       >  2:3831247/1‑93 (MQ=255)
                                                                                 cTGGATGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGACaaa       >  1:1130867/1‑93 (MQ=255)
                                                                                     aTGGCATCGCCTGGCTTTACGGGCTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGACAAACAGa   <  2:3395357/93‑1 (MQ=255)
                                                                                      tGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGACAAACAGAc  >  2:4312929/1‑93 (MQ=255)
                                                                                      tGGCATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGACAAACAGAc  <  1:3390736/93‑1 (MQ=255)
                                                                                         |                                                                                         
AGCCCATTCCACCACACCTGTTGGTTTTTCTTCTGCGCTACGCTGTTTTTTCAGCGTCTTGCGCTTACAGGTGAGACGCACCTGGATGGTATCGCCTGGCTTTACGGGTTCGATAAAACGCAAGCTTTCCAGCCCGTAGTTAGCAATGACCGGACCGACACCGGCATCGACAAACAGAC  >  NC_000913/1449711‑1449889

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