Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,500,673 C→T T109T (ACG→ACA yi41 ← IS4 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,500,6730CT85.9% 822.7 / 94.5 283T109T (ACG→ACAyi41IS4 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base T (112/131);  total (130/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 297 total aligned reads displayed.

GTGTGCGCCAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCAT  >  NC_000913/4500582‑4500765
                                                                                           |                                                                                            
  gtgCGCCAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTtgtg                                                                                           >  1:281355/1‑93 (MQ=255)
   tgCGCCAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGa                                                                                          <  2:6466842/93‑1 (MQ=255)
      gcCAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAaca                                                                                       >  1:805743/1‑93 (MQ=255)
       ccAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAacac                                                                                      <  1:3606620/93‑1 (MQ=255)
       ccAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAacac                                                                                      >  2:4902439/1‑93 (MQ=255)
       ccAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAacac                                                                                      <  2:5672693/93‑1 (MQ=255)
        cAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACAcg                                                                                     >  2:5848705/1‑93 (MQ=255)
         aGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACAcgc                                                                                    >  1:1836178/1‑93 (MQ=255)
          gAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACAcgcg                                                                                   <  2:6254109/93‑1 (MQ=255)
           aaCACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACgcgg                                                                                  <  1:6497583/93‑1 (MQ=255)
            acacacCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACgcggc                                                                                 <  2:4097616/93‑1 (MQ=255)
              acacCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGcggcgg                                                                               >  1:3700303/1‑93 (MQ=255)
              acacCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGcggcgg                                                                               <  2:6505816/93‑1 (MQ=255)
               cacCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGa                                                                              >  2:1623975/1‑93 (MQ=255)
                acCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGAc                                                                             <  2:5961646/93‑1 (MQ=255)
                 ccATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACa                                                                            >  1:5802558/1‑93 (MQ=255)
                 ccATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACa                                                                            <  1:5741628/93‑1 (MQ=255)
                 ccATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACa                                                                            >  2:4265451/1‑93 (MQ=255)
                 ccATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACa                                                                            <  2:4770410/93‑1 (MQ=255)
                 ccATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACa                                                                            >  1:1723375/1‑93 (MQ=255)
                 ccATCGATGGCCAGCAGATTCAGGCCGTCCCAGTGCGGATGCGGCGTGGCGTTATGCCAGCGCTGCGCTGTTTTTGTGAACACGcggcgggcc                                                                            >  1:1803757/1‑90 (MQ=255)
                    tCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGcc                                                                         <  1:4814100/93‑1 (MQ=255)
                     cGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCt                                                                        <  2:2922253/93‑1 (MQ=255)
                     cGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCt                                                                        <  1:1023949/93‑1 (MQ=255)
                        tGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCAc                                                                     <  1:2875818/93‑1 (MQ=255)
                        tGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGGTATGCCAGAGCTGCGCCGTTTTTGTGTACACGCGGCGGACAGCCTCAc                                                                     >  2:4579545/1‑93 (MQ=255)
                         ggCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACt                                                                    >  1:2648450/1‑93 (MQ=255)
                          gCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACtt                                                                   >  1:891685/1‑93 (MQ=255)
                          gCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACtt                                                                   >  1:2540066/1‑93 (MQ=255)
                             agcagGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTccc                                                                >  1:1897516/1‑93 (MQ=255)
                             agcagGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTccc                                                                <  1:2100193/93‑1 (MQ=255)
                             agcagGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTccc                                                                <  1:5201801/93‑1 (MQ=255)
                              gcagGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCa                                                               <  1:6623155/93‑1 (MQ=255)
                               cagGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAg                                                              >  2:5241915/1‑93 (MQ=255)
                                agGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAgg                                                             >  2:5524290/1‑93 (MQ=255)
                                  ggTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGgcg                                                           >  2:3994811/1‑93 (MQ=255)
                                   gTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGgcgc                                                          >  2:1446203/1‑93 (MQ=255)
                                   gTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGgcgc                                                          <  2:6375411/93‑1 (MQ=255)
                                   gTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGgcgc                                                          >  1:5343322/1‑93 (MQ=255)
                                    tCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCt                                                         <  1:4156508/93‑1 (MQ=255)
                                    tCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCt                                                         >  1:1472266/1‑93 (MQ=255)
                                    gcAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCt                                                         <  2:5578115/92‑1 (MQ=255)
                                     cAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTg                                                        >  1:5451593/1‑93 (MQ=255)
                                     cAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTg                                                        >  1:1335124/1‑93 (MQ=255)
                                       ggCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGc                                                      >  1:4095174/1‑93 (MQ=255)
                                       ggCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGc                                                      <  2:1240877/93‑1 (MQ=255)
                                       ggCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGc                                                      >  2:2878419/1‑93 (MQ=255)
                                          cGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCggg                                                   <  1:3363852/93‑1 (MQ=255)
                                           gCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGTGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGc                                                  >  2:6058195/1‑93 (MQ=255)
                                           gCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGc                                                  <  2:3490357/93‑1 (MQ=255)
                                            cACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGcc                                                 <  2:2747607/93‑1 (MQ=255)
                                            cACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGcc                                                 >  1:1070336/1‑93 (MQ=255)
                                              ccAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTg                                               >  1:5672355/1‑93 (MQ=255)
                                                aGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGaa                                             >  1:727963/1‑93 (MQ=255)
                                                  tGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATa                                           >  1:1132871/1‑93 (MQ=255)
                                                    cGGATGGGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATaag                                         >  2:2638791/1‑92 (MQ=255)
                                                    cGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAAc                                         <  1:3180426/93‑1 (MQ=255)
                                                    cGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAAc                                         <  2:4694330/93‑1 (MQ=255)
                                                    cGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAAc                                         >  1:2304591/1‑93 (MQ=255)
                                                      gATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACgg                                       >  2:1192127/1‑93 (MQ=255)
                                                      gATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACgg                                       <  1:4979501/93‑1 (MQ=255)
                                                       aTGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGCGCCTGAATAACgga                                      >  2:5911858/1‑92 (MQ=255)
                                                         gcggcgTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCAc                                    >  1:6346093/1‑93 (MQ=255)
                                                         gcggcgTGGCGTTATGCCAGAGCTGCGCTGGTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCAc                                    <  2:2588858/93‑1 (MQ=255)
                                                          cggcgTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACt                                   <  2:3269746/93‑1 (MQ=255)
                                                          cggcgTGGCGTTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCAGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACt                                   >  1:2993035/1‑93 (MQ=255)
                                                            gcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgg                                 >  2:3678595/1‑93 (MQ=255)
                                                            gcgtggcgtTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgg                                 >  2:706117/1‑93 (MQ=255)
                                                             cgtggcgtTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggg                                <  2:2348297/93‑1 (MQ=255)
                                                              gtggcgtTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                               <  1:3920593/93‑1 (MQ=255)
                                                              gtggcgtTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTgggg                               >  2:4132526/1‑93 (MQ=255)
                                                               tggcgtTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTggggg                              >  2:831106/1‑93 (MQ=255)
                                                                ggcgtTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGc                             >  2:3840359/1‑93 (MQ=255)
                                                                   gtTATGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAAc                          <  2:3178764/93‑1 (MQ=255)
                                                                   gtTATGCCAGAGCTGCGCTGTTTTCGTGACCACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAAc                          <  2:4182833/93‑1 (MQ=255)
                                                                      aTGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGaa                       <  1:2426197/93‑1 (MQ=255)
                                                                      aTGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGaa                       <  1:3760996/93‑1 (MQ=255)
                                                                      aTGCCAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGaa                       <  2:2152187/93‑1 (MQ=255)
                                                                         ccAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAggg                    >  1:4700105/1‑93 (MQ=255)
                                                                          cAGAGCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAgggg                   <  2:4412269/93‑1 (MQ=255)
                                                                           agagCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGc                  >  1:6086409/1‑93 (MQ=255)
                                                                           agagCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGc                  <  2:5488888/93‑1 (MQ=255)
                                                                             agCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGa                <  2:5273304/93‑1 (MQ=255)
                                                                             agCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGa                <  2:757847/93‑1 (MQ=255)
                                                                              gCTGCGCTGTTTTTGTGAACGCGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGAt               <  1:3284879/93‑1 (MQ=255)
                                                                              gCTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGAt               >  1:2042506/1‑93 (MQ=255)
                                                                               cTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGAtt              >  2:381691/1‑93 (MQ=255)
                                                                               cTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGAtt              <  1:5928570/93‑1 (MQ=255)
                                                                               cTGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGAtt              <  2:1509332/93‑1 (MQ=255)
                                                                                tGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGGTTg             <  2:1963084/93‑1 (MQ=255)
                                                                                tGCGCTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTg             <  2:785762/93‑1 (MQ=255)
                                                                                 gcgcTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGc            <  2:5028363/93‑1 (MQ=255)
                                                                                 gcgcTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGc            >  2:2253572/1‑93 (MQ=255)
                                                                                   gcTGTTTTTGTGAACACGCGGCGGACACCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCGCTGGGGGCAACGAAGGGGCGATTGccc          <  2:5097101/93‑1 (MQ=255)
                                                                                    cTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCg         >  2:3873620/1‑93 (MQ=255)
                                                                                    cTGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCg         >  2:1020114/1‑93 (MQ=255)
                                                                                     tGTTTTTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCgg        <  2:862594/93‑1 (MQ=255)
                                                                                       ttttTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGgca      <  1:4380415/93‑1 (MQ=255)
                                                                                       ttttTGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGgca      >  1:6141962/1‑93 (MQ=255)
                                                                                          tCGTGAACACGCGGCGGACAGCCTCCCTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGgcagca   <  1:1302594/93‑1 (MQ=255)
                                                                                           |                                                                                            
GTGTGCGCCAGAACACACCATCGATGGCCAGCAGGGTCAGGCCGCACCAGTGCGGATGCGGCGTGGCGTTATGCCAGAGCTGCGCTGTTTTCGTGAACACGCGGCGGACAGCCTCACTTCCCAGGCGCTGGCGGGCCTGAATAACGGCACTGGGGGCAACGAAGGGGCGATTGCCCGGCAGCAT  >  NC_000913/4500582‑4500765

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