Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,490,510 C→T S6L (TCA→TTA)  trg → methyl‑accepting chemotaxis protein III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,490,5100CT99.2% 840.0 / ‑5.1 237S6L (TCA→TTA) trgmethyl‑accepting chemotaxis protein III
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (105/128);  total (108/129)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 241 total aligned reads displayed.

TTGGTCATAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTCACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCATTAAGT  >  NC_000913/1490421‑1490598
                                                                                         |                                                                                        
   gTCATAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCg                                                                                    <  1:4099987/93‑1 (MQ=255)
   gTCATAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCg                                                                                    >  2:7134604/1‑93 (MQ=255)
    tCATAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGa                                                                                   >  1:4330367/1‑93 (MQ=255)
     cATAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGAt                                                                                  <  2:1459697/93‑1 (MQ=255)
      aTAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGAtt                                                                                 >  2:6320228/1‑93 (MQ=255)
       tAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGGCGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTa                                                                                >  2:3420115/1‑93 (MQ=255)
       tAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTa                                                                                >  1:1016980/1‑93 (MQ=255)
       tAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTa                                                                                >  2:4245720/1‑93 (MQ=255)
       tAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTa                                                                                >  2:5591815/1‑93 (MQ=255)
       tAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTa                                                                                <  2:7435073/93‑1 (MQ=255)
        aaGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAg                                                                               <  1:6108931/93‑1 (MQ=255)
        aaGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAg                                                                               <  2:276363/93‑1 (MQ=255)
        aaGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAg                                                                               <  2:7843956/93‑1 (MQ=255)
         aGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAgg                                                                              >  1:6378155/1‑93 (MQ=255)
         aGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAgg                                                                              >  2:731700/1‑93 (MQ=255)
          gTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGt                                                                             <  2:4350731/93‑1 (MQ=255)
           tAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGtt                                                                            >  2:2558468/1‑93 (MQ=255)
            aaTTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGttt                                                                           <  1:654035/93‑1 (MQ=255)
             aTTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGtttt                                                                          <  1:7435245/93‑1 (MQ=255)
              ttACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGttttt                                                                         <  2:250560/93‑1 (MQ=255)
                 ccGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATATCACTCCCTTACAGCGATTAGGTTTTTTGc                                                                      >  1:4265679/1‑93 (MQ=255)
                 ccGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGc                                                                      <  1:6733558/93‑1 (MQ=255)
                  cGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGGCGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGca                                                                     <  1:4977190/93‑1 (MQ=255)
                  cGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGca                                                                     <  1:4821288/93‑1 (MQ=255)
                   gTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTCACAGCGATTAGGTTTTTTGcat                                                                    >  1:7755517/1‑93 (MQ=255)
                    tCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGcatc                                                                   >  1:5275939/1‑93 (MQ=255)
                      aaGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGcatcac                                                                 >  2:1265144/1‑93 (MQ=255)
                      aaGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGcatcac                                                                 <  2:1775780/93‑1 (MQ=255)
                          gCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGcatcacatca                                                             <  1:2757900/93‑1 (MQ=255)
                           ccGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAg                                                            >  1:4647857/1‑93 (MQ=255)
                           ccGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAg                                                            <  2:7868089/93‑1 (MQ=255)
                           ccGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAg                                                            <  2:54840/93‑1 (MQ=255)
                            cGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAgg                                                           <  1:3435887/93‑1 (MQ=255)
                            cGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAgg                                                           <  2:4157753/93‑1 (MQ=255)
                             gATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAggt                                                          >  1:7534010/1‑93 (MQ=255)
                              aTGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAggtt                                                         <  2:5309632/93‑1 (MQ=255)
                              aTGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAggtt                                                         >  1:6409763/1‑93 (MQ=255)
                                gACTTTCTATCAGGAGTAAACCTGGGCGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAggttgg                                                       <  2:1528755/93‑1 (MQ=255)
                                gACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAggttgg                                                       <  2:3664153/93‑1 (MQ=255)
                                gACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAggttgg                                                       >  1:2713796/1‑93 (MQ=255)
                                 aCTTTCTATCGGGAGGAAACCTGGGCGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAggttggt                                                      >  2:5094471/1‑93 (MQ=255)
                                        aTCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTa                                               <  2:6738281/93‑1 (MQ=255)
                                        aTCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTa                                               >  2:1158437/1‑93 (MQ=255)
                                        aTCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTa                                               <  2:2982275/93‑1 (MQ=255)
                                        aTCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTa                                               <  1:1029237/93‑1 (MQ=255)
                                         tCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTAt                                              <  1:4585644/93‑1 (MQ=255)
                                         tCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTAt                                              >  2:1195048/1‑93 (MQ=255)
                                              aGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGcc                                         <  1:1466830/93‑1 (MQ=255)
                                               gTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCt                                        >  1:316932/1‑93 (MQ=255)
                                               gTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCt                                        >  1:1548867/1‑93 (MQ=255)
                                                tAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTg                                       >  1:7910447/1‑93 (MQ=255)
                                                  aaCCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCa                                     <  1:7688509/93‑1 (MQ=255)
                                                   aCCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCAt                                    >  1:6753919/1‑93 (MQ=255)
                                                   aCCTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCAt                                    >  2:6539320/1‑93 (MQ=255)
                                                    ccTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCAtt                                   >  2:757843/1‑93 (MQ=255)
                                                     cTGGACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTc                                  <  2:3638459/93‑1 (MQ=255)
                                                      tGGACGAGAGACAACGGTAATGAATACCACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGACTGCATTcg                                 >  2:2898278/1‑92 (MQ=255)
                                                        gACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAg                               <  1:3685842/93‑1 (MQ=255)
                                                        gACGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAg                               <  2:2278519/93‑1 (MQ=255)
                                                         aCGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAgg                              <  2:1416086/93‑1 (MQ=255)
                                                         aCGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAgg                              <  2:6561705/93‑1 (MQ=255)
                                                          cGAGAGACAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGc                             >  1:4672834/1‑93 (MQ=255)
                                                           gagagaCAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCg                            >  1:3830076/1‑93 (MQ=255)
                                                           gagagaCAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCg                            >  2:1191362/1‑93 (MQ=255)
                                                            agagaCAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCgg                           <  2:5148472/93‑1 (MQ=255)
                                                             gagaCAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGt                          >  1:5488315/1‑93 (MQ=255)
                                                              agaCAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTa                         >  2:622648/1‑93 (MQ=255)
                                                              agaCAACGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTGTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTa                         <  1:4026618/93‑1 (MQ=255)
                                                                  aaCGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCtt                     <  1:1322724/93‑1 (MQ=255)
                                                                   aCGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTa                    <  1:5252514/93‑1 (MQ=255)
                                                                   aCGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTa                    >  2:2203967/1‑93 (MQ=255)
                                                                   aCGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTa                    >  1:5725070/1‑93 (MQ=255)
                                                                    cGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAg                   <  2:4611775/93‑1 (MQ=255)
                                                                    cGGTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAg                   >  2:2428157/1‑93 (MQ=255)
                                                                      gTAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGtt                 >  2:6155562/1‑93 (MQ=255)
                                                                       tAATGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTc                >  1:1839060/1‑93 (MQ=255)
                                                                        aaTGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCt               >  1:116463/1‑93 (MQ=255)
                                                                         aTGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTa              <  1:4703285/93‑1 (MQ=255)
                                                                         aTGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTa              <  2:516819/93‑1 (MQ=255)
                                                                         aTGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTa              <  1:2433264/93‑1 (MQ=255)
                                                                         aTGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTa              >  1:5157916/1‑93 (MQ=255)
                                                                          tGAATTCAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTAt             >  2:1873621/1‑93 (MQ=255)
                                                                          tGAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTAt             >  1:2791830/1‑93 (MQ=255)
                                                                           gAATACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTAtt            <  1:2209133/93‑1 (MQ=255)
                                                                            aaTACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTc           <  2:3468698/93‑1 (MQ=255)
                                                                            aaTACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTc           >  2:1307799/1‑93 (MQ=255)
                                                                            aaTACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTc           <  2:7075415/93‑1 (MQ=255)
                                                                              tACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGc         >  1:7803091/1‑93 (MQ=255)
                                                                              tACAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGAATCTTAGTTCTATTCGc         >  1:5073244/1‑93 (MQ=255)
                                                                               aCAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTCTCTTAGTTCTATTCGCa        <  1:796782/93‑1 (MQ=255)
                                                                               aCAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCa        <  2:2574314/93‑1 (MQ=255)
                                                                               aCAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCa        <  1:5576313/93‑1 (MQ=255)
                                                                               aCAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCa        <  1:3042894/93‑1 (MQ=255)
                                                                               aCAACTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCa        >  1:5927045/1‑93 (MQ=255)
                                                                               aCAACTCCCTTACAGCGATTAGGTTGGTTGCATCACAGCAGGTTGGTTCCGGTAGTCGCCTGCATTCTAGGCGGTATCTTAGTGCTATTCGCa        >  2:7020503/1‑93 (MQ=255)
                                                                                   cTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCATTaa    <  2:6929622/93‑1 (MQ=255)
                                                                                   cTCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCATTaa    <  2:6032828/93‑1 (MQ=255)
                                                                                    tCCCTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCATTAAg   >  2:2591270/1‑93 (MQ=255)
                                                                                     cccTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCATTAAGt  <  1:2394394/93‑1 (MQ=255)
                                                                                     cccTTACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCATTAAGt  >  2:494849/1‑93 (MQ=255)
                                                                                         |                                                                                        
TTGGTCATAAGTAATTACCGTCAAGTGCCGATGACTTTCTATCAGGAGTAAACCTGGACGAGAGACAACGGTAATGAATACAACTCCCTCACAGCGATTAGGTTTTTTGCATCACATCAGGTTGGTTCCGTTATTTGCCTGCATTCTAGGCGGTATCTTAGTTCTATTCGCATTAAGT  >  NC_000913/1490421‑1490598

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