Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,269,236 G→C A59G (GCT→GGT)  b3124 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,2360GC84.3% 702.2 / 110.9 249A59G (GCT→GGT) b3124hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/24);  new base C (99/111);  total (114/135)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

Only 100 of 255 total aligned reads displayed.

TTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCC  >  NC_000913/3269145‑3269328
                                                                                            |                                                                                            
 tCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGcc                                                                                            <  2:3785257/93‑1 (MQ=255)
   aaTAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACcccc                                                                                          <  2:7275269/93‑1 (MQ=255)
    aTAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAcccccc                                                                                         >  1:2450919/1‑93 (MQ=255)
    aTAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGccccc                                                                                         <  1:5203765/93‑1 (MQ=255)
    aTAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCAGCCCTGTAACCCAGGCGTGACGTTCAcccccc                                                                                         <  2:7113956/93‑1 (MQ=255)
      aaaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTg                                                                                       <  2:2015386/93‑1 (MQ=255)
      aaaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTg                                                                                       >  2:3466223/1‑93 (MQ=255)
       aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTgg                                                                                      <  1:7215464/93‑1 (MQ=255)
       aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTgg                                                                                      <  1:5651888/93‑1 (MQ=255)
        aCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTggg                                                                                     >  2:4045423/1‑93 (MQ=255)
        aCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTggg                                                                                     <  1:3786005/93‑1 (MQ=255)
         cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCAGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGt                                                                                    <  1:2580405/93‑1 (MQ=255)
          gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTg                                                                                   >  1:1281049/1‑93 (MQ=255)
          gcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTg                                                                                   >  2:7142165/1‑93 (MQ=255)
                ttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCa                                                                             >  1:7364330/1‑93 (MQ=255)
                 tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCaa                                                                            >  2:2799303/1‑93 (MQ=255)
                   ccATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAatc                                                                          >  2:6316451/1‑93 (MQ=255)
                   ccATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCAGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTCCAatc                                                                          <  1:4702742/93‑1 (MQ=255)
                    cATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAatca                                                                         <  1:6035960/93‑1 (MQ=255)
                     aTCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAatcat                                                                        <  2:2406812/93‑1 (MQ=255)
                      tCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAatcatc                                                                       >  2:1588616/1‑93 (MQ=255)
                      tCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAatcatc                                                                       >  1:3939087/1‑93 (MQ=255)
                       cGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCg                                                                      <  2:3189202/93‑1 (MQ=255)
                        gCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCAGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGc                                                                     >  1:4465826/1‑93 (MQ=255)
                          cGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCtt                                                                   <  2:2300016/93‑1 (MQ=255)
                           ggAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTc                                                                  >  2:5024010/1‑93 (MQ=255)
                            gagaTCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCGCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTcc                                                                 >  1:6221423/1‑93 (MQ=255)
                            gagaTCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTcc                                                                 <  2:750334/93‑1 (MQ=255)
                            gagaTCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTcc                                                                 >  1:3289290/1‑93 (MQ=255)
                             agaTCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCa                                                                <  2:2248376/93‑1 (MQ=255)
                              gaTCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCAc                                                               <  2:390969/93‑1 (MQ=255)
                              gaTCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCAc                                                               >  2:45484/1‑93 (MQ=255)
                                tCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCg                                                             <  1:5496211/93‑1 (MQ=255)
                                   cccAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTc                                                          >  1:4934003/1‑93 (MQ=255)
                                   cccAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTc                                                          >  1:2210956/1‑93 (MQ=255)
                                    ccAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTcc                                                         >  1:1607135/1‑93 (MQ=255)
                                     cAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCt                                                        >  2:5365684/1‑93 (MQ=255)
                                      aaCTGGCATTCACTTTCTCGGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCt                                                        >  1:6158151/1‑93 (MQ=255)
                                      aaCTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCtt                                                       >  2:1294168/1‑93 (MQ=255)
                                       aCTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTc                                                      >  2:2897060/1‑93 (MQ=255)
                                       aCTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTc                                                      <  2:3305172/93‑1 (MQ=255)
                                         tGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGc                                                    <  2:5805224/93‑1 (MQ=255)
                                         tGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGc                                                    >  2:3682877/1‑93 (MQ=255)
                                         tGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGc                                                    <  1:3570321/93‑1 (MQ=255)
                                            cATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCAc                                                 >  1:5069385/1‑93 (MQ=255)
                                             aTTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCAcc                                                >  1:4755135/1‑93 (MQ=255)
                                              ttCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCg                                               >  2:6788231/1‑93 (MQ=255)
                                              ttCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCg                                               >  1:3313392/1‑93 (MQ=255)
                                              ttCAATTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCg                                               <  1:5728163/93‑1 (MQ=255)
                                               tCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGt                                              <  2:1757005/93‑1 (MQ=255)
                                                 aCTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCg                                            <  2:6869034/93‑1 (MQ=255)
                                                 aCTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCg                                            >  2:5147649/1‑93 (MQ=255)
                                                   tttCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGc                                          >  2:2330448/1‑93 (MQ=255)
                                                     tctc‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCaa                                        >  2:3639923/1‑93 (MQ=255)
                                                     tctc‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCaa                                        >  2:5447391/1‑93 (MQ=255)
                                                        c‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCg                                     <  1:2662927/93‑1 (MQ=255)
                                                        c‑GCCCAGCAGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCg                                     <  2:2162216/93‑1 (MQ=255)
                                                          gCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCgg                                    >  2:4366111/1‑93 (MQ=255)
                                                          gCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCgg                                    >  1:3728735/1‑93 (MQ=255)
                                                          gCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCgg                                    <  1:4510377/93‑1 (MQ=255)
                                                            ccAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                  <  1:3498519/93‑1 (MQ=255)
                                                             cAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAAc                                 <  1:4282836/93‑1 (MQ=255)
                                                             cAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAAc                                 <  2:448442/93‑1 (MQ=255)
                                                             cAGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAAc                                 <  1:1580561/93‑1 (MQ=255)
                                                              aGCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACa                                <  2:914211/93‑1 (MQ=255)
                                                               gCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAg                               >  1:6869149/1‑93 (MQ=255)
                                                               gCGGCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAg                               >  2:6125155/1‑93 (MQ=255)
                                                                  gCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGaaa                            >  1:6503147/1‑93 (MQ=255)
                                                                  gCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGaaa                            <  2:5693557/93‑1 (MQ=255)
                                                                  gCCCTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGaaa                            <  1:7129771/93‑1 (MQ=255)
                                                                  gCCCTGAAACCCAGGCGTGACGTTCTCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGaaa                            <  1:6827128/93‑1 (MQ=255)
                                                                   cccTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAAc                           >  2:2547301/1‑93 (MQ=255)
                                                                     cTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGt                         <  1:209543/93‑1 (MQ=255)
                                                                     cTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGt                         >  1:4218836/1‑93 (MQ=255)
                                                                     cTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGt                         <  2:4156424/93‑1 (MQ=255)
                                                                     cTGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGt                         >  1:2019370/1‑93 (MQ=255)
                                                                      tGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                        <  1:4098851/93‑1 (MQ=255)
                                                                      tGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                        <  2:4852513/93‑1 (MQ=255)
                                                                      tGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                        >  2:5091876/1‑93 (MQ=255)
                                                                      tGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                        <  2:5563172/93‑1 (MQ=255)
                                                                      tGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                        >  1:2831799/1‑93 (MQ=255)
                                                                      tGTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                        >  2:7316105/1‑93 (MQ=255)
                                                                       gTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTAc                       <  2:2634339/93‑1 (MQ=255)
                                                                       gTAACCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTAc                       <  2:6543799/93‑1 (MQ=255)
                                                                         aaCCCAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACtg                     >  2:3118180/1‑93 (MQ=255)
                                                                           cccAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACtgtg                   >  2:6285338/1‑93 (MQ=255)
                                                                             cAGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCa                 <  1:6884103/93‑1 (MQ=255)
                                                                             cAGGCGTGACGTTCACCCCCCCGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGCCGGCAACCGGAACAGAAACGTACTGTGCa                 >  1:1761145/1‑93 (MQ=255)
                                                                              aGGCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCAt                <  2:1856788/93‑1 (MQ=255)
                                                                                gCGTGACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCa              <  1:806490/93‑1 (MQ=255)
                                                                                 cGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAg             <  1:2979793/93‑1 (MQ=255)
                                                                                    gACGTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaa          >  2:7401841/1‑93 (MQ=255)
                                                                                     aCGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaaa         <  1:342017/93‑1 (MQ=255)
                                                                                      cGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaaaa        <  2:1060110/93‑1 (MQ=255)
                                                                                       gTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaaaaa       <  2:601426/93‑1 (MQ=255)
                                                                                       gTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaaaaa       <  1:953116/93‑1 (MQ=255)
                                                                                       gTTCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGaaaaa       >  1:1123740/1‑93 (MQ=255)
                                                                                        ttCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAAt      >  2:4610062/1‑93 (MQ=255)
                                                                                         tCACCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAAtt     <  1:5446732/93‑1 (MQ=255)
                                                                                            gCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTcc  >  2:5852664/1‑93 (MQ=255)
                                                                                            |                                                                                            
TTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTC‑GCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCC  >  NC_000913/3269145‑3269328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑