Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,958,261 Δ1 bp intergenic (‑830/‑31) ppiC ← / → rep peptidyl‑prolyl cis‑trans isomerase C/rep helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,958,2610T.75.8% 268.6 / 116.2 127intergenic (‑830/‑31)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  major base . (27/67);  minor base C (12/18);  total (41/86)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 226 total aligned reads displayed.

GGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCTCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGGTGC  >  NC_000913/3958174‑3958352
                                                                                          |                                                                                           
tgtcttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccc                                                                                           <  1:4418981‑M2/5‑1 (MQ=255)
 gtcttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatccccccccc                                                                                          >  2:4017505‑M2/88‑93 (MQ=255)
 gtcttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatccccccccc                                                                                          >  1:3860524‑M2/88‑93 (MQ=255)
 gtcttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatccccccccc                                                                                          <  1:5560728‑M2/6‑1 (MQ=255)
 gtcttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatccccccccc                                                                                          >  1:909017‑M2/88‑93 (MQ=255)
  tcttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccc                                                                                         >  2:6797187‑M2/87‑93 (MQ=255)
  tcttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccc                                                                                         >  1:2464082‑M2/87‑93 (MQ=255)
  gcttctctctctgatttaaacatccgcagccagccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccc                                                                                         <  2:7000622‑M2/7‑1 (MQ=255)
   cttctctctctgatttaaacccccccagccaaccggttaggggcttccccgcggtcacattcaaatgcgattctgctacaacccccccccccc                                                                                        <  2:5543939‑M2/8‑1 (MQ=255)
    ttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccccg                                                                                       <  1:7702664‑M2/9‑2 (MQ=255)
    ttctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccccc                                                                                       >  1:1672736‑M2/85‑93 (MQ=255)
    ttctctctctgatttaaaaatccgcagccaaccggttagcggcttacacgcggtcacattcaaattcgattctgctacaatctccccccccat                                                                                       >  1:4260099‑M2/85‑91 (MQ=255)
    ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TACCCCCCCCGtt                                                                                    >  1:3655201/1‑93 (MQ=255)
     tctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccccgc                                                                                      >  2:3217542‑M2/84‑91 (MQ=255)
      ctctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccccccg                                                                                     >  1:2294093‑M2/83‑91 (MQ=255)
       tctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccccccccccccgt                                                                                    >  2:5336010‑M2/82‑90 (MQ=255)
       tctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCCCtt                                                                                    >  2:6931421‑M2/82‑93 (MQ=255)
       tctctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCCCCCCCGTTCGa                                                                                 >  1:4328790/1‑93 (MQ=255)
        ctctctgatttaaacatccgccgccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGa                                                                                 <  1:1564872‑M2/13‑1 (MQ=255)
        ctctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCg                                                                                  <  1:6315338‑M2/13‑1 (MQ=255)
        ccttctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcccattcaaatgcgattctgcttcaaacccCCCCCCCCGTTCg                                                                                  <  2:89882‑M2/13‑1 (MQ=255)
         tctctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCCCTtcc                                                                                  >  2:5758004‑M2/80‑92 (MQ=255)
           tctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCCTTCGAAg                                                                               >  1:3486558‑M2/78‑93 (MQ=255)
            ctgatttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGAt                                                                             <  2:6594945‑M2/17‑1 (MQ=255)
             tGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAC‑‑TCCCCCCCCCGTTCGAAGATTg                                                                           <  1:1779126/93‑1 (MQ=255)
              gATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCCCCCCCGTTCGAAGATTGa                                                                          >  1:6615113/1‑93 (MQ=255)
                tttaaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCCTCCGCAGATtgg                                                                          >  1:9418767‑M2/73‑92 (MQ=255)
                   aaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCa                                                                       <  1:8624318‑M2/24‑1 (MQ=255)
                   aaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCCTTCCACgattggaca                                                                       >  1:4561550‑M2/70‑89 (MQ=255)
                   aaacatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCaa                                                                      <  2:742600‑M2/24‑1 (MQ=255)
                    aacatccgccgccaaccggttagcggcttacacgcggtcacattcaaatgcggttctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCaa                                                                      <  1:1065023‑M2/25‑1 (MQ=255)
                     acatccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAAt                                                                     <  1:1436907‑M2/26‑1 (MQ=255)
                     aCATCCGCAGCCAACCGGTTAACGGCCTACCCGCGCTCACATTCAAATGCGATTCTGCCAAAA‑‑TCCCCCCCCCCGTTCGAAGATTGAGCAATa                                                                    <  2:102749/93‑1 (MQ=255)
                      catccgcagccaacgggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGGTCGCAGATTTggcaaaaa                                                                   >  2:7501632‑M2/67‑90 (MQ=255)
                      aatccccagcccacccgttagcggcttacacggcgtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATa                                                                    <  2:528114‑M2/27‑1 (MQ=255)
                       atccgcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTCCGAGGATTGAGCAATac                                                                   >  2:1651089‑M2/66‑93 (MQ=255)
                         ccgcagccaaccggttagcggcttacccgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATacac                                                                 <  1:5683583‑M2/30‑1 (MQ=255)
                         ccccAGCCAACCCGTTTGCCGCCTTCCCGGCGTCACAATTAAATTCGATTCTGCTACAAACCCCCCCCCCCGTTCGAAGATTGAGCAATacac                                                                 <  2:7342103/90‑1 (MQ=255)
                           gcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcgcCCCCCCCCGTTCCATGATTACGCACTAAacca                                                               >  1:9863444‑M2/62‑92 (MQ=255)
                           gcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTa                                                              <  1:8107813‑M2/32‑1 (MQ=255)
                           gcagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctactacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTa                                                              <  1:4708545‑M2/32‑1 (MQ=255)
                           ccAGCCAACCGGTTAGGGGGTTCCCCGGGGCCCAATTAAAATGGGTTTTGGTAAAAACCCCCCCCCCCGTTCGAAGATTGAGCAATACACCTa                                                              <  2:2048274/92‑1 (MQ=255)
                            ccagccaaccggttaggggcttacacccggtcacattcaaatgcgattctgctaaaacccCCCCCCCCCGTTCGAAGATTGAGCAATACACCt                                                               <  2:1315960‑M2/33‑1 (MQ=255)
                            cagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTa                                                              <  2:3664626‑M2/33‑1 (MQ=255)
                            cagccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTAt                                                             <  2:483309‑M2/33‑1 (MQ=255)
                             gccaaccggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatccccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTa                                                              <  2:1578511‑M2/34‑1 (MQ=255)
                              gcgcaaccggttcgcggcttacacgggttcacattcaaatgcgtttctgctacaacccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTAt                                                             <  2:5658480‑M2/35‑1 (MQ=255)
                              gccaaccgggtagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATg                                                            <  1:6476733‑M2/35‑1 (MQ=255)
                                  accggttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACCCGTATGAGTCt                                                       >  2:2997471‑M2/55‑93 (MQ=255)
                                      gttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCACCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTaa                                                     <  1:9590341‑M2/43‑1 (MQ=255)
                                      gttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGACGAGTGGGGAATACACCTATGGGTCTaaa                                                    >  2:2753458‑M2/51‑93 (MQ=255)
                                      gttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTaaa                                                    <  1:6284342‑M2/43‑1 (MQ=255)
                                       ttagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAAcc                                                  >  2:829100‑M2/50‑93 (MQ=255)
                                       gtagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAAc                                                   <  2:1150976‑M2/44‑1 (MQ=255)
                                        tagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGACGATTGAGCAATACACCTATGCGCCTAAAccc                                                 >  2:6113260‑M2/49‑93 (MQ=255)
                                        tagcggcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAAccc                                                 <  2:8706114‑M2/45‑1 (MQ=255)
                                        tAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCCCCCCCCCTTTTAAGGATTACCAATAAACCTATGCGGTTAAcccc                                                 >  2:6427257/1‑93 (MQ=255)
                                        tAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAAcccc                                                <  1:5443746/93‑1 (MQ=255)
                                         agcggcttacccggggtcacattcaaatgccattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAAccc                                                 <  1:4180084‑M2/46‑1 (MQ=255)
                                         agcggcttacacgcggtcacatttaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAAccc                                                 <  1:4673532‑M2/46‑1 (MQ=255)
                                             gcttacacggggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGACGATTGAGCAATACACCTATGCGTCTAAACCCCGGcc                                            >  2:6732199‑M2/44‑93 (MQ=255)
                                             gcttacacgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGcc                                            <  1:3010520‑M2/50‑1 (MQ=255)
                                                tacacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaa                                          <  2:5000836‑M2/53‑1 (MQ=255)
                                                 tcacgcggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaac                                         <  1:3372139‑M2/54‑1 (MQ=255)
                                                  cacgcggttacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaacaa                                       <  1:3272606‑M2/55‑1 (MQ=255)
                                                  cacgcggtccccttcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaaca                                        <  1:7811873‑M2/55‑1 (MQ=255)
                                                   acgcggtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaacaac                                      <  1:7434421‑M2/56‑1 (MQ=255)
                                                    ccgcggtcacattcaaatgcgattctgctacaacccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaacaa                                       <  1:3211655‑M2/57‑1 (MQ=255)
                                                      cggtcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCcaacaacaa                                    <  2:4675230‑M2/59‑1 (MQ=255)
                                                      cgGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCCGCCAACAACAAGc                                  >  1:9912845/1‑93 (MQ=255)
                                                        gtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCt                                 <  1:7026253‑M2/61‑1 (MQ=255)
                                                        gtcacattcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGAAATACACCTCTGCTTCCAAACCACGGCCAACAACAAGCt                                 >  2:637214‑M2/33‑93 (MQ=255)
                                                        gTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTg                                <  1:9807485/93‑1 (MQ=255)
                                                         tcacattcaaatgcgattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCt                                 <  1:1282625‑M2/62‑1 (MQ=255)
                                                          cacaataaaaagcggttcttccacaatcccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCt                                 <  2:8037568‑M2/63‑1 (MQ=255)
                                                            cattcaaatggcattcttccacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTc                              <  2:1770762‑M2/65‑1 (MQ=255)
                                                            cattcaaatgctattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTc                              <  1:4087391‑M2/65‑1 (MQ=255)
                                                              ttcaaatgcgattctgctacaatcccC‑‑CCCCCCGTTCGCAGATTGAGCAACACACCTATGCGTCTAAACCCCGGCCCACAACAAACTGTCGaa                           >  1:6225552‑M2/27‑93 (MQ=255)
                                                                 caaatgcgattctgctacaacccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAAtt                         <  2:7157703‑M2/70‑1 (MQ=255)
                                                                 aaatgcgattctgctacaaacccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAAtt                         <  1:6366238‑M2/70‑1 (MQ=255)
                                                                   atgggattctgctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCg                       <  2:3926401‑M2/72‑1 (MQ=255)
                                                                   atgcgattctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGt                      <  2:7567549‑M2/72‑1 (MQ=255)
                                                                    agatattcttctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAGCCCCGGCCAACAACAAGCTGTCGAATTCGtt                     <  2:6285136‑M2/73‑1 (MQ=255)
                                                                      ccgttctgccaaaaacccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTa                    <  2:4433663‑M2/75‑1 (MQ=255)
                                                                       gattctgttacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAc                   <  1:6667132‑M2/76‑1 (MQ=255)
                                                                        attctgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACCCCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCg                 <  1:8754726‑M2/77‑1 (MQ=255)
                                                                         ttctgctacaaacccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGTTGTCGAATTCGTTACCg                 <  2:3102697‑M2/78‑1 (MQ=255)
                                                                         atgtgcaacaaccccCCCCCCCCAGTTCAAAGATTGAGCAATACACCTATGCTTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAcc                  <  2:5594756‑M2/78‑1 (MQ=255)
                                                                           tctgctacaacccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGc               <  1:3424805‑M2/80‑1 (MQ=255)
                                                                            tgctacaatcccC‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAACCTGCCCAATTCGTTACCGGccc             >  1:3310647‑M2/13‑93 (MQ=255)
                                                                             cctccaaacccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGcc              <  2:4354739‑M2/82‑1 (MQ=255)
                                                                              ctacaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGcccc            <  1:3972083‑M2/83‑1 (MQ=255)
                                                                              ctacaaacccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGcccc            <  2:2125322‑M2/83‑1 (MQ=255)
                                                                                ccaaccccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCcct           <  1:9286680‑M2/85‑1 (MQ=255)
                                                                                acaatcccCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCcctg          <  1:8748104‑M2/85‑1 (MQ=255)
                                                                                  ccaaccCCCCCCCCCGTTCGAAGATTGAGCAATACCCCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCcctgc         <  2:7981559‑M2/87‑1 (MQ=255)
                                                                                    atccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCcctgcct       <  2:4713267‑M2/89‑1 (MQ=255)
                                                                                    aatcCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCcctgcct       <  1:2857189‑M2/89‑1 (MQ=255)
                                                                                     cccCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCcctgcctg      <  2:179460‑M2/90‑1 (MQ=255)
                                                                                     accCCCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGTTGTCGAATTCGTTACCGGCCcctgcctg      <  1:5166070‑M2/90‑1 (MQ=255)
                                                                                          |                                                                                           
GGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCTCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGGTGC  >  NC_000913/3958174‑3958352

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