Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,900 G→A G113G (GGC→GGT yieO ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,9000GA82.1% 735.2 / 118.5 302G113G (GGC→GGTyieOputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base A (142/105);  total (172/130)
Fisher's exact test for biased strand distribution p-value = 8.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 313 total aligned reads displayed.

CCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACT  >  NC_000913/3937809‑3937989
                                                                                           |                                                                                         
   gCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcgccgc                                                                                       >  1:4576308/1‑93 (MQ=255)
   gCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcgccgc                                                                                       >  1:98689/1‑93 (MQ=255)
   gCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcaccgc                                                                                       >  2:7397976/1‑93 (MQ=255)
      tGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTa                                                                                    >  2:2997997/1‑93 (MQ=255)
      tGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTa                                                                                    >  2:5464068/1‑93 (MQ=255)
      tGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTa                                                                                    <  2:5074559/93‑1 (MQ=255)
        gCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATc                                                                                  >  2:7181965/1‑93 (MQ=255)
        gCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATc                                                                                  <  2:6398140/93‑1 (MQ=255)
          gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATccc                                                                                >  1:5108629/1‑93 (MQ=255)
          gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATccc                                                                                >  2:5610407/1‑93 (MQ=255)
           aCAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATcccc                                                                               <  1:3239832/93‑1 (MQ=255)
           aCAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATcccc                                                                               >  1:653388/1‑93 (MQ=255)
           aCAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATcccc                                                                               >  2:5316229/1‑93 (MQ=255)
            cAAAATTCAATCCTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCt                                                                              <  1:6013641/93‑1 (MQ=255)
            cAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCt                                                                              >  2:5821140/1‑93 (MQ=255)
            cAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCt                                                                              >  2:7516142/1‑93 (MQ=255)
            cAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCt                                                                              <  1:5697073/93‑1 (MQ=255)
               aaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGaa                                                                           >  1:7291763/1‑93 (MQ=255)
               aaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGaa                                                                           >  1:4293669/1‑93 (MQ=255)
                 ttCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATa                                                                         >  1:5403300/1‑93 (MQ=255)
                 ttCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATa                                                                         >  2:4130968/1‑93 (MQ=255)
                        cTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGa                                                                  <  2:1816272/93‑1 (MQ=255)
                         tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGaa                                                                 >  1:3993172/1‑93 (MQ=255)
                          ggAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAg                                                                >  1:6297723/1‑93 (MQ=255)
                           gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGa                                                               <  1:2568134/93‑1 (MQ=255)
                           gaagaaGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGa                                                               >  2:1531643/1‑93 (MQ=255)
                            aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGAc                                                              >  1:2489713/1‑93 (MQ=255)
                            aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGAc                                                              <  2:6901542/93‑1 (MQ=255)
                            aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGAc                                                              >  1:3329257/1‑93 (MQ=255)
                             agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACa                                                             >  1:6079922/1‑93 (MQ=255)
                             agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACa                                                             >  2:2902674/1‑93 (MQ=255)
                             agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACa                                                             >  1:839317/1‑93 (MQ=255)
                             agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACa                                                             <  2:4221780/93‑1 (MQ=255)
                             agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACa                                                             <  1:4827517/93‑1 (MQ=255)
                               aagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAAc                                                           <  1:210468/93‑1 (MQ=255)
                                agTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAAcc                                                          >  2:4270507/1‑93 (MQ=255)
                                 gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGTAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCa                                                         >  1:6216147/1‑93 (MQ=255)
                                    cATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCt                                                      >  2:5755002/1‑93 (MQ=255)
                                    cATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCt                                                      <  2:2461456/93‑1 (MQ=255)
                                     aTTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtg                                                     >  1:6717878/1‑93 (MQ=255)
                                     aTTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtg                                                     >  2:1374144/1‑93 (MQ=255)
                                      ttACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtgt                                                    <  2:1481933/93‑1 (MQ=255)
                                         cGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGt                                                 >  2:3666081/1‑93 (MQ=255)
                                           aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAg                                               >  2:4869073/1‑93 (MQ=255)
                                           aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAg                                               >  2:4687164/1‑93 (MQ=255)
                                            ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              >  2:1200153/1‑93 (MQ=255)
                                            ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              <  1:3268125/93‑1 (MQ=255)
                                            ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              <  2:4363681/93‑1 (MQ=255)
                                            ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              >  2:4754730/1‑93 (MQ=255)
                                             gATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCg                                             >  1:1012251/1‑93 (MQ=255)
                                              aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGa                                            >  2:6820261/1‑93 (MQ=255)
                                               tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGaa                                           >  2:1980558/1‑93 (MQ=255)
                                               tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGaa                                           <  2:3505364/93‑1 (MQ=255)
                                               tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGaa                                           <  2:2629651/93‑1 (MQ=255)
                                                aaGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAt                                          <  2:5176156/93‑1 (MQ=255)
                                                aaGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAt                                          <  2:6705055/93‑1 (MQ=255)
                                                  gcgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTa                                        <  2:5277996/93‑1 (MQ=255)
                                                  gcgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTa                                        <  1:6990797/93‑1 (MQ=255)
                                                  gcgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTa                                        <  1:4915034/93‑1 (MQ=255)
                                                   cgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAg                                       <  1:2800835/93‑1 (MQ=255)
                                                    gcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGa                                      >  1:2670706/1‑93 (MQ=255)
                                                    gcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGa                                      >  1:3406479/1‑93 (MQ=255)
                                                    gcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGa                                      <  2:3299954/93‑1 (MQ=255)
                                                      gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGaaa                                    <  1:1425134/93‑1 (MQ=255)
                                                      gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGaaa                                    <  2:1960880/93‑1 (MQ=255)
                                                       cAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAg                                   <  2:5905100/93‑1 (MQ=255)
                                                        aGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                  <  1:162624/93‑1 (MQ=255)
                                                         gTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTg                                 <  1:6933301/93‑1 (MQ=255)
                                                          tAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgc                                >  2:605624/1‑93 (MQ=255)
                                                            aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              <  1:4767753/93‑1 (MQ=255)
                                                            aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              >  1:290526/1‑93 (MQ=255)
                                                             ggCCAGCCGAGCAACGGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCa                             >  1:5550550/1‑93 (MQ=255)
                                                              nncAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAg                            <  1:1743701/91‑1 (MQ=255)
                                                              gCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAg                            <  1:2264926/93‑1 (MQ=255)
                                                                cAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAgnn                          >  1:1975110/1‑91 (MQ=255)
                                                                cAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGc                          >  1:3549888/1‑93 (MQ=255)
                                                                    cGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                      <  2:2212087/93‑1 (MQ=255)
                                                                    cGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                      >  2:7338476/1‑93 (MQ=255)
                                                                     gAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCagag                     >  1:4406841/1‑93 (MQ=255)
                                                                       gCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGaa                   <  1:3034346/93‑1 (MQ=255)
                                                                       gCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGaa                   >  1:5873434/1‑93 (MQ=255)
                                                                         aaCAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAcc                 >  2:2042363/1‑93 (MQ=255)
                                                                          aCAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                >  1:4729728/1‑93 (MQ=255)
                                                                             ggCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  1:6404908/1‑93 (MQ=255)
                                                                             ggCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  2:6218216/1‑93 (MQ=255)
                                                                             ggCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  1:7405761/1‑93 (MQ=255)
                                                                             ggCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  1:1540049/1‑93 (MQ=255)
                                                                             ggCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  2:4573533/1‑93 (MQ=255)
                                                                              gCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtg            <  2:3162516/93‑1 (MQ=255)
                                                                              gCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtg            <  2:4024692/93‑1 (MQ=255)
                                                                               catcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgt           <  2:3855305/93‑1 (MQ=255)
                                                                                 tcataatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGa         <  1:3653830/93‑1 (MQ=255)
                                                                                  catcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGaa        >  2:5229020/1‑93 (MQ=255)
                                                                                   atcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAAc       >  2:2400089/1‑93 (MQ=255)
                                                                                    tcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACa      >  2:1110727/1‑93 (MQ=255)
                                                                                      atTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGa    <  1:735703/93‑1 (MQ=255)
                                                                                      atTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGa    >  2:7000212/1‑93 (MQ=255)
                                                                                       ttGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGAc   >  1:5783594/1‑93 (MQ=255)
                                                                                       ttGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGAc   >  2:6587249/1‑93 (MQ=255)
                                                                                        tGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACt  >  1:4181281/1‑93 (MQ=255)
                                                                                           |                                                                                         
CCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACT  >  NC_000913/3937809‑3937989

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